,OCEEDIN GS 
AYSICAL | ‘SOCIETY. 


_ Section A 


4 January 1951 


“contents: 


[ Cascade heory ee teageanen bog | Ceo aye 1 
K. H. Banker andl Dy c c. holves. The neler Inter action baneth of the Particles ia 


ey Penetrating ‘Cosmic-Ray Showers. . hie He raeicolae gue 4 

Mr. 16; Rippiror. — nye he Importance of Gs Scattering 4 in ipa tel cacao Bega 10 
Dr P, poe ‘CayaNacu, ‘Mr. J. F. Turner, Mr. D. V. Booxer and Mr. H. os DUNSTER. Aan 05 

Investigation of the Secondary Electron Spectrum, oft FA. 4 13 


W. BREMNER. > A ‘Method for Determining Uranium and T hori tum in Rocks ie the 
Nuclear Photographic PITAL cn eRe We sass 


Dr. P. “ToRKINGTON. The Planar Whimioks of Biliviens ane Ae bees ethylene A 


Mie Critical Analysis of the Potential Function ‘ hae Cais : é \ aa abs 
Mr. Suraj NN. Gupta. On the Elimination of Divergencies from Classical tecerodiyniacoies Fei I0 
Dr. J.C. Warp. On the Renormalization of Quantum Electrodynamics PASC epee) a oe 
Dr. A. PAPAPETROU. _ Equations of Motion in General Relativity ngs FAL agai Wa a Ge oy 
Prof. 'P. B. Moon. Resonant Nuclear Scattering of nein tage A ee: and poetry 

; Experiments a 3 : : 76 


Ee: j. V. JELLEY. petection of Dito hd other Fast Charged Pareioles | in Chemid Radi- 
ation, by the Cerenkov Effect 1 in Distilled Water he ‘ : : : ivicearsyA 


M4 Letters to the, Editor : ee : 


alt JD. Jottey and Dr. F. C. Cuampron. The Angular Distribution of Protons | 

Emitted in some F (a, p) Ne Resonance Reactions. . 88 
“Mr. K. R. WirkInson and Mr. J. Witxs. The Thermal Conductivity of Solid Helium . 89 
- Miss J. M. M@Atister and Mr. D. W. Kream. Factors Involved in the Accuracy and 


Reproducibility of Depth Measurements on Nuclear Research Emulsions . 91 
- Mr. F. K. Gowarp and Mr. J. J. Witxins. Measurement of Gamma-Ray Momenta and 
os the Thresholds of the Photo-disintegrations 12C > 34He and 160 > 44He . ; 93 
Mr. F. K. Gowarp and Mr. J. J. Witxins, The Photo-disintegration of Oxygen into 
Two ®Be Nuclei : 94 


- Mr. F. K. Gowarp, Mr. E. J. “Jones, Mr. H.H. H. WaTSON and Mr. D. hs LEES. Deter- 
mination of Fission and Neutron Yields, and the Average Neutron Energy in the Photo- 


disintegration of Uranium . 95 
‘Mr. P. Goopman, Mr. K. P. NicHotson and Dr. H. D. Ratuceser. The Ionization of 
aap ‘Cosmic-Ray Particles : . 96 
~ Mr. W. H. Dove. and Dr. R. F. Barrow. A New Ultra- Violet Band- System of SiF Anal chs: 
NE ae ‘G. STEPHENSON. Calculation of the Oscillator Strengths for certain Band-Systems of 
Ny and C, : ih y é ‘ F A j A j : ! 5 
A evidys of Books. | SG i b 4 : ; ; . 101 
, - The Institution of Electrical Uae a Prot and Reprint Service . - 2 ir LOD 
a Contents for Section B : t : t : : : : POAT : ea va 
_ Abstracts for Section B.. i y : ‘ j p 3 5 3 . 103 


Laie Price to non-members 10s. net, by post 6d. extra. Annual subscription: £5 5s. 
| Composite subscription for both Sections A and B: £9 9s. 
Published by 


THE PHYSICAL SOCIETY 
1 Lowther Gardens, Prince Consort Road, London 8.W.7 


‘ SS SS TIT MA aA a OE 


THE PROCEEDINGS OF THE PHYSICAL SOCIETY 


PSC A IIT! 


PROCEEDINGS OF THE PHYSICAL SOCIETY 


: The Proceedings is now published monthly in two Sections. 


ADVISORY BOARD 


Chairman: The President of the Physical Society (L. F. Bares, D.Sc., Ph.D., F.R'S:). 


E. N. pa C. ANDRADE, Ph.D., D.Sc., F.R.S. Sir CuarLes DaRWIN, K.B.E., M.C., M.LA., 
Sir Epwarp APPLETON, G.B.E., K.C.B., Se.D., F.R.S. 
D.Sc., F.R.S. N. FreaTuEr, Ph.D., F.R.S. 
P, M.S. Buacxett, M.A., F.R.S. G. I. Fincu, M.B.E., D.Sc., F.R.S. 
Sir Lawrence Brace, O.B.E., M.C., M.A., D. R. Harrrex, M.A., Ph.D., F.R.S. 
ee Se ERS. N. F. Mort, M.A., D.Sc., F.R.S. 
Sir JAMES CuaDwIck, D.Sc., Ph.D., F.R.S. Mi OirEin Phe Dee ERS: 
ee oh ee MA RED. F. E. Simon, C.B.E., M.A., D.Phil., F.R.S. 
if Joun € T. Situ, M.A., F.R.S. 
ia ES al Saar peri Sir Gzorcr THomson, M.A., D.Sc., F.R.S. 


Papers for publication in the Proceedings should be addressed to the Hon. Papers Secretary, 
Dr. H. H. Hopkins, at the Office of the Physical Society, 1 Lowther Gardens, Prince 
Consort Road, London S.W.7. ‘Telephone: KENsington 0048, 0049. 


Detailed Instructions to Authors were included in the February 1948 issue of 
the Proceedings; separate copies can be obtained from the Secretary-Editor. 


inet nae = ated 


BULLETIN ANALYTIQUE 
Publication of the Centre National de la Recherche Scientifique, France 


The Bulletin Analytique is an abstracting journal which appears in three parts, Part A 
covering scientific and technical papers in the mathematical, chemical and physical sciences and their 
applications, Part 2 the biological sciences and Part 3 philosophy. 

The Bulletin, which started on a modest scale in 1940 with an average of 10,000 abstracts per 
part, now averages 35 to 45,000 abstracts per part. The abstracts summarize briefly papers in scientific 
and technical periodicals received in Paris from all ovér the world and cover the majority of the more 
important journals in the world scientific press. The scope of the Bulletin is constantly being enlarged 
to include a wider selection of periodicals. 

The Bulletin thus provides a valuable reference book both for the laboratory and for the individual 
research worker who wishes to keep in touch with advances in subjects bordering on his own. 

A specially interesting feature of the Bulletin is the microfilm service. A microfilm is made of 
each article as it is abstracted and negative microfilm copies or prints from microfilm can be purchased 
from the editors. 

The subscription rates per annum for Great Britain are 4,000 frs. (£4) each for Parts 1 and 2, and 
|- 2,000 frs. (£2) for art 3. Subscriptions can also be taken out to individual sections of the Bulletin 


as follows : frs. 
Pure and Applied Mathematics—Mathematics—Mechanics 550 14/6 
Astronomy—Astrophysics—Geophysics .. ee a 700. = 18/- 
General Physics—Thermodynamics—Heat—Optics—Elec- 
tricity and Magnetism .. ee He ae oe 900 22/6 
Atomic Physics—Structure of Matter ae ai aie 325 8/6 
General Chemistry—Physical Chemistry .. 325 8/6 


Inorganic Chemistry—Organic Chemistry—Applied 


Chemistry—Metallurgy 1,800 45/- 
Engineering Sciences ie me a a oe 1,200 30/- 
Mineralogy—Petrugraphy—-Geology—Palaeontology rh 550 = 14/6 
Biochemistry—Biophysics—Pharmacology ae ee 900 22/6 
Microbiology—Virus and Phages  .. ey a as 600 15/6 
Anima] Biology—Genetics—Plant Biology . . xe a 1,800 45/- 


Agriculture—Nutrition and the Food Industries .. 550 8 14/6 


Subscriptions can be paid directly to the editors : Centre National de la Recherche Scientifique, 
18, rue Pierre-Curie, Paris 5¢me (Compte-chdque-postal 2,500—42, Paris), or through Messrs. H. K. 
Lewis & Co. Ltd., 136, Gower Street, London W.C. 1. 


HANDBOOK 
OF THE 


Pity SICAL. SOCIETY’S 
35.h EXHIBITION 


OF 


) SCIENTIFIC INSTRUMENTS 
| AND APPARATUS 


1951 


5s.; by post 6s. 


To be published at the 
beginning of March 


Orders, with remittances, to 
THE PHYSICAL SOCIETY 


1 Lowther Gardens, Prince Consort Road. 
London S.W.7 


PTICAL GROUP 


THE PROCEEDINGS OF THE PHYSICAL SOCIETY i 


THE HANDBOOK OF THE 
PHYSICAL SOCIETY’S 34th EXHIBITION 
OF SCIENTIFIC INSTRUMENTS 
AND APPARATUS, 

1950 


5s.; by post 6s. 


Orders, with remittances, should be sent to 


THE PHYSICAL SOCIETY 
1 Lowther Gardens, Prince Consort Rd., London S.W.7 


BINDING CASES 
for the 
PROCEEDINGS OF THE 
PHYSICAL SOCIETY 


Binding cases for Sections Aand B (separate) for 
Volume 63 (1950) may be obtained for 7s. each, 
post free, from the Offices of the Society. 


The Proceedings may be bound in the Society’s green 
cloth for 13s. 6d. each. Journals for binding should be 
sent direct to Messrs. Taylor and Francis, Ltd., 
Red Lion Court, Fleet Street, London E.C.4. 
Remittances should be sent to the Physical Society. 


PHYSICAL SOCIETY SPECIALIST GROUPS 


The Physical Society Optical Group exists to foster interest in and development of all branches of 
Dtical science. ‘To this end, among other activities, it holds meetings about five times a year to 
iscuss subjects covering all aspects of the theory and practice of optics, according to the papers offered. 


OLOUR GROUP 

The Physical Society Colour Group exists to provide an opportunity for the very varied types of 
orker engaged on colour problems to meet and to discuss the scientific and technical aspects of their 
ork. Five or six meetings for lectures and discussions are normally held each year, and reprints of 
apers ate circulated to members when available. A certain amount of committee work is undertaken, 
id reports on Defective Colour Vision (1946) and on Colour Terminology (1948) have already been 
iblished. 
OW TEMPERATURE GROUP 

The Low Temperature Group was formed to provide an opportunity for the various groups of 
sople concerned with low temperatures—physicists, chemists, engineers, etc.—to meet and become 
miliar with each other’s problems. The group seeks to encourage investigations in the low temperature 
ld and to assist in the correlation and publication of data. 
COUSTICS GROUP 

The Acoustics Group was formed to meet the long-felt need for a focus of acoustical studies in 
reat Britain. The scope includes the physiological, architectural, psychological, and musical aspects 
acoustics as well as the fundamental physical studies on intensity, transmission and absorption of 
und. The Group achieves its object by holding discussion meetings, by the circulation of reprints 
id by arranging symposia on selected acoustical topics. 


Further information may be obtained from the Offices of the Society : 
1 LowTHER GaRDENS, PRINCE Consort Roap, Lonpon S.W.7. 


iv 


THE PROCEEDINGS OF THE PHYSICAL SOCIETY 


PHILIPS 
X-ray” 


iffraction 
Unit 


A product cf N.V. Philips, 
Glocilampenfabrieken, Eindhoven, Holland. 


This apparatus constitutes the ideal flexible diffraction machine. 
With a power unit rated at 60 kvp. 20 ma. any of the modern 
techniques can be readily employed. Among the many appli- 
cations for this machine are:— 


CHEMICAL IDENTIFICATION - PHASE ANALYSIS - SOLID SOLUBILITY 


STUD A =ssINGLEGRYSTALSSTUDIES PARTICLE SIZE MEASURE- 
MENTS ? METALLOGRAPHIC STUDIES : MOLECULAR WEIGHT 
DETERMINATION - STRUCTURE IDENTIFICATION - STRESS ANALYSIS 

FIBRESSTUDIES.. a -PROCESS CONTROL vax ORIENTATION STUDIES 


PRECISION MADE - GENEROUSLY RATED - ALWAYS READY FOR USE 


Fully descriptive catalogue on request 


PHILIPS ELECTRICAL 


LIMITED 
Manufacturers of : 
INDUSTRIAL X-RAY EQUIPMENT : LAMPS AND LIGHTING EQUIPMENT : SOUND AMPLIFYING 
EQUIPMENT U MAGNETIC FILTERS : BATTERY CHARGERS AND RECTIFIERS ¢ CAPACITORS 
AND MOTRONIC EQUIPMENT : HIGH FREQUENCY GENERATORS 


PHILIPS ELECTRICAL LIMITED, CENTURY HOUSE, SHAFTESBURY AVENUE, LONDON, W.C.2. 


(XB 731B) 


THE PROCEEDINGS OF 
THE PHYSICAL SOCIETY 


Section A 


VoL. 64, Part 1 1 January 1951 No. 373 A 


Cascade Theory including Ionization Loss 


By L. JANOSSY* anp H. MESSEL 
Dublin Institute for Advanced Studies 


MS. received 24th July, 1950 


ABSTRACT. The Bhabha—Chakrabarty solution of the cascade equations with 
ionization loss is discussed, and a simplified method of numerical evaluation suggested. 


§ 1. 

ANY treatments of the cascade theory have been given, both in ‘ approxi- 
\ / | mation A ’, i.e. using asymptotic cross sections and neglecting ionization 
loss, and also using asymptotic cross sections and including effects of 
ionization loss; for a summary of the older literature see for example Janossy 
(1950). Recently two important articles have appeared (Snyder 1949, Bhabha and 
Chakrabarty 1948): the latter is an elaboration of the original paper by Bhabha and 

Chakrabarty (1943). 

The treatment of the ionization term by Bhabha and Chakrabarty has, in our 
opinion, the great advantage that the first order correction of the effects of ionization 
loss is expressed as a shift of the energy spectrum by an amount of the order of the 
critical energy; this correction gives a direct physical insight into the mechanism 
of the cascade. ‘The Bhabha-Chakrabarty treatment was criticized and it was 
alleged that the first approximation did not give a sufficiently accurate description of 
the phenomenon. Inthe meantime it has been shown by Bhabha and Chakrabarty 
(1948) that the latter criticism need not be taken too seriously: in such regions of 
energy where the cross sections can be relied upon, the deviations between their 
own treatment and the mathematically more exact treatment of, for instance, 
Snyder are not of importance. In our own opinion the treatment of Bhabha and 
Chakrabarty gives sufficient accuracy for practical purposes and has the great 
advantage of simplicity. For small depths, however, Bhabha and Chakrabarty 
sacrifice this simplicity and evaluate the average numbers of particles by a 
rather elaborate procedure; the object of the present paper is to show that this 
complicated procedure is unnecessary. 


Se2e 
The first approximation to the average number of electrons of energy greater 
than E along a path zg to which a primary energy E) gives rise can be written 


J yotico ES y-1 dy 
= —————$___ K(y,z) —> Cis ne es 1 
T; (Eo, E, 2) ne (ow) (y, rei ( ) 


* Now at the Institute of Physics, Eétvés I.orant University, Budapest, Hungary. 
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where g(y, 2) is a slowly varying function defined by Bhabha and Chakrabarty. 
B is the rate of ionization loss and 


exp (—a@2). eee (2) 


2 D 
K(y,2)= ay 


For B=0 equation (1) represents the solution in approximation A. Integrals of 
the type (1) can be evaluated by the saddle-point method. 
In some of the treatments of the cascade problem T; is split into two terms 


T; = TS, + T?, © eke enete (3) 


according to the terms of K given in (2) T4, is evaluated by the saddle-point 
method. T?, cannot be evaluated in this way as it possesses no suitable minimum, 
but 7?, is assumed to be small as compared with 71, for sufficiently large values of 
z and is usually neglected. 

It was pointed out by Janossy and Tzu (1946) that it is quite unnecessary to 
split the integral into two parts according to (3); the integral (1) can be evaluated by 
using the unabridged function K(y, z) in (1). The use of the unabridged 
function K(y, z) has the advantage that saddle points are obtained right down 
to 2=0; for z=0 one finds, for E, >, 2(9, 0)—Ovand Ky, 0) —1. 


ay Ay 
—As 


1 


1 (exact) 
T(E, E,0) = x | 3)" ey ee (4) 
2m jJ\E veo Ta = 1-08 (saddle-point method), 
T 


thus for x =0 the error of the saddle-point method amounts to 8%. Further it was. 
shown that for large z, i.e. whenever 


aiz<as2 (at the saddle-point)— | = eee (5) 


the terms arising from the second term of K(y, 2) have a very slight effect only on 
either the position or height and width of the saddle; thus, provided (5) holds, 

Tp Tt Re ee eee (6) 
The latter consideration shows that one is indeed justified in neglecting T?, 
provided (5) is satisfied. Thus for small depths the method of Janossy and Tzu 


gives errors of about 8% while for large depths it gives results essentially in agree- 
ment with the older procedure. 


§ 3. 

The method proposed by Janossy and Tzu was criticized by Chakrabarty (1946) 
and therefore it appears worth while to give some justification. 

From the point of view of saddle-point integration it is certainly justified to use 
the whole function K(y, 2) rather than the first exponential only. In fact the 
function K(y, 2) is analytic in the whole complex y-plane with the exception of 
poles along the negative y-axis; the separate terms of K(y, 2) given in (2) are more 
complicated—they are analytic only in the complex plane slit along the negative 
y-axis. Thus K(y, 2), as it must do, shows in general the analytic behaviour to be 
expected for the Mellin transform of T, while the separate terms have a more 
cumplicated behaviour. This consideration weighs in favour of using the whole 
function K(y, 2). 

No argument against the use of K(y, ) can be deduced from the fact that the 
second term given in (2) cannot be evaluated by a saddle; every function having 
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) aminimum can be split into two components one possessing a minimum and the 
| other not. 
Thus from general grounds there seems to be no justification for the splitting 
of 7; according to (3). 

§ 4. 
We have thus evaluated T 3 (Ep, E, z) from (1) for E=0 (total number of part- 
@ icles), In E,)/B=6, z=0:5—> 20, using the full expression for K(y, z). In course of 
{ computation it is seen immediately that without loss of accuracy it is passinle to 
_ factorize the integrand into a suitable ‘ slowly varying function’ and a ‘ rapidly 
varying function ’, to determine the saddle point by the latter function and replace 
the former by its alter at the saddle thus defined. We have thus taken 


2 D—-a 1 
Ty (Bo E, 2) ={e(¥o2)} wi (14 myee ls oe) an 


«|(Z)" Dat esp (- a) oe ica (7) 
with G ae + {in ; = ae ai WS ee os (8) 
ANY 
mi] ae Spares 
exp (— 4,2) WISE 
s EG o—1)- *+(In—A1) -a,"2h |" parca (9) 


In the dotted curve of the Figure we have given results obtained from (7), (8) and 
(9). The full line on the same graph gives the results of Bhabha and Chakrabarty 


Log ¥,(0,t) 


0 2 4 6 8 10 12 14 16 18 20 
t (cascade units) 


A plot of log Ny (0, t) against the depth ¢ in cascade units for In E)/B = 6. The full curve 
is that obtained by Bhabha and Chakrabarty, the dotted one obtained by us. 
obtained from (1) employing a rather elaborate method of computing the complex 
integral. The agreement between the two curves is excellent and therefore we 
conclude that the parametric representation (7), (8) and (9) gives a satisfactory 
method of evaluating the integral (9). 
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§ 5. 
For most practical purposes a satisfactory approximation of T, can be obtained 
with the aid of tables of 7, and g, using the relation 


Ts (Eo, £, 2) = Ty (Eo, E+ Bg(¥o,2),2%) ve ee (10) 


Detailed tables of T,, using the method of Janossy and Tzu and of the corresponding 
g-values have been computed by the present authors and are reproduced by 
Janossy (1950, Appendix IT, p. 399-407). 
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The Nuclear Interaction Length of the Particles in 
Penetrating Cosmic-Ray Showers 


By K. HO GARKERANDIC] GCSB LER 
The Physical Laboratories, University of Manchester 


Communicated by P. M.S. Blackett; M/S. received 10th July 1950 


ABSTRACT. ‘The work of Rochester and Butler, published in 1948, and that of Butler, 
Rosser and Barker, published in 1950, has been continued. In particular, the interaction 
properties and the spectrum of the penetrating particles are considered. The interaction 


length of the shower particles with momenta greater than 10° ev/c. is found to be 18 cm. 
lead (+ 20%). 


§1. INTRODUCTION 


HE work of Rochester and Butler (1948 a, b) and of Butler, Rosser and 
Barker (1950)* has been extended using the same apparatus. The cloud 
chamber was operated in a field of 7,500 gauss ; the Geiger counter arrange- 

ment was a seven-fold coincidence set similar to the P-set in the arrangement B 
of B.R.B. Seventy-five penetrating showers were observed under 25 cm. of 
lead and twenty-two under 25 cm. of paraffin wax. These sea-level observations 
are now complete and the total data obtained during all the experiments have been 
analysed. A total of 164 penetrating showers has been observed which shows 240 
particles penetrating the 3-4 cm. lead plate in the chamber. In order to obtain the 
interaction length of the secondary, or shower, particles it is necessary to measure 


* Subsequently referred to as B.R.B. 
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| the momentum spectrum of these particles. The total data for this differential 


momentum spectrum are given in § 2 and the interaction length is discussed in § 3. 
_ Four unusual showers obtained under 25 cm. of paraffin wax are shown in the Plate 
_ and discussed in the Description of Plate (see end of issue). 


§2. THE DIFFERENTIAL MOMENTUM SPECTRUM OF THE 
PENETRATING PARTICLES 

The experiments with either 25 cm. of lead or 25 cm. of paraffin wax as transi- 
tion material provide considerable additional data for the momentum spectrum of 
the penetrating particles. ‘That the counter selection provided a strong bias towards. 
high-energy events is evident from the fact that 80°% of the photographs obtained 
under 25cm. of lead showed electrons, because Fretter (1949) shows that, in 
general, only very energetic interactions give rise to an appreciable electronic 
component. The data obtained in all the sea-level experiments are collected in 
Table 1; no correction has been made for ionization loss in the transition material. 


dablewt 
Momentum range ( x 108 ev/c.) 1-5 5-10 10-20 20-30 > 30: 
No. of particles under lead 16 34 50 21 14 
No. of particles under paraffin wax 0 6 5 Z 7 
sbableyZ 
Momentum range ( x 108 ev/c.) 5-10 10-20 =) 
Cloud chamber (under lead) 3:5+0:-7 2:9+0:7 2:°5+0°8 
Photographic plate 8:8+1:-9 3:-4+0-9 — 


A positive excess of 3:0 + 0-6 is found among the penetrating particles under 
lead and a positive excess of 2-0 + 1-0 under paraffin wax. It is generally accepted 
that a considerable number of the shower particles are 7-mesons. On the assump- 
tion that the showers are produced equally by both neutral and charged primaries. 
(cf. Janossy and Rochester 1943 and Walker 1950), the number of positive and 
negative mesons produced should be about equal, and so the observed positive 
excess must be attributed to the presence of protons. Identifiable protons are 
strongly absorbed in the 3-4cm. lead plate ; consequently very few identified 
protons appear in the spectrum of Table 1. 

The investigations of Fowler (1950) and Camerini et al. (1950), using the 
photographic-plate technique, suggest that 80% of the shower particles 
(1/Imin<1:5) are 7-mesons. Thus it is interesting to compare the cloud-chamber 
and photographic-plate results for the positive excess among the secondary 
particles. 

In Table 2 the discrepancy between the cloud-chamber and photographic- 
plate data in the lowest momentum range is due to the strong absorption of protons. 
in the transition material and lead plate of the cloud chamber. For this reason the 
percentage of 7-mesons among the shower particles in the chamber cannot be 
determined and, therefore, asimple comparison with the photographic-plate data 
cannot be made. Between 10%ev/c. and 2 x 10% ev/c. agreement is good, but no 
comparison is possible in the highest momentum range although the cloud-chamber 
data provide some evidence for an appreciable component of protons among the 
very energetic shower particles. Individual showers of very high energy, photo- 
graphed in the cloud chamber, show a very large positive excess, so it seems certain 
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that a process must exist whereby several protons, and presumably neutrons, can 
be accelerated to relativistic velocities. 

Beyond momenta of 10° ev/c. the observed differential momentum spectrum of 
the penetrating particles from lead is of the form 


N(E) dE o E-@5+03dR, 


and is thus similar in form to the spectrum of mesons at formation. 

The data obtained using paraffin wax as transition material are very meagre, but 
relatively more high energy particles are found under paraffin wax than under lead. 
This is to be expected, as fewer secondary interactions within the nucleus, made by 
energetic recoil nucleons from the first collision, are likely to take place in paraffin 
wax than in lead. 

A number of primary and secondary interactions, mostly of high energy, have 
been observed in the lead plate. Some of the shower particles produced in the 
plate are suitable for measurement, and the differential momentum spectrum of 
44 particles is given in Table 3. 


Tabless 
Momentum range ( x 108 ev/c.) 1-5 5-10 10-15 15-20 >20 
No. of particles 19 9 9 5 1 


The overall positive excess of these particles is 3:1+1-2 and the maximum 
intensity occurs at a much lower momentum than is the case for the penetrating 
particles. This is to be expected, since the selection is made from a much wider 
angular distribution of particles than was possible in the results of Table 1, and 
also the shower particles are not required to penetrate more than about 1 cm. of 


lead. 


§3, THE NUCVEAR INTERACTIONS OF THE SHOWER PARTICLES 


The nuclear interactions of the shower particles observed in the lead plate fall 
into three main categories: (1) lightly ionizing particles undergoing a large-angle 
scattering considerably in excess of the calculated coulomb scattering appropriate 
to the particle momentum and the thickness of the lead plate; (ii) interactions of 
lightly ionizing or neutral particles with the production of both heavily and lightly 
ionizing secondaries; (ii1) lightly ionizing particles other than electrons stopping 
in the lead plate. 

Thirty-eight secondary interactions have been observed, consisting of 8 of type 
(i), 23 of type (ii), of which 14 have ionizing primaries, 5 neutral primaries and 
4 either ionizing or neutral primaries, and 7 of type (ili). 

The interaction length of the particles constituting penetrating showers can be 
deduced from the observed numbers of interactions and penetrating particles. 
B.R.B. found a value of 35 cm. of lead for this interaction length. Using all the 
data now available, the value is 26cm. of lead. It is necessary, however, to assess 
the efficiency of detection of the three types of nuclear interactions in order to 
correct this figure. 


(a) Events of type (t). Anomalous scattering. 

The majority of the eight particles showing large-angle scattering in the lead 
plate cannot be protons because two emerge from the plate with momenta less than 
6 x 108 ev/c. and with minimum ionization, and four of the remaining six particles 


The Interaction Length of Shower Particles 7 


are negative. ‘The variation of the scattering angle with the particle momentum is 
shown in the Figure; the curve corresponds to five times the probable multiple- 
coulomb scattering calculated using Williams’ formula. 

The eight examples of anomalous scattering may represent the nuclear 
scattering of 7-mesons. On this assumption, and on the assumption that half of 
the shower particles are 7-mesons, the cross section for this processis 2-8 x 10-2? cm? 
per nucleon ( + 50°) or approximately 50cm. of lead. This is considerably larger 
than the cross section found by Fretter (1949) and by Brown and McKay (1950), 
but is in good agreement with that calculated by Heitler on the basis of damping 
theory. ‘The present very limited data do not permit an examination of the 
variation of the cross section with energy. The value of the interaction length of 
all the shower particles, i.e. protons and 7-mesons, indicates that 7-mesons must 
initiate interactions of types (ii) and (iii) as well as of type (1). 
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Events of type (i) are detected easily if the scattering is in the plane of the 
photograph, but there is considerable uncertainty when the scatter is towards or 
away from the camera. In addition, the geometrical form of the lead plate 
introduces a bias against very large-angle scattering, and if the scattering is small 
there is confusion with multiple-coulomb scattering. The probability of 
observing this type of interaction cannot be calculated, and we have to make the 
assumption that no events of this type are omitted from the analysis. 


(b) Events of type (i). Nuclear interactions in the lead plate. 

Events of type (ii) will be detected with a high probability if lightly ionizing 
particles are emitted, but only with a considerably lower probability if such par- 
ticles are not present, due to the fact that, for moderate excitation energies, the 
majority of the energy is carried by particles of intermediate ionization which are 
strongly absorbed in lead. In this connection the data of Brown et al. (1949) are 
useful; the average number and energy of the grey tracks, associated with ‘ stars’ 
of various total energies occurring in silver and bromine nuclei, are estimated. If 
we assume that these data apply to stars produced in lead nuclei, a rough calcu- 
lation suggests that the probability of detecting an interaction with an excitation 
energy of 500 mev. is about one half in the present investigations. Below this 
energy the probability of detection falls quite rapidly, whilst above 500 Mev. 
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it rises towards unity at the energy where particle production is a frequent process, 
which is likely to be in the range (1-2) x10%ev. Thus, for the purpose of 
estimating an interaction length, we shall confine our attention to the interacting 
particles with momenta greater than 10%ev/c., corresponding to a proton 
energy of about 600Mev., and assume that all the observed interactions of 
type (ii) are due to such particles, and that the probability of detecting such inter- 
actions is unity. 


(c) Events of type (ui). Fast particles stopped in the plate. 


All the particles which produced interactions of type (ii) have measured 
momenta greater than 10%ev/c. Electrons with similar momenta would have a 
high probability of producing electron cascades in the plate. The stopping 
particles were observed immediately below 5 cm. of lead and, since they are not 
accompanied by many soft electrons, are definitely not parts of cascade showers. 

If the stopping particle has a measured momentum less than 10% ev/c. then 
in many cases it is less easy to be sure that itis not anelectron. In fact we have not 
found any such stopping particles which cannot be electrons. Events of this 
kind must occur, for example, due to the presence of lightly ionizing protons with 
kinetic energy less than 500 Mey. If such a proton interacts in the plate, then a 
simple star will be formed and in most cases all the heavy tracks will be absorbed 
inthe lead. ‘Thus the star will be invisible and the event should be classified as. 
type (iil). 

Adequate correction can be made for unobserved events of type (iii) by con- 
fining attention to energetic stopping particles, that is, with momenta greater than 
10° ev/c. 


(d) The calculation of the interaction length of the shower particles. 


To calculate the interaction length of the ionizing shower particles we neglect 
all events with incident momenta less than 10° ev/c. and determine the thickness of 
lead penetrated by the particles with momenta greater than 10%ev/c. From Table 1 
we find that 64° of the penetrating particles have momenta greater than 10° ev/c., 
so that of the 240 penetrating particles, 153 of at least this momentum penetrate the 
3:4cm. lead plate, accompanied by 27 interactions of ionizing particles. The 
interaction length is thus 18cm. of lead(+20%), corresponding to 1:3 times 
the geometrical cross section in lead. ‘lhe error quoted is a purely statistical one 
based on the observed number of interactions. ‘There is a further small uncer- 
tainty due to undetected nuclear interactions (mainly scatterings) of particles with 
momenta greater than 10%ev/c. ‘Therefore it appears likely that the interaction 
length calculated on the foregoing assumptions will tend to be too great. 

The positive excess of the shower particles with momenta above 10° ev/c. is 
2:7+0-7. Assuming that this positive excess is due to the presence of protons 
and that the proton interaction length is 14cm. of lead, then the overall inter- 
action length of 18cm. of lead indicates that the high energy 7-mesons possess 
an interaction length of the same order as that of fast protons. 

Camerini et al. (1950) have made a determination of the interaction length of 
shower particles (//J,,;,<1:5) in photographic emulsion and find a value 1-1 times 
the geometrical cross section. 

Several workers have investigated the interactions of secondary penetrating 
shower particles in cloud chambers containing several thin lead plates (Fretter 
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1949, Lovati et.al. 1950, Brown and McKay 1950). The values obtained for the 
‘interaction length, after applying corrections for invisible interactions, are all 
: about the same and equal to 30 cm. of lead. This value is considerably larger than 
that calculated above, and a possible explanation is that a number of interactions, 
due to particles with energies of several hundred Mev., with the emission of heavily 


ionizing secondaries only, were not effectively corrected for by the procedure 
adopted. 


§4. CONCLUSIONS 


(1) The shower particles with momenta greater than 10®ev/c. must consist 
mainly of 7-mesons and protons. ‘The number of protons is likely to be between 
25 and 50% of all the shower particles. 


(2) The interaction length of the shower particles with momenta greater than 
10° ev/c. is found to be 18cm. lead (+ 20%). 
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ABSTRACT. Inelectron and proton synchrotrons there exists a minimum useful vertical 
aperture of the vacuum chamber, and for apertures less than this the beam intensity falls 
to zero. It is considered, on the basis of a comparison of theory and experiment, that this 
phenomenon is in the main due to free oscillations which arise from gas scattering. The 
final beam size in proton synchrotrons may be substantial for this reason also, but for 


-cyclotrons the large spread which is known to exist must be attributed to some other 
phenomenon. 


and Courant (1948, 1949). The effect of this scattering is to set up vertical 
and horizontal free oscillations, whose mean square amplitude rises 
sharply to a maximum when the particles have reached four times their injection 


-energy, and then falls asymptotically to zero. It may be shown that this maximum 
value is, in the vertical direction, 


B?=8-17 x10 


T= theory of gas scattering in a synchrotron has been given by Blachman 


2 Ree 

VT Sanaee (1) 
with R the orbit radius in centimetres, P the air pressure in mm. Hg, and Vand T; 
the energy gained per revolution and injection energy respectively, both in electron 
volts. Equation (1) gives the obvious result that the vacuum should be as good as 
possible, and injection energy high. It is no doubt a good reason why betatron 
output increases so rapidly with injection voltage. 

If (B12. is of the order of the vertical aperture of the vacuum box, many 
particles will hit the walls of the doughnut and so be lost. Nowit has been reported 
by Elder, Langmuir and Pollock (1948) and also by Lawson (1948) that an investi- 
gation of the variation of synchrotron (or betatron) beam intensity, as a function of 
the vertical aperture of the accelerating space, shows the existence of a minimal 
useful aperture of the order of 1:0 cm. in the centre of the doughnut; for apertures 
less than this the output from the machine falls to zero. McMillan (unpublished 
lecture) also has reported that a similar effect occurs in the Berkeley 300 mev. 
electron synchrotron. Further, Lasich and Riddiford (1947) have noticed that for 
the Melbourne betatron (now a 14 Mev. synchrotron), pressures less than 
2x10-° mm. Hg are necessary for efficient operation. This point was also 
stressed by Kerst (1948) (see also Fry et al. 1948). 

So far no reasonable explanation of this minimal useful aperture has been given 
although the answer has been sought in terms of resonances between the various 
particle oscillations (Lawson 1948). However, the reason seems clear if we 
substitute some figures into equation (1). The results for several accelerators 
are summarized in Table 1 for an operating air pressure of 1 x 10->mm. Hg. 

It should be noted that, assuming no loss of particles due to wall collisions, the 


distribution about b=B,,,,/1/2 is a Rayleigh distribution, i.e. the probability 
that a particle has a value of B between B and B+ dB is 


(6?)"1Bexp{—(B?/2b3)}dB. ae, (2) 
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The integral of this distribution has a shape similar to the curves presented by 
Elder et al. and by Lawson. But equation (2) is not correct for the case in which 
particles are being continuously lost to the walls. _ If one calculates the fraction of 
particles P not scattered to the walls as a function of the available vertical semi- 
aperture from Blachman and Courant’s theory, it is found that the curve has the 
‘same general shape as the integral of the Rayleigh distribution of equation (2), but 
that P is practically zero for values of B less than about B,,,, as defined by 
‘equation (1). This accounts for the minimal useful aperture. 

A minimal aperture due to magnetic field inhomogeneities is not possible, and 
-since the particles move in phase with the field disturbances, the curve showing the 
variation in beam intensity with aperture need not be symmetrical about the 
‘median plane, as it is found to be in practice. Further, with the type of gun used 
in electron accelerators, a wide spread of injection oscillations about the ideal orbit 
is to be expected, so that no minimum useful aperture would occur due to this type 


Vablesi 
2B Minimal 
Accelerator IRS). We T; maxX Aperture 

ky oh) 
Melbourne 3 Mev. betatron 75 20 I SOY DED — 
po ey 300 Mev. synchrotrons 100 1500 80x10 1:7 = 
‘Cornell 
Malvern 30 Mev. synchrotron 10 DS 1-0 x 104 iil 0-6 
GEC. 70 Mev. synchrotron 29 200 S305) SC OE isi] 1-0 
Oxford 16 Mev. synchrotron 20 20 2°5 x 104 D8 = 


* During betatron acceleration. 


-of free oscillation alone. Of course both of these phenomena occur in practice, 
and an investigation of the manner in which the intensity—aperture curve depends 
-on the azimuth would permit of separation of the phenomena, in any case in which 
they are of comparable magnitude to the gas-scattering oscillations. Such an 
investigation will be an invaluable aid to efficient operation of the proton 
‘synchrotron. 

It is worth while making two further points which have been investigated. 
‘The Blachman—Courant calculations, based on current design figures, indicate 
that proton synchrotron acceleration is quite practicable at physically realizable 
pressures. ‘The operation of a model at Berkeley substantiates this conclusion 
(Sewell et al. 1950). However, gas scattering appreciably affects the final size of 
the beam, and may consequently make efficient extraction of the beam quite 
difficult. This is because the effects of damping of the free oscillations 
are being continuously annulled by further gas scattering throughout the acceler- 
ationcycle. For the Birmingham synchrotron (Oliphant, Gooden and Hide 1947) 
‘injecting at 500 kev. on a radius of 450.cm., with V = 200 volts, it is estimated that 
the final beam size due to free gas-scattering oscillations will be 4-4 cm. horizontally 
and 2:5 cm. vertically, if the particles fill the vertical aperture of the vacuum box 
when 7,=4T,, which they do unless the air pressure is less than 3-5 x 10-6 mm. Hg 
(see Table 2). Considering phase oscillations only, the final beam width is 
‘0-9cm. This does not, however, mean that the severe restrictions upon the amount 
that the radio frequency can slip into error in a time of the order of one phase- 
oscillation period can be relaxed, for although most of the particles are finally 
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geometrically outside the region of stable phase, they are still within that region 
when it is represented as a velocity-phase diagram. ‘Table 2 shows the air 
pressure p (in mm. Hg) at which the vertical semi-aperture a of the vacuum box is 


equal to the value of (B?)}/2. from equation (1), for the three proton synchrotrons at 
present under construction. The advantages of high injection energy and rapid 
acceleration are apparent. ‘The Table also includes figures for the model 
accelerator. In fact it is found that, for an operating pressure of 2 x 10- mm. Hg, 
the beam intensity in this machine is very small when the vertical aperture is. 
reduced to 10cm. by insertion of screens. 


Table 2 
Laboratory Energy (10° ev.) R (cm.) V T;(Mev.) a(cm.) p (mm.Hg) 
Birmingham 1:3 450 200 0:5 5-0 3255c 10 
Brookhaven 3-0 1000 1040 3:0 8-9 BOP alm 
Berkeley 6:5 1525 1750 10-0 14-6 ESP GhOme 
15-0 Ng Ses (=P 
Berkeley (model) 0-006 350 45 0-6 4 5-0 1-9 10-8 


Now a recent paper from Berkeley (Henrich et al. 1949) has stated that for 
200 mev. deuterons accelerated in the 184 in. synchro-cyclotron, the final beam 1s. 
about 10cm. wide and 3cm. high. Further, only 2% of the beam has been 
successfully extracted (Powell et al. 1948). Since the phase oscillation amplitude 
at full energy is only 4-4 cm., it remains to explain the presence of free oscillations of 
horizontal amplitude 2-8 cm., and vertical amplitude 1-5cm. In view of the above 
considerations, and the fact that the ions are formed at thermal energies, it seemed 
possible that gas scattering was the cause. However, an investigation based on 
adapting the Blachman—Courant theory to synchro-cyclotrons shows this not to 
bethecase. But initial oscillations of the order of several centimetres arising from 
gas scattering will explain at least part of the factor of ten increase in beam current 
obtainable by shaping the central region of magnetic field (Henrich et al. 1949). 
It is possible to convince oneself that the final beam spread is not due to injection 
oscillations and certainly the cause is not magnetic field inhomogeneities. 
Possibly the vertical oscillations arise from electrostatic defocusing although this 
is unlikely since in a synchronous accelerator most of the phase stable particles 
lie on the falling part of the wave, in which case the electrostatic action is a 
focusing one. A solution for the horizontal oscillations particularly is most 
desirable from the point of view of accelerators of the cyclo-synchrotron type, 
as proposed by Oliphant (1950). 
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ABSTRACT. ‘The secondary electron spectrum of 1*8%Au has been studied using very 
strong sources, up to 0-3 curie, in a magnetic lens spectrometer. In addition to the 
well-known y-ray at 0:411 Mev., evidence of two much less intense higher energy y-rays 
at 0-67 Mev. and 1:09 mev. has been found, both from the Compton and photoelectric 
‘spectra. The intensities of these two y-rays with respect to the one at 0:411 mev. have 
been measured. ‘Two photoelectric peaks at low energy, in about the region expected 
for Ka and Kf x-radiation from elements in the neighbourhood of gold, have been 
examined with radiators of different atomic number. The radiation responsible for these 
peaks has been subjected to critical absorption measurements in the spectrometer. Its 
behaviour is compatible with the Ka and Kf x-radiation of gold, and no evidence of 
platinum x-radiation, such as would suggest K-capture, has been found. 


§1. INTRODUCTION 


HE 2-7-day half-life activity of 1**Au has been studied very extensively in 

recent years. It is an ideal nucleus for study because it can be obtained 

easily, free from all other activities, and in a form of high specific activity. 
It is formed by the reaction 1*’Au (n,y) *8Au, and not only is 1®7Au the sole 
stable isotope of gold, but also it has a large cross section for slow neutrons, of the 
order of 100 barns. Other possible reactions, such as (n, p), (n,2p), etc., appear 
to have very small cross sections by comparison. It is surprising therefore that 
a larger measure of agreement has not been obtained by the various investigators. 
There is complete unanimity as to the presence of a y-ray of energy 0-411 Mev., 
as evidenced both by the internal conversion lines and by the secondary 
electron spectrum, this being in coincidence with f-rays of maximum energy 
0-960 + 0-005 Mev. However, Levy and Greuling (1948), using a semi-circular 
focusing instrument, find that internal conversion lines appear, corresponding 
to y-energies of 0-158 and 0-208 mev., with an intensity of about 15°% of the main 
y-ray. In addition, what is believed to be an x-ray at 70 kev. has been observed 
and investigated by various workers, notably by Feather and Dainty (1944). 
Very few workers have undertaken an accurate measurement of the secondary 
electron spectrum; the most recent is that of Siegbahn and Hedgran (1949) 
who, using a double-focusing spectrometer, find evidence only for the 0-411 kev. 
y-ray, and conclude that if y-rays at 0-158 and 0-208 Mev. exist, they must be 
present to less than 2°% of the main branch. 

The present investigation of the secondary electron spectrum of '8Au has 
been carried out using a short magnetic lens spectrometer, essentially similar 
in design to that described by Deutsch, Elliott and Evans (1944). Investigations 
have also been made on the primary 6-spectrum and on fy and y-y coincidences. 
These will be reported at a later date. 
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§2. THE SPECTROMETER 
The magnetic lens was used with a focal length of 25 cm. and unit magnification,. 
and is able, with the present design of coil, to focus electrons with a maximum. 
energy of 5 Mev. Current for the coil is supplied by a5 kw. amplidyne generator, 
and is stabilized by conventional methods to a few parts in ten thousand. For y-ray 
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Figure 1(6). Radiator arrangement for critical absorption experiment. 


work heavy lead bafHes are used to reduce the effect of direct y-radiation at the 
counter, and also to prevent the scattering of secondary electrons into the counter 
window. ‘The vacuum chamber is made of aluminium, and the baffles are 
aluminium clad also, to reduce electron scattering. The y-ray source is placed in 
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a small cylindrical aluminium box, which is pushed into a blanked-off aluminium 
tube lying on the axis of the spectrometer. ‘The tube enters the vacuum chamber 
via a Wilson seal, so that the axial position of the source may be varied. A three-- 
position radiator holder is mounted before the source, and operated from outside 
the vacuum through a Wilson seal (Figure 1(a)). The radiators, which varied 
in thickness from 1 to 58 mg/cm? of lead and from 1 to 2 mg/cm? of tin and gold, 
were, with the exception of the heaviest lead radiator, made by evaporating layers 
of these metals on to 5 mg/cm? mica discs. The heavy lead radiator was of foil 
rolled to the required thickness. 

For work at higher energies, the 6mm. diameter counter window was of 
2 mg/cm? mica, while for the investigation of the x-rays 0-5 mg/cm? nylon windows 
wereused. ‘The spectrometer was used with a transmission of 0:2 °% anda resolving 
power of about 3° at lower energies. At higher energies, where more than the 
outer radial section of the coil was required to focus the electrons, the resolution 
was somewhat worse. 


$3. THE SOURCE 
The gold used was in the form of foil 0:13 mm. thick, and discs of 5mm. 


diameter were punched from it. The gold was subjected to spectrographic and 
chemical analysis before irradiation, and was found to contain not more than a few 
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Figure 2. Lifetime measurements on gold. 


parts per million of any impurity. No chemical separation was made after irradia-- 
tion, and there remains the possibility that some other nuclear reaction may have 
occurred. However, the decay of an irradiated sample was followed over seven 
half-periods and compared with the uranium standard in a standard jig. The 
decay was found to be simple, and with half-period, obtained by least squares 
fitting, of 63-86+0-13 hours (Figure 2). In the initial investigations sources of 
the order of 10 millicuries were used. These were produced by the (n, y) reaction 
in the smaller Harwell pile. More detailed work at high energies required source 
strengths of up to 300 millicuries. For work at low energies we were able to obtain 
sources with an order of magnitude of greater specific activity from the larger 
Harwell pile. 
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§4. THE EXPERIMENTAL MEASUREMENTS AT HIGH ENERGIES 

The entire secondary electron spectrum of !%Au has been covered several 
times, with varying conditions, with experimental points spaced at the most at Pied. 
intervals of momentum on a logarithmic scale, and in the regions of interest, at 
considerably closer spacing. All points have been taken to a statistical accuracy 
of within 1-5°%, and most to1°%. During the investigations a total of well over ten 
million counts was recorded. An initial run over the spectrum of secondary 
electrons ejected from aluminium gave a main Compton spectrum which appeared 
to be due to a single y-ray of energy about 0-41 Mev. ‘There was also evidence 
of photo-electrons ejected from aluminium by the same y-ray. However, for 
higher values of the momentum, the counting rate did not at once fall to the back- 
ground value, but appeared to drop down in two steps. ‘The first was partly 
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Secondary electron spectrum from 1°8Au. 

obscured by the photoelectric line due to the 0-411 Mev. y-ray, and was at 0-49 mey. 
The second was at 0:87 Mev. ‘These results were confirmed when stronger 
sources, up to 300 millicuries, were used. A polystyrene radiator was used in 
place of the aluminium one, and with this no photoelectric line was observed. 
It was evident that the two steps represented the edges of Compton spectra due 
to y-rays of energies (given by the Compton edges) of 0-67 and 1-09 mev. (Figure 3). 
Concurrent measurements were made on the photoelectrons ejected from lead 
radiators of thickness varying from 1 to 58 mg/cm?. The K and LM photoelectric 
lines due to the y-rays of 0-411 and 0-67 Mey. were clearly resolved, and those due 
to the 1:09 Mev. partially resolved. In the measurements, alternate readings with 
and without the lead radiator were made by rotating the radiator holder. The decay 
of the higher energy y-rays was followed over nearly four half-periods, check runs 
being made with the 0-411 Mey. y-rays. The half-periods, as determined by 
a least squares fit to the data, were, for the 0:67 and 1-09 Mev. y-rays respectively, 
63-4 + 1-2 hours and 65:8 + 1-6 hours (Figure 2). It is therefore extremely likely 
that they are y-rays originating from the decay of Au. 
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§5. THE DETERMINATION OF THE ENERGY OF THE y-RAYS 


To obtain sufficient intensity in the photoelectric peaks due to the higher 
energy y-rays thick radiators had to be employed, up to58 mg/cm?. Experiments 
were made with the 0-411 Mev. y-ray with radiators of various thicknesses. For 
thin radiators the spectral lines were symmetrical, and of width given by the 
expected resolving power. As the thickness of radiator was increased, the line 
increased in height and developed a tail on the low energy side. However, the 
high energy side remained practically unchanged. For still greater thicknesses of 
radiator the increase in height was less marked, as also was any change in the shape 
of the line. Finally, for a thickness of radiator about equal to the range of photo- 
electrons no further change in height or line shape occurred. This we refer to as 
the saturation line shape, and experimentally the peak of the line, under these 
circumstances, was found to move down from the posicion of the peak for a thin 
radiator by only 1% in energy. Moreover the saturation line shape was very 
nearly the same for all three y-rays, and a correction of 1 % in energy was therefore 
made for the higher energy y-rays. ‘The spectrometer was calibrated by replacing 
the lead radiator by a source of Th (B+C) of the same diameter, formed by the 
standard technique, and determining the position of the F-line, the energy of which 
has been accurately determined by Ellis. A correction was made for the effect 
of the horizontal component of the earth’s field. The y-ray energies were 
determined by adding the K-binding energy of lead, 88-2kev., to the electron 
energies so determined. The values were 0-410 + 0-006 Mev., 0-671 + 0-009 mev., 
1-092 +0-013 mev. ‘The energies were also obtained, less accurately, from a 
consideration of the Compton edges, correction being made for the finite resolution 
of the instrument, and the results were 0-412 + 0-007 Mev., 0-668 + 0-015 mev., 
1-075 + 0-025 Mev. 


§6. THE DETERMINATION OF THE RELATIVE INTENSITIES 
OF THE y-RAYS 

Estimates of the relative intensities of the y-rays have been made both from the 
relative numbers of electrons in the Compton spectra and from the heights of the 
photoelectric lines, with good agreement. It is assumed that the relative number 
of Compton electrons of all energies, per y-ray, as a function of the y-ray energy 
will be given closely by the relative quantum efficiencies of a Geiger counter having 
a cathode of the same material as the radiator, provided that the Compton effect is 
very much greater than the photoelectric effect, as it is in this case. The relative 
efficiencies of a brass cathode counter at various quantum energies have been 
determined by Bradt et al. (1946), and the ratio of the efficiencies of counters with 
cathodes of various materials is given by the formula of van Droste (1936). From 
these we may deduce that the relative efficiencies of an aluminium-walled counter 
for quantum energies of 0-41, 0-67, and 1-09 Mev. are in the ratios 1: 2:1: 3-7. In 
order to determine the relative numbers of electrons in the three Compton spectra, 
those corresponding to the higher energy y-rays must be extrapolated back over 
the greater part of their range to zero energy. ‘This extrapolation has been made 
on an energy similarity basis for Compton spectra, following the work of Bleuler 
and Ziinti, but the results are not critically dependent on the actual form of extra- 
polation used. On this basis the ratios of the intensities of the three y-rays are 
1: 0-015: 0-004. The measurements were repeated with thinner sources, down 
to 250 mg/cm?, but the results obtained were within a few per cent of those with 
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the thicker source, indicating that the effect of differential absorption of the y-rays 
within the source material is negligible with our geometry. 

Our measurements of the heights of the photoelectric lines from the three 
y-rays may be interpreted on the basis of an expression for the variation of line 
height with radiator thickness given by Deutsch, Elliott and Evans. This gives 
ratios of intensities of 1: 0-012: 0-004. 

No evidence was found for y-rays of 0-2 and 0-3 Mev. reported by other workers, 
and we can set an upper limit on their intensities, if present, of 2% and 4%, 
respectively, in agreement with Siegbahn and Hedgran. 


§7. THE INVESTIGATION OF THE LOW ENERGY 
SECONDARY ELECTRONS 


In the initial investigation of the secondary electron spectrum, using a source of 
superficial density 10 gm/cm?, a 1-0 mg/cm? lead radiator, and a counter window of 
2 mg/cm?, two photoelectric peaks at energies of 52 and 63-5 kev. were observed. 
They were of nearly equal heights, 60-70°%, of the Compton background. 
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Figure 4. External conversion of low energy radiation from 1°Au in lead, gold and tin, 


Measurements in this region were continued using sources of much higher 
specific activity, obtained from the larger Harwell pile, and of superficial 
density 1-25 gm/cm®, together with a counter window of 0:5 mg/cm? nylon. 
Again peaks were observed, using the same radiator, at the same energies, 
of approximately equal intensities, but now almost equal in height to the back- 
ground. ‘The separation of the lines is approximately equal to the LM separa- 
tion in the gold-lead region. It was first necessary to establish whether 
conversion took place in the K or L shells, especially in view of the conversion 
line corresponding to a y-ray of energy 0:158mev., reported by Levy and 
Greuling. ‘The K-conversion of this y-ray would place it at 70 kev., just above 
the position of the second of the two peaks observed here. The measurements 
were repeated using a 1-0 mg/cm? gold radiator, and only a very small shift, of 
the order of 1 kev., was observed for both lines (Figure 4). The second line 
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appeared to be somewhat narrower with the gold than with the lead radiator. 
Had either of the two lines been converted in the K shell the peak would 
have moved to a higher energy by 8 kev. in changing from the lead to the gold 
radiator. ‘The readings were carried on sufficiently far past the second peak to 
ensure that no larger part of it could be ascribed to the K-conversion of a y-ray at 
0:158mev. ‘There was some evidence in both spectra of a third small peak lying 
above the second, and separated from it by about the same energy difference as 
between the first and second. This third peak, at 74 +2kev., is rather above the 
position expected for the K-converted 0-158 Mev. y-ray and, since its presence is. 
not nearly so well marked in the case of the gold radiator, there is perhaps some 
possibility that it may have moved in the manner expected from K-conversion. 
However, it is much more likely that it corresponds to the M-conversion of a 
quantum radiation, the L-conversion peak of which forms part of the second 
spectral line. 

It is evident that the first photoelectric peak is due to the L-conversion of 
quantum radiation of energy 66:5 + 1 kev., while the second could be due to the 
M-conversion of the same radiation. In view of the presence of the third low 
intensity peak, and also the height of the second peak in relation to the first, it 
seems likely that the second peak is composite, consisting, in addition, of the 
L-conversion of asecond quantum radiation of energy 78-0 + 1 kev., while the third 
peak is due to the M-conversion of the same radiation. A preliminary study of 
the 6-spectrum did not show any internal conversion peaks at low energies, which 
would be expected if low energy y-rays were emitted. On the other hand the 
energies of the two quantum radiations indicated by the secondary electron 
spectrum are compatible with their being Kx and Kf x-rays of an element in the 
neighbourhood of gold. ‘These x-rays would not be expected to give internal 
conversion in the 6-spectrum, even if they were emitted as a result of the internal 
conversion of one of the higher energy y-rays, because the Auger effect is small for 
high atomicnumbers. Ifthe x-rays were emitted as a result of external conversion, 
there is even less likelihood of their giving rise to conversion peaks in the 
B-spectrum. 

The emission of X-rays can arise in a number of ways. It is known that the 
0-411 Mey. y-ray is internally converted to the extent of 2-5-3°%, and this will give 
rise to the same fraction per disintegration of mercury x-rays. ‘The source is 
comparatively thick for the 0-411 Mev. y-rays, and a calculation of the absorption of 
this y-ray within the source, based on the model of a point source embedded in the 
middle of a disc of absorbing material, and using the photoelectric cross section 
given by Heitler (1936), gives a value of 15 °% of gold x-rays per disintegration. ‘This 
is unlikely to differ by more than 30° from the result for a distributed source. 
Gold x-rays will therefore be present to about six times the intensity of mercury 
x-rays. ‘There remains the possibility that 1*Au disintegrates by an alternative 
mode of decay, viz. K-capture to !°8Pt, in which case the emission of platinum 
x-rays would be expected. This latter possibility has been investigated by Feather 
and Dainty (1944) using critical absorption and coincidence absorption methods. 
They concluded that K-capture did not take place in more than 15% of the dis- 
integrations. However, in view of the fact that they reported a nuclear y-ray of 
energy 65 kev., present to the extent of about one per disintegration, it was thought 


desirable to re-examine the problem. 
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§8. CRITICAL ABSORPTION MEASUREMENTS ON **Au X-RAYS 
IN THE SPECTROMETER 

It was thought that critical absorption measurements could be made much more 
accurately in conjunction with the spectrometer, because not only is the Compton 
background on which the absorption is superimposed reduced very considerably, 
but also x-rays which differ in energy by more than a few kev. will be resolved, and 
their absorption followed separately. Moreover, the secondary x-radiation from the 
critical absorber may also be resolved and studied separately. Figure 5 shows the 
K-absorption edges for tantalum and tungsten, together with the Ka,, Ka, KB,, 
KB, x-ray lines for platinum, gold and mercury, as well as for tantalum and 
tungsten. Let us consider just the Ka, and Ka, lines. Those for platinum at 66-9 
and 65-2kev., respectively, both lie below the K-absorption edge of tantalum at 
67-5 key., and hence will not be critically absorbed in either tantalum or tungsten. 

The Ka, of gold at 67-1 kev. will similarly not be critically absorbed in either, 
but the Ka, of gold at 68-8 kev. will be critically absorbed in tantalum but not in 
tungsten. Both mercury x-rays at 68-8 and 70-6 kev. will be critically absorbed in 
tantalum, but only the Ka, in tungsten. This difference in behaviour of the 


lau Kar] Ka] | iKa, j 2 ra 
Pt Korg] Kay] | Ka Kg. 


w Ka, jt Ka, “elliceage | 


Ta Ka, Kay Ka; | K,|[K-edge | 
(tet nnel eer en revel e rene  eerh ees 


60 70 
Energy in keV. 


Figure 5. X-ray energies of Hg, Au, Pt, W and Ta. 


K x-rays of platinum, gold and mercury should enable us to distinguish between 
them, or even to separate out the effects of each, if they are present to a comparable 
extent. 

A set of piles of absorbers of the same thickness, 0-3 mm., and of diameter 5 mm. 
was made up with varying proportions of aluminium and tantalum, from all 
aluminium to all tantalum. A similar set of aluminium and tungsten absorbers 
was made up. 

The absorber piles were stuck on to the centres of one-inch mica discs of weight 
5:0mg/cm?. ‘These were placed in turn behind the 5 mm. diameter 1 mg/cm? lead 
radiator, which was also mounted on the front of a similar mica disc, so as to shield 
it from the source. ‘The sandwich was mounted in the radiator holder, and the 
aluminium tube containing the source was pushed forward until it came in contact 
with, and slightly bent forward, the mica discs (Figure 1(5)). This was to ensure 
that both absorbers and radiator always took up the same position in relation to 
source. ‘The 5 mg/cm? mica backing to the radiator is sufficient to stop any low 
energy photoelectrons from the absorber. Usually two absorbers with two 
radiators of nominally the same thickness were placed in two of the positions of the 
radiator holder and alternate readings taken. ‘The two radiators were afterwards 
calibrated against each other. An attempt was made to measure the Compton 
background by placing a plain mica disc in the third position of the radiator holder, 
but this was unsatisfactory because the curve so obtained, although similar in 
shape, did not fit the spectrum obtained with the radiator, away from the neighbour- 
hood of the conversion peaks. ‘This led to an uncertainty in the actual heights of 
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/ radiator, the backing of which was sufficiently thick to stop the low energy photo- 
| electrons from the radiator, a background spectrum could be obtained which was 
identical with that for a radiator in its normal position, away from the conversion 
‘6¢ peaks. 
Cie The results of this set of experiments are shown in Figure 6, where spectra are 
““— given for the region of the first photoelectric peak, for the thickness of tantalum of 
0, 84, 168, and 500 mg/cm?. Besides reducing the height of the photoelectric 
peak, the introduction of the absorber also results in the appearance of a small peak 
at 42:5 +0-8kev., which may be ascribed to the L-conversion in lead of the K« 
radiation of tantalum. It is also to be expected that there should be peaks due to 
L-conversion in lead of tantalum Kf and the M-conversion in lead of the tantalum 
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Figure 6. Absorption in tantalum of low energy radiation from 1*8Au 
(1:1 mg/cm? lead radiator). 


Ka. Both of these are centred about 53 kev., i.e. almost exactly at the position of 
the photoelectric peak under investigation at 52 + 1 kev. Itis evident therefore that 
the true absorption of this peak is greater than would be deduced immediately from 
the heights of peaks at various absorber thicknesses. A correction for the effect 
of the secondary x-rays from tantalum may be made in the following way: if we 
are correct in assuming that the second of the two primary photoelectric peaks is 

made up in a similar way to the second of the <econdary photoelectric peaks due 

to tantalum, viz. as an L-conversion of KB plus the M-conversion of Ka, then the 

ratio of the two primary photoelectric peaks will be very nearly the same as that for 
the secondaries, since all the energies are very much higher than the L-absorption 

edge of lead. So the height of the first secondary peak may be used to deduce the 
height of the second. ‘They are of nearly equal heights. For large thicknesses of 
absorber, however, the height of the secondary peak is as much as two-thirds that of 
the primary, and consequently the difference becomes rather inaccurate. 

It was decided to repeat the experiments using a tin radiator, so that the K- 
conversion lines could be studied. This would have the advantage of moving the 
lines to lower energies, so improving the height-to-background ratio. An initial 
spectrum was measured under the same conditions as for the lead and gold 
radiators, but using a 2mg/cm? tin radiator. This is shown in Figure 4. Again 
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two conversion lines were observed in the region investigated, at 38 +0-°8 and 
48 + 1 key., but the second line was now only about 50% of the height of the first. 
They may be ascribed to the K-conversion in tin of the Kx and K8 lines from 
elements in the neighbourhood of gold. The use of the tin radiator therefore splits 
up the composite second photoelectric peak. Critical absorption experiments were 
then carried out in the same manner as before, using 0, 84, 168 and 335 mg/cm? 
tantalum absorbers and 196 and 390 mg/cm? tungsten absorbers. ‘The results are 
given in Figure7. Again the introduction of the absorber results in the appearance 
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Figure 7. Absorption in tantalum and tungste2 of low energy radiation from 1%Au 
(2 mg/cm? tin radiator). 

of lines at 28 + 0-6 kev. in one case and 29-7 + 0-6 kev. in the other, which may be 
ascribed to the K-conversion in tin of the Ka radiation of tantalum and tungsten 
respectively. Other smaller peaks at approximately 54 kev. could be ascribed to 
the L-conversion of the same K« radiation. Again a peak from the K-conversion 
of the Kf radiation is to be expected at 37 kev. for tantalum and 38 kev. for 
tungsten, i.e. just at the position of the primary peak, but this time it is only 50% 
of the first secondary peak, and the correction factor for the height of the primary 
peak is never large. It is evident from an inspection of the curves that the absorption 
in tantalum is very much greater than thatin tungsten. This at once suggests that 
gold or mercury x-radiation very much predominates over platinum x-radiation, 
or that of any other nucleus higher in atomic number than mercury. 

We must now consider the spectra in more detail. Since we are working not 
far from the cut-off of the counter window at 15 kev., considerable corrections must 
be applied for window absorption. However, this will only affect us as far as the 
estimation of the height of the secondary peak, due to the K-conversion in tin of the 
tantalum or tungsten Kf, is concerned. ‘The correction is greater with tin than 
with the lead radiator. It was applied on the basis of an exponential law of attenu- 
ation, with the exponent inversely proportional to the fourth power of the electron 
velocity.* Although the attenuation of the lowest energy peak, viz. the K-conver- 
sion of the Ka of the absorber, due to this cause is large (a factor of nearly three in 
the case of the tantalum Kz), the effect on the estimated height of the tantalum 

W. Bothe, Handbuch der Physik, 1933, Vol. 22/2, p. 41. 
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KB peak is such as to reduce the height of the primary photoelectric peak by only 
about 10° in the worst case. All the peaks contain more than one unresolved 
component, and to compare relative average intensities, the areas of the peaks 
should be used, divided by the mean Hp value, to take into account the resolution of 
the instrument. With these corrections the relative intensities of the two primary 
photoelectric peaks for the tin radiator become 3:1, which is just what would be 
expected on the basis of their being Ke and K8 x-radiations. Applying the same 
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Figure 8. Absorption in tantalum of low energy radiation from !*%Au. 


corrections to the lead radiator results, on the basis of the second peak being compo- 
site, the components being the L-conversion of a Kf and an M-conversion of a 
Ka x-ray, the ratio of conversions in the L and M shells is 2:1. 

Attenuation curves calculated on the above basis are shown in Figures 8 and 9, 
where the results for both the lead and the tin radiators are shown, together with 
attenuation curves expected on the basis of the radiation being the Ka, and Kx, of 
platinum, or gold, or mercury alone, in tantalum and tungsten, the latter being 
obtained from standard textbooks. Let us consider first the absorption in tungsten. 

The experimental points lie somewhat above the absorption curve for gold 
x-rays, or even for a frequency just at the foot of the K-absorption edge. However, 
it is possible that the experimental determinations on which the published 
absorption curves are based may be in error. The absorption in tungsten therefore 
gives a strong indication of gold or, perhaps, platinum K x-radiation. Mercury 
radiation can only be present to a small extent. ‘The experimental points for 
tantalum, however, lie just about on the absorption curve expected for gold 
K x-radiation and very considerably different from eit-er platinum or mercury 
K x-radiation. Again the absorption is rather smaller than would be expected on 
the basis of a 6:1 mixture of gold and mercury x-radiation, arising from external 
and internal conversion of the 0-411 Mev. y-ray, but nevertheless it is evident 
that platinum K x-radiation, if present, must be much less intense than the gold 
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K x-radiation. The upper limit which can be set on the K-capture process, 
from this evidence, is 6%. 
§9. CONCLUSIONS 

1. In addition to the 0-411 Mev. y-ray there are two higher energy y-rays at 
0-670 + 0-008 Mev. and 1-089 + 0-012 Mev. present to the extent of 1-5 and 0-4% 
respectively. The sum of the energies of the first two y-rays is approximately 
equal to that of the third. ‘This would suggest that there are energy levels of the 
198119 nucleus at 0-411 Mev. and 1:09 Mev. 

2. There is no evidence for y-rays in the region of 0-2 and 0-3 Mev. as reported by 
other workers. The upper limit which may be placed on their intensities, if 
present, is 2 and 4% respectively. 
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Figure 9. Absorption in tungsten of low energy radiation from }°SAu. 


3. There is no evidence for the reported y-ray at 0-158 Mev. 

4. There in no evidence for the nuclear y-ray at 65 kev. reported by Feather and 
Dainty. 

5. There is no evidence for K-capture, in agreement with Feather and Dainty. 
‘The maximum branching ratio for such a mode of disintegration compatible with 
our results is 6%, 
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ABSTRACT. A nuclear photographic emulsion which has been in contact with material 
containing traces of uranium and thorium shows on processing «a-ray stars arising from 
atoms which have undergone several successive disintegrations during the exposure. The 
relative frequency of the various multiple stars, which depends on the uranium/thorium 
ratio, is here evaluated on simple assumptions to show the possibility of determining these 
elements by star counting. 


HE use of the nuclear photographic plate in determining small quantities 
of uranium and thorium in rocks has been suggested (Joliot-Curie 1946), 
and subsequently results of the technique have been given (Poole and 
Bremner 1948, 1949, Coppens 1950). 
From the «-track density on the plate, we have an equation for cy and cy, 
the concentrations of uranium and thorium (Evans 1934): 


ud 
number of «-tracks per cm? per sec. = aud uR)cy + ES R) Crp 
4h ou 4h th 


where ny, mp, are the total numbers of «-particles emitted per second per gramme 
from uranium I and thorium, d is the density of the rock, 


y. =(range of «-particles in air)/(range in rock), 


and & R, Xu R are the sums of the ranges of the «-particles from the uranium and 
oe th 

thorium families. Radioactive equilibrium is assumed, and the actinium family 

isneglected. The necessary modification to the equation to allow for the existence 

of a minimum observable range p is obtained by replacing each R in the 

summations by the effective range R—p. 

Another equation is required. This can be obtained by finding the density 
of 2,3, 4,5 stars. These stars arise from atoms in or near the surface which have 
undergone several successive disintegrations in the time of exposure of the plate. 
The ratio of the density of, say, 3 stars to the track-density gives an equation 
for Cy, Cy, independent of (1). 

As a first step in obtaining such an equation we require the number per cm* 
of atoms which have emitted 2, 3, 4 particles in time 7. -particles may be 
regarded as a-particles of range zero. We consider the general case of a primary 
element A whose amount may be taken as constant during the exposure giving 
rise successively to elements B, C, ....N...., the respective decay constants 
being A, A,, Ax, ....A,.... and the elements being present in amounts v, 1, 
Vane yo.e eatoms per cm’. 

The number of A atoms breaking up in the interval di is vA dt, and if they have 
broken up at time ¢ (within the interval 0 to 7), the number of these which break 
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up further before time T is vA dt{1 — exp ({— A,(T—1t)}], therefore the total number 
of at least double emissions from A atoms is 


1 exp(—A,T) 
if va[1 —exp {-A,(T—t)}] a= sd(T — y a rl a F 
The number of double emissions from A formed in dt at time ¢ is 


d 1 exp( —A,#) ye ae 
aa(t-x + ny dt =v\{1 — exp (—Ayt)} dt. 


Therefore the total number of at least triple emissions from A in time T is 
T 
| v1 — exp (—A,f)][1 — exp {—A,(T—2)}] ae 
0 


atl v 1 
=vA< T—-— —- —-+ —~ ee =A) exp (— aa ee 
aq Lh SD ee ee ee 


2 


and by induction it can be shown that the number of at least r+ 1 emissions 


ne ie es RiNghson aan 
E,=v4T- 4+ i 5 OoOO m eXP (- nm. 
Me Sener 
The term (A;—A;) does not appear in the difference product. 
The number of 7 emissions per cm® from A is therefore 


vr r oe mee 
BaE vil Pau Go) a ee 


If the th member of the family is stable, so that there can be no more than 
(J— 1) «-particles emitted from an A atom, the number of at least (J— 1) emissions 
E)_» 1s the number of (/— 1) emissions per cm? from A. 

This approximates to vAT for large T’ as would be expected. ‘The expres- 
sions (2), on the other hand, rise exponentially to a constant value which will 
eventually be negligible in comparison with vAT. 

So far only multiple emissions provided by A have been considered. At the 
beginning of the exposure, however, there will be present the v,, v,.... atoms 
per cm? of B, C..... These will gradually disappear, giving rise to multiple 
stars ; new atoms of B, C.... will form to replace them, but the break-up of such 
an atom of B, for example, is a further disintegration of an A atom. Similarly, 
if the number of multiple emissions from A and the original v, atoms of B present 
at the beginning of the exposure have been found, we have only to take account of 
the multiple emissions from the original v, atoms of C, for new C atoms can only 
be formed from B or from A via B, and so on. In an exposure long compared 
with the half-lives of B, C...., such multiple emissions will be negligible in number 
compared with those from A. 

In the case of B, the number of at least double emissions is found by the same 
procedure as before: the number of atoms of the original v, which decay in dt 
at time t=v,A, exp (—Aj?) dt. 

The number of at least double emissions by the end of time T is 


a, | ; exp(—A,/)[1 —exp {—A,(T—2)}] dt. 


By a repetition of the process, the number of at least triple emissions by time 
T is found. 
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The number of double emissions is equal to the number of at least double 
minus the number of at least triple, i.e. 
= vhf exp(—AT) exp(--T) | _exp(-A7) \ 
(A2—Ay)(As— Ax) (Ar =Aa)(Ag3— Aa) (Ay = Ag)(Az—As)J 
In general, the number of (m—1) emissions corresponding to the number of 
N atoms formed in time T from "1 atoms of B, is 


exp (—A,T) 

Be nN UNA)” «ae iw 
the difference products as before not involving a (,—) term, as can be proved 
by an inductive argument. 

It should be said that many of the above results have been obtained (Bateman 
1910) by solving in the general case the simultaneous differential equations for the 
rate of growth of the members of a radioactive family by the method of the Laplace 
transform. It is necessary for the present theory, however, to terminate his 
equations by one expressing the formation of a stable member. 

We now turn to consider the problem: if there are N atoms per cm? which 
have emitted «-particles of range in air R,, R,....R,in exposure time T, how 
many of the various multiple stars per cm® are to be expected from these atoms? 
‘These ranges are supposed to be placed in ascending order of magnitude, so the 
‘subscripts 1,2....do not denote association with elements A, B.... . Letthelayers 
lying between depths 0 and R,, R, and R,, R, and R,,.... be labelled layers 
1, 2, 3,.... . It will also be convenient in the following to use R,, R,....R, for 
the ranges in the rock (=~! x range in air). 

Take layer 1: it is possible for all the 7 particles emitted from each atom to 
-emerge from the rock, that is, be registered as an 7 star. The probability of 
emergence of the R, particle at depth x intermediate between 0 and R, is 
4(1—x/R,) (Evans 1934). The number of 7 stars due to the atoms between 
x and x+ dx is, therefore, 


N x x x 
LT soa SE (Rela ae es cl 
(1 z)(! = (: 7 : 


-so that the number of 7 stars due to these atoms per cm? is 


N x x x 
= ——- ———' Ja, .("L.—-— | 
7, (i- Re) Rm) 


eee i Ree Re ad 


Boal jee eas eee ete 
Ube ie y ies 3 1 R;R; 4 1 R;R;R;, 
(i<j) (i<j <k) 
(-1) Rae 4 
Se oe (4) 


To find the density of (r — 1) stars from these atoms due to layer 1, we calculate 


the chance at depth x that the R,, Rs, Ry,....R,particles or the R,, Ry 1 Sate Aer 
OMtNCr eatin: coi, 50 2.0L the: R,, jm Re Ree oe ik 4, proceeding on the 
same atom shall pee " Integrating as above over x= 0 6 R,, we obtain 
N r sip RY Ss peer Ry" S ne 
Die | Ck Ce 2 ay x aR: a Cys 2 iF. RR; pa 
t (i<3) 
RR,’ 2 R: 
—1)y2 42S Ser ee os ) 
pies 7s ARR Re Gok: ey 
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This gives the density of (r—1) stars; it also includes, however, the + stars,. 
each of which is counted "C,,_, times in this process. So to obtain the density 
of (r—1) stars from layer 1, we must subtract from this "C,_, times the density 
of r stars. 

Using the same process for the (r—2) star density, we obtain 


N r r—-1 Ri? 1 r—2 Ri? nity, 
7] Ce oR, — Cy-3 Lom 452 =i Ca 3 z RR; Oe 
(i<3) 
z RR. 9 
SS \—2 Pi Caled] 
Beet ER Bees E BRE she (6) 
(i<j) 


which gives the density of (r —2) stars when we subtract "-1C,_, times the (7 — 1) 
star density, and "C,_. times the 7 star density. 

The same process can be followed down to the 2 stars. 

Layer 2 can now be considered: the R, particle cannot emerge from this layer, 
so that the maximum possible is an (r—1) star. The integration is now carried 
over the range of x from R, to R,; thus the density of (r—1) stars from layer 2 
is found by writing N/2"-1 outside the braces in (4) in place of N/2’, replacing 
R," outside the summations by R,” — R,” and allowing the summations to range 
from 2 to r. If we make the last two changes in (5) and (6) also, and replace 
N/2'-1 and N/2"-? outside by N/2’-? and N/27-% respectively, we obtain two 
expressions appropriate for calculating the (r — 2) star and the (7 — 3) star densities. 
From the first we must subtract ""!C,_, times the (r—1) star density, giving 
the (r—2) star density, from the second "-?C,,_, times the (r —2) star density and 
7-1C__, times the (r—1) star density, giving the (r—3) star density. We may 
proceed in this way down to the 2 stars. 

The (r—2), (r—3).... star densities due to lower strata can be found by the 
same method. 

If, say, two B-particles are among the 7 particles emitted, on writing R, = R, =0, 
the above expressions become just what would be obtained if the atoms had been 
treated from the start as emitting (r—2) particles, ranges R3, R,....R, in order 
of magnitude. The only part played by the B-particles is that they, to some 
extent, determine N, the number of (r—2) emitters per cm*, by reason of the time 
which elapses before they are emitted. In the following, in calculating these 
range functions, f-particles will therefore be ignored. 

All that is required has now been obtained: if we wish to calculate the ratio of, 
say, 3 star density to the total «-track density, we collect the 3, 4, 5.... emitters, 
multiply each by its appropriate range function and add, obtaining 


BND TAN Aa coe to Rapes Os 


where f; and /,’ are functions calculated as above. The total «-density is 4nATXR, 
where UR is the sum of all the ranges. The ratio is therefore independent of 
nand Aas is to be expected. Also the dimensions of all the f are those of length, 
so that if the ratio is multiplied above and below by ,, it is seen to be unaltered by 
writing everywhere air ranges for rock ranges. Correction for minimum observ- 
able range may be made by replacing these air ranges by the effective air ranges. 

This theory cau now be applied to uranium and thorium. In the case of 
uranium, no chain of quickly disintegrating «-particle emitters occurs except 
following radium. ‘This may be taken as ending as RaD for, though this is not 


A Method for Determining Uranium and Thorium in Rocks 29 


stable and will supply single «-particles, its half-life (22 years) is so long that the 
chance that any RaD atom formed during the exposure will disintegrate before 
the exposure is completed is negligible. Similar reasons justify our neglecting 
all members except those in the Ra-RaD chain. Further, the branching ratio 
at RaC is so much in favour of RaC’ that all the atoms may be regarded as dis- 
integrating by this route; also the half-life of RaC’ (1:5 x 10-4sec.) is so short 
that RaC may be regarded as turning to RaD by the simultaneous ejection of an 
a- and a f-particle. We therefore have the following scheme, treating radium 
as the primary element, which gives a maximum possible emission per atom of 
six (four «- and two f-particles) 


Re ee ny ep RT) 


The values for Aj, As, As, Ay calculated from the half-lives given in a recent review 
(Seaborg 1944) are as follows: A, =0-174 days, A, =328 days, A, =37-3 days-1, 
A,=50-7 days"'. The Table shows the density per cm® of the various kinds of 
emitter for different exposures, and the total number of «-particles emitted per 
cm?; it is calculated for nA=2-78 x 10° days“! corresponding to a decay constant 
for UI of 4-21 x 10-% day4, a concentration of 10-*gm. uranium per gramme 
of rock, and a rock-density of 2-6gm/cm*. In the sense used in this Table, the 
almost simultaneous ejection of an «- and a -particle from RaC is a single emission. 


Weeks 1 » 3 5 10 20 
Double emissions per cm® from 

Ra 5-98 7°75 8:26 8-46 8-49 8-49 

Rn 222 6:54 2:01 0-17 0 0 

RaA 0 0 0 0 0 0 

RaB 74:6 74-6 74:6 74-6 74-6 74:6 
Triple emissions per cm? from 

Ra 52°) 68-0 Wes 74-4 74-6 74:6 

Rn 16:3 4-83 1:49 0-12 0 0 

RaA 8-49 8-49 8-49 8-49 8-49 8-49 
Quadruple emissions per cm* from 

Ra 38-6 50-1 53-4 54-7 54-9 54-9 

Rn dfe3°< 102 ==1-46 x 10% 7 1°56 « 1041-60 10% +=1-60'x'10* 1-60: 10* 
Quintuple emissions per cm* from 

Ra S2060>5<10°1 62-42% 10% 4-269¢10? 8:12 «104 1°738><10° 3:73'x 10% 


‘Total number 


of particles 
emitted per cm? 1:56X10® 3:11x105 467x105 7-78x10® 1:56x10® 3:11 10° 


It is obvious that Rn is the only one of Rn, RaA, RaB, RaC with a small enough 
decay constant to play an appreciable part; the atoms of it present at the beginning 
of the exposure supply many quadruple emissions since its equilibrium propor- 
tion is large, and by its delay in disintegrating, it holds up the production 
of quintuple emissions from Ra. 

In the case of thorium, we have the following system, where half-lives are 
written below each element: 


Th—> MsTh, —> MsTh, —> RdTh —> ThX —> Tn —> ThA —> 


1:39 x 10!° yr. 6-7 yr. 6:13 hr. 1:9 yr. 3:64 days 5:45 soc. 0-168 sec, 
ala 
3x10 7sec. 
ibaa Ene Php: 


3:1 min, 


10-6 hr. 60-5 min. ne Tha” wa 
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We take RdTh as the primary element: it will be seen that this is justifiable 
even though its decay constant is quite large. Such a primary atom can emit 
seven particles, five « and two f, and, applying the principle seen to operate in 
the case of the uranium series, that for the exposures considered an element with 
half-life less than about half-an-hour may be taken as disintegrating as soon as it 
is formed, we have the following system: 


Rah > Th he ee 


The branching need only be taken into account now in the calculation of the 
range functions. 
We may reduce this further to 


RdTh—> ThX —> ThB > ThD. 

The error caused by this in the number of emissions in 1 week is 1:6% for 
quintuple emissions and less than 1°, for quadruple or less. It is also easily 
shown that less than 1-36 x 10-* of the RdTh atoms formed during a 20 weeks 
exposure have come from Th, and of these very few will have time to disintegrate 
further. That is, the decay of RdTh is made good almost entirely by B-emission 
from MsTh, even in long exposures, so that the objection of large decay-constant 
in the case of RdTh as primary element has been met. ‘This gives the following 
results for 4a and 5« emitters: 


Weeks 1 2 3 5 10 20 


(no. of transitions per cm? from ThX to ThD)/rA 3-68 = 4:85 5:16 5:26 5:27 -5:27 
(no. of transitions per cm® from RdTh to ThD)/nA 2:68 8:51 15:2 29:1 64:1 134-1 


For thorium, the branching ratio has to be taken into account in calculating the 
range function; 65°% of the atoms traverse one branch and 35% the other in any 
exposure; if f, is any range function for the set of ranges corresponding to the 
first branch, and f, the same function for the set corresponding to the second, the 
correct range function will be 0-65 f, + 0-35 fy. 

We adopt the following values for the ranges required (Joliot-Curie 1946) 
Ra 3-30 cm., Rn 405cm., RaA 4-66 cm., RaC’ 6-91 cm., RdTh 3-98cm-s 
ThX 4:28cm., Tn 5:00cm., ThA 5-64cm., ThC 4:66cm., ThC’ 8:57 cm. 


Thus the ratios to the total track density are, for the uranium family alone: 


Weeks 1 2 3 5 10 20 


4 stars per cm? (x 10 ~°) 2:78 4-16 4-89 Sy 6:13 6:42 
3 stars per cm? (x 10 ~*) 2°81 3:26 3-50 Sie 3-90 4-00 
2 stars per cm? ( x 10 —+) 1-02 1:05 S07 1-09 1-10 oil! 


and for the thorium family alone: 


Weeks 1 2 3 5 10 20 
5 stars per cm? ( x 10 ~$) 1-93 BE07) 3°67 4-20 4-62 4-84 
4 stars per cm? (x 10 ~*) 1:50 1-93 2:14 D385 2°50 DAS 
3 stars per cm? (x 10 -?) S78} 6:46 6:78 7:07 730) 7-41 
2 stars per cm? ( x 10 +) Ios 1-40 1-43 1-45 1:47 1:47 


No correction for minimum visible range is made. 
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The required equation additional to (1) is now easily obtained: if, for instance, 
a three weeks exposure is used, the total number of «-tracks per cm? in this time 
from the uranium is, from (1), 


2 
Ny x 60? x 24 x 21 xd. p eS 
4p U 
so that, referring tothe Table above, the number of 2 stars per cm? due to uranium is 


"4 2 
0-107 (Sees s R) a 
4p U 


Similarly, the number of 2 stars per cm? due to thorium is 


2 
0-143 Nyy, X 60 x24 x21 xd 5p ae 
4 Th F 


and we add these and equate them to the observed density of 2 stars. 

One can make use of the 3, 4 and 5 star densities to provide checks. 

The effect of allowing for the minimum observable range would be to decrease 
the ratios, but would affect the thorium ratios rather less than the uranium 
because on the whole the relevant thorium ranges are greater; thus there will be 
an increased accuracy in distinguishing uranium from thorium. 

The stars can only be recognized easily when the track-density is low, so that 
the calculation as it stands applies to rocks poor in uranium and thorium; even for 
such rocks the method would not be suitable if the radioactive elements were 
localized in very small grains, as often happens. For markedly inhomogeneous 
rocks, the method in combination with petrographic technique (Poole and 
Bremner 1949) would have its virtues: for example, one of the extra star-density 
equations could be used to eliminate p. 
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ABSTRACT. A careful analysis of the planar modes of vibration of ethylene and 
ethylene-d, is carried out in order (a) to determine finally the value of the C=C stretching 
force constant in these key molecules, (6) to find what interactions, particularly those 
between the two ends of the molecules, are real, and (c) to determine the vibrational forms 
and the potential energy distributions in the normal modes, and how they are affected 
by altering the potential function. Solutions are obtained both with and without correction 
for anharmonicity, and estimates are made of the effects of other corrections which may be 
necessary. The limits found for the C=C stretching constant are 9-336, 9-388 x 10° dyne/cm., 
with a most probable value of 9-343, the carbon atoms being uncorrected for anharmonicity. 
The zero-order constant is unlikely to be greater than 9:57. There is detectable real inter- 
action between the two ends of the molecule, for both stretching and bending coordinates. 
Interactions are discussed. 


§1. INTRODUCTION 


ECENTLY (Torkington 1949a, 1950 a) it was shown that the non-planar 
R vibrations of substituted ethylenes could only be explained satisfactorily 
if the C-H bending force constants varied with the substituents. Some 
evidence has been obtained for similar variations in the planar C—H bending 
modes (Torkington 1949b), and it is intended to apply a similar treatment to 
the skeletal vibrations of substituted ethylenes ; asa first step the planar vibrations 
of the parent compound tetra-chloro-ethylene were recently subjected to a 
complete analysis (Torkington 1950b). Now it becomes necessary to fix the 
value of the C=C stretching constant in ethylene itself, this involving a careful 
analysis of the potential function of ethylene and ethylene-d,, because the required 
constant 1s sensitive to some of the various small interactions which are theoretically 
possible. Ethylene has previously been treated quite fully, but none of the 
investigations is satisfactory, for reasons mentioned below. To avoid repeating 
the present calculations, the problem is approached from several different angles 
and necessary key results given. The procedures are outlined so that any 
required stage can be rapidly checked or adjusted. 


§2. PREVIOUS INVESTIGATIONS 


Brief summaries of previous investigations will be found in publications by 


Herzberg (1945) and Wu (1946); a review of work up to 1939 was given by 
de Hemptinne and Manneback (1939). 


The experimental measurements are summarized below: 


Levin and Meyer (1928) CoH,” Rea Gas 
Daure (1929) C,H, R:~ Liquid 
Dickinson, Dillon and Rasetti 

(1929) CA ahomeaGas 


Bhagavantam (1936) C,H, R. © Gas 
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Bonner (1936) C,H, R. Liquid 
Sutherland and Conn (1937) CD72 ic Re Gas 
Glockler and Renfrew (1938) CAE enh eelotcuid 
de Hemptinne, Jungers and cae $ Fa vs 
Delfosse (1938) oD, iquid (and partial spectra of 
all compounds C,H,,D,_,) 
Conn and Sutherland (1939) Ce Deel eGas 


Gallaway and Barker (1942) oa. ae ce 
2 4 . . 


Photographic infra-red bands have been reported by Badger and Binder (1931), by 
Ganswein and Mecke (1936), and by Thompson (1939); a weak infra-red band at 
approximately 800 cm™! was reported by Rasmussen (1942) and confirmed by 
Thompson and Harris (1944). The specific heat of ethylene was measured by 
Eucken and Parts in 1933, and by Burcik, Eyster and Yost in 1941. From measure- 
ments on a photographic infra-red band Badger (1934) deduced the dimensions of 


C,H, to be: rg q=1-04a., r¢2g=1-37a., HCH=126°. Thompson (1939) 
obtained 1-085, 1-331 a. and 118°. From an analysis of the rotational structure of 
infra-red absorption bands obtained with high resolving power in the fundamental 
region, Gallaway and Barker (1942) obtained the dimensions now accepted for 
Reet Pande, yea pep! -07 LAL, roog = 135345 HCH =119°55’ tom 

The vibrational assignments for ethylene and ethylene-d, now accepted as 
correct (Herzberg 1945) are based largely on the high-resolution work of Gallaway 
and Barker (1942), though the assignments of Conn and Sutherland (1939) were 
essentially the same. Previously, only the A, and B;, modes were certain. The 
B,, and B,, (planar) rocking modes were usually both taken to be in the region 
of 950 cm! (for C,H,), though it is the non-planar modes which are of this 
magnitude.* Earlier assignments for ethylene are discussed by Mecke (1932) 
and by Teller and Topley (1935); the former had assigned the non-planar modes 
correctly. Correlation with C,D, is discussed fully by Conn and Sutherland 
(1939). 

There have been numerous investigations of the potential function of 
ethylene. Mecke (1932), using a simplified valence force field, obtained a value 
of 8-0 x 10° dyne/cm. (quoted as 52 volts), for the C=C stretching force constant. 
Sutherland and Dennison (1935) treated the molecule as two coupled CH, groups; 
the potential function for each group was not of valence type. ‘The only cross-term 
corresponded to Arg_y.A@cy, interaction. Only the parallel, A, and By, 
vibrations were treated, and the force constants were assumed the same in both 
symmetry factors. The frequencies were adjusted to fit a relation deduced 
from this assumption, the adjustments (to the C-H stretching modes) being 
3019 3039 and 2988-2968. The value 9-79 x 10° dyne/cm. was deduced for 
the C=C stretching constant. Delfosse (1935) used a standard valence-type 
potential based on the symmetry coordinates of Brester as developed by Manneback 
(1935), which, together with the nomenclature, was carried over essentially 
unchanged by later European workers. He obtained a value 8-64 for kg. 
Bonner (1936) employed a potential due to Wilson. Here the bending coordinates 
were expressed in terms of the acute angle between a C-H bond and the C=C 


* Wu (1946) still attempts this assignment. 
PROC. PHYS. SOC. LXIV, I—A 3 
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axis. Interactions between the stretchings of adjacent C-H bonds and the 
deformations of adjacent angles were included, but none involving C=C 
stretching. Because of this omission, no real roots could be found for the cubic 
factor, but a best solution was reached by a variational method (which was not 
explained) giving kg_~=8-2. The treatment included all twelve fundamentals, 
but wrong assignments were employed for the B,,, B,, and B,, modes; the A, 
(twisting) mode was correctly assigned from the Raman line at 1656 cm”, 
reported here for the first time. Manneback and Verleysen (1936, 1937), using 
Delfosse’s potential function with eleven constants for the nine planar vibrations, 
obtained kg_¢ =8-674, this being chosen as the best value from the range 8-37 
to 8-97; negative constants were used for Arg_c.Arg_y and for Arcg_c.A%cx, 
interactions. Predictions of the frequencies of deuterated ethylenes were made, 
and the forms of the vibrations calculated. Considering the complexity of the 
potential function, the predicted spectra were not in good agreement with the 
values observed later by de Hemptinne, Jungers and Delfosse (1938). ‘Thompson 
and Linnett (1937), using a valence force field with one interaction term for 
Aro-c- Adcox, interaction, obtained a value kg_~. =9°8, associated with a negative 
interaction constant. Calculations of the frequencies of C,D, by Conn and 
Sutherland (1939) showed that this potential function was not as satisfactory as 
previous ones with a smaller C=C stretching constant. Fox and Martin (1938) 
used in addition a constant for interaction between the stretchings of adjacent 
C-H bonds, but gave no value for kg_g, which was eliminated from their 
equations in order to study the CH, groups. Lemaitre, [chang and Manneback 
(1937) and Tchang (1938), using the maximum number of interaction terms and 
one constant to correct for anharmonicity, were able to predict the frequencies 
of C,H, and of the six deuterated ethylenes, in so far as their spectra were complete, 
with an error less than 0:5°% in general and always less than 1:5%. The 
anharmonicity correction was introduced because it was obvious that the different 
degrees of anharmonicity in motions involving hydrogen and deuterium atoms 
were the chief cause of the discrepancies in the predictions of Manneback and 
Verleysen. In the method of Tchang, the masses my and my are replaced by 
effective masses my’ =My, p+ Cmy, »'”, where C is an anharmonicity constant 
for the isotopic pair; it is given a value such that the isotope product rule for the 
frequencies of each symmetry factor is obeyed. ‘The formula can be deduced from 
the expression for the energy levels of an anharmonic oscillator (see Kilpatrick 
and Pitzer 1947). A review of this work on ethylene and deuterated ethylenes 
was published by de Hemptinne and Manneback in 1939; the value quoted there 
for the C=C stretching constant was 9-356. In view of the fact that the B,, 
and B,, bending modes were assigned incorrectly, and that the dimensions of 
Badger were used in the calculations, it must be supposed that the relatively 
high accuracy of frequency prediction attained by Tchang is due at least in part 
to accidental coincidence. Only one analysis has been carried out with the 
correct assignments and dimensions (Kilpatrick and Pitzer 1947)*. Results for 
C,H, and C,D, were obtained both with and without Tchang’s correction for 
anharmonicity, using only the one cross-term employed by Thompson and 
Linnett. ‘The anharmonicity constant is calculated by fitting the A, factors 


* Apart from preliminary calculations by the author (Torkington 1949 b); see also Linnett, 
Heath and Wheatley (1949) and § 7. 
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perfectly, but its use does not greatly improve the fits for the other symmetry 
factors. ‘The values obtained for kg_ were 9-66 and 9-23 respectively, with 
and without the correction for anharmonicity. Two obvious criticisms of the 
treatment are as follows: (i) If it is granted that small differences between 
corresponding constants in different symmetry factors are to be allowed, then there 
is no justification for assuming zero for all except one of the off-diagonal elements. 
of the symmetry force constant matrix, for a difference between two corresponding 
~ constants is equivalent to an interaction constant in the expanded valence-force 
potential function. (ii) The Raman data for C,H, refers to the gas, that of C,D, 
to the liquid; some correction is surely necessary before evaluating the 
anharmonicity constant. 


$3: THE VIBRATIONAL ASSIGNMENTS AND POTENTIAL FUNCTION 


The assignments given by Herzberg were accepted as a working basis; they 
are summarized in Table 1; all values are in cm™!. 


Table 1. Vibrational Assignments for Ethylene and Ethylene-d, 


(Planar Modes) 
Ag Big Bou Buu 
CeH, CD, CoHy C2D, CoHy C2Dy C2H, CD, 
Vo=C 16232 Sree oeS 
pete B19: sme 251 s0 75-2304. = 3405-5 23454, 2999-5. 220082 
ioe Den 1342-4 981 1443-5 1077-9 
Deca oocrn 1050 883 995 7A 


All the frequencies of C,H, are taken from measurements of the gas, with the 
exception of the B,, stretching frequency; the value for this mode adopted here 
is taken from a weak Raman line of the liquid. This latter value is preferred 
to 3272-3 cm“, which occurs as a very weak Raman line of the gas, among other 
reasons because it allows a reasonable fit with the isotope product rule. The 
B,,, bending modes of both compounds are deduced from an infra-red combination 
band, measured for the gases. ‘The B,, and B,,, frequencies of C,D, are taken 
from measurements of the gas; the B,,, bending mode is in some doubt. ‘The 
value given here is that deduced by Herzberg by application of the isotope relations 
to the corresponding frequency for C,H,. Various maxima have been reported 
in this region, at 723-4, 727 and 740 cm“ (Herzberg 1945); the very strong B,y 
maximum at 720 cm™ is too close for a certain value to be obtained. The B,, 
and A, frequencies of C,D, (with the exception of the B,, bending mode, see 
above) are taken from the Raman spectrum of the liquid. There are no data 
available on liquid-vapour shifts for C,D,, but a certain amount for C,H, has 
been summarized by Wu (1946). The probable Raman frequencies of liquid 
5H, are 1620, 3008) 1341 (4A,) cm. 

The coordinates used are identical with those employed for tetra-chloro- 
ethylene, (Torkington 1950b), and the same order of frequencies is adopted. 
In general, formulae will not be repeated and the paper should be referred to. 

The atomic masses used were: my =1-00813, mp =2-01473, mg =12-00398. 
Conversion factor: 472c?m, =0-70691, gm.cm? sec. 


~ 
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§4. PRELIMINARY ANALYSIS 
Results for the three second-degree factors, assuming the off-diagonal elements 
of the force constant matrix to be zero, have been given previously (‘Torkington 
1949 b) and will not be repeated here. Comparison of the two sets of constants 
obtained (a) by solving explicitly and (6) by factoring off the C-H and C-D 
stretching modes by the method of Wilson (1941) suggests that the error intro- 


12 
(@)C2Hg I 


< 
+ 
Ss 


d,, 
aly 


oo 
nn “SN wo wo 


p- 


(dyne/cm x10") 
So 
dyne/cm. x 10°) 


22 


g 
012 


d 
oS 
oS 
oe 

i= 
Ss 


0:07 0-08 


0:06 
-05 =0:4 0:30 0:2 0d =0;3) 50:61 0:4 e072 0 
dy, (dyne/cm.x 10°) 


Figure 1. Allowed solutions for the force constants of (a) C3Hy, (6) C2D4 
from the second-degree approximation to the Ag factors. 
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Figure 2, Distributions of the potential energy in the second-degree approximation to the 
Ag factors for (@) C.Hy, (6) C.D. 
duced by factoring off these modes in the A, factor will be quite small. To 
obtain insight into the various allowed potential functions for each molecule, 
therefore, the A, factors will be treated as second-degree equations. The pro- 
cedure has been standardized (Torkington 1949c); the results of the analysis 
are conveniently summarized in Figures 1, 2 and 3. By correlating the two sets 
of constants it is seen that d;, must be Peereed to its range in ethylene, d,, to 
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the range 0-098 to 0-071, and hence d,, to 8-5 to 10:5. The potential energy 
distributions and displacements should be noticed. No common solution of the 
two second-degree equations is to be expected, since there will be some error 
introduced in factoring off vg_p, but we may obtain a preliminary fit for ethylene 
by correlating the A, and B,, factors; assuming no interaction between the two 


(b) CD, 
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Figure 3. The directions of the displacement vectors in the second-degree approximation to tne 


Ag factors for (a) CyH,, (6) C.D4. 


(B;= =arc tan (y,/21)(%i=), with the numbering and coordinates as for C,Cl, (Torkington 
1950b); || and 1 signify along and perpendicular to rc—u, p.) 


Coordinates: A,, C=C stretching; A», CH, (CD,) deformation. vi, and vz, are 
the modes of vibration in which d,, has respectively the greater and the smaller value for 
a given value of d,.. (See Torkington 1949 c.) 


ends of the molecule, the CH, deformation constant dj. will be the same in 
both factors. Carrying over the factored constant for dgq,, we obtain solutions at 


(a) djg=—0161, d,,=9-784 
ed g—=— 317, d= 10-86 


but for C,D, the B;, CD, deformation constant is less than the minimum A, 
constant, and there is no real solution. If we attempt to fit the two second-degree 
equations with common constants d,,, d,,, we obtain the two solutions 
(@)d,.= —0-1223, Hal 3 oD) nos = 0-0 7778 (CoE) 8 007527, (ESD, 
(b) d,.= —0-4206, Gi 10:53,.) tag = 0-07842(CH1,)3) ~ 0:07639(C.D,): 
It would be expected that the greater anharmonicity in a mode involving displace- 
ment of hydrogen relative to the same mode in the deuterium compound would 
lead to the uncorrected constant having the higher value in the latter; data on 
the increments have been given previously (‘Torkington 1949b). Inspection of 
solutions (a) and (b) above shows that the factoring of vg_p» must introduce an 
error which more than compensates for this effect. That this is so is obvious 
from the following: 


C,H, C,D, 
X,AcoA;/| A| 0-84.26 0-8701 (from cubic equation) 
AjA./| A’ 0-6489 0:-6275 (from second-degree equation) 


yo = (A,Aq/ | A’ | a d49”)/dy1, 
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where | A| is the determinant of the matrix A, = 1286, and |A’| is the determinant 
of the reduced matrix A’, = 64e2/(1+ 2c), (e=mc/my,p). Referring to the energy 
curves for ethylene, we see that the component in the potential energy from the 
cross-term is negative in vg_g and positive in dgy, from d,,= —0-12 to — 0:23, 
positive in vg_g and negative in doy, from d= —0-34 to —0-45, and negative 
in both modes from d,,=—0-23 to —0-:34. Over the same range in C,D, the 
cross-term is always negative in v¢_g and only positive in dep, for |dy2.|<0-24. 
But it should be remembered that these results will be more sensitive to the 
factoring approximation than those for ethylene. Ifit is argued that the cross-term 
must always tend to reduce the potential energy, then dj, will be restricted to the 
range —0-24 to —0-34. 

It can be shown that the approximate corrections to be applied to constants 
obtained from a reduced matrix A’ where this is of second degree are as follows: 


Aj," rs ! , Ay; ( A; ) } 
Nip 2 Mga da, 
wast Ay Age & = 7 au - Ay,A33 ry 4 r3 2 


A,,.A r r 

Adis = UT aes |] 2 a 3 |e / 
1 2A Aas LQ Az) Ar—As)_} 

where the d,,’ refer to the reduced factor and other quantities to the initial 


cubic factor. For C,H, and C,D, these equations yield the following results 
expressed as percentages: 


Cee Cen 
Ad.(Ro_c) 1:64 6-41 
Adys( fon.) 0-11 0-365 
Adj 1-32 4-11 


Note that the actual sign of the increment to the interaction constant is negative. 
Clearly, factorization of vg_y gives a good approximation for ethylene, but, as 
expected, the approximation for C,D, is not nearly as close. A factored high 
frequency is in general associated with interaction constants which retain the 
factored coordinate rigid in the residual modes (Torkington 1948, 1949 d); 
in C,H, and C,D, there are such arbitrary constants for Arg_¢. Arg_y, p interaction 
of magnitudes 0-39 and 0-68 respectively. The corrections to be applied to 
Re y,p When these are reduced to zero are proportional to these values, and 
inversely proportional to the element in the kinetic energy matrix for C-H, D 
stretching, so that the correction for C,D, will be about four times that for C,H,. 


§5. THE CUBIC Ag FACTORS 


Proceeding to the unreduced cubics, we first attempt solutions in which 
constants are carried over from the B;, factors. "The equations have been given 
for C,Cl, (Torkington 1950 b, equation (20)); here the constants k, f and d are for 
C-H stretching (C—D stretching), CH, (CD,) deformation, and for interaction 
between these two coordinates, and dj), d,, and d,; are for C=C stretching, 
Argo. Arc_y interaction, and for Arg_¢.A@oy, interaction (and for the corres- 
ponding interactions in C,D,). Using observed frequencies from Table 1, the 
constants on the right-hand side of the equations take the following values: 


C,H, (Hb). 
Cy = 31 (AL + Ag + Ag) — Ay’ +A’) 8-953 8-207 
Cy! =(1/8e)[(AyAg + AgAg + AgAg) — Ay/Ag’] | 13-943 13-724 
Cy =(1/128e?)A, AoA; 0-84.26 0-8701 
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with the force constants in units of 10° dyne/em. With d,.=d,;=0 the 
corresponding values obtained for kg_¢( = dj) are: 


C,H, C,D, 
From C,' 8-95 8-21 
From C,’ 9-01 8:58 + x 10° dyne/cm. 
From C,’ 9-60 9-39 


_ The values from C;,' will be quite good approximations, since the left-hand side 
of the corresponding equation contains only squares of dj, and dj3. If it is 
permissible to carry over the B;, constants, then the differences (C,4—C,?) 
will arise from anharmonicity; for zero-order frequencies there would be only 
one value for each constant C,’. The term A,A,/|A| in the B,, factors (corres- 
ponding to C,’) takes the values 0:08781 (C,H,), 0-09262 (C,D,), from which it 
appears that there is considerably more anharmonicity in the B,, than in the 
A, modes. If the anharmonicities were the same in both C-H (C-D) stretching 
modes, which will have the greatest effect, they would be eliminated in the 
constant C,’. Since (C,’"—C,’”) is so appreciable, and the carrying over of the 
B;,, constants will be at least a good approximation, apart from considerations of 
anharmonicity, the only possible conclusion is that the degree of anharmonicity 
is much greater in the B;,, than in the A, modes. ‘The larger negative correction 
to C,'# to reduce to zero-order frequencies would then bring C,’# and C,’? toa 
common value. Confirmation, if considered necessary, is obtained by adjusting 
the B,, frequencies of C,D, only, so that the constants for both molecules are 
identical (and the isotope product rule satisfied). ‘The frequencies are then 
corrected to ‘C-H anharmonicity’. The increments are: Avg_p=—39-3, 
Adcp,= —9:3cem. Using the corrected frequencies, the corresponding incre- 
ments to the A, frequencies can be found, using the first two equations. ‘They are: 
Avg_p = + 12-4, Adcp,=—18-6cm-. These values give AC,'? = —0-0231, which 
is satisfactory. ‘he degree of anharmonicity will be proportional to the difference 
(d,;> —d,;") (referred to subsequently as the anharmonicity increment) and hence 
will be greater the greater the negative increment Av,_y required to bring 
(d,;?—d;;") to zero. The absolute frequency increments given above are probably 
of no great significance, but the results can be interpreted as implying that there 
is much more anharmonicity in the B,, stretching mode than in the corresponding 
A, mode, though probably not in the bending mode. 

With the observed vibration frequencies for ethylene, and substituting 
k=1-2712, f=0-06907, d=0 from the B,, factor, we obtain the solutions 
d,,=9-962, d,z=—0-158,; and d,,=11:081, d,,=—0-320), both going with 
di, = —0:044,. The solution with the smaller values of d,, and d,; may be 
compared with the approximate solution obtained from the reduced equation and 
the factored B,, constant. Because of the large component from vg_y which occurs 
in the cubic, and the unknown anharmonicity, it seems better to take the solution 
from the reduced equation as a standard for A,/B,,, fitting. If we attempt to 
solve the cubic for C.D, similarly, no real roots can be found. This does not 
necessarily mean that carrying over the B,,, constants is not a reasonable approxima- 
tion, because of the following considerations. By eliminating d,,, the following 
two second-degree equations in dj, and d,; are obtained: 


(dy2—k)* fh (dis + 29/3f)? ee: (diy — hk)? 4 (d,3 + 2/3f) at 
Ky K,/4 ‘ K, (f/k)K 
where K,=h?+48f2+(k+4f)Cy'—Cy'; K,=R? + 12kf+kCy' —(C3’/f). 
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These are equations of ellipses having a common centre; the condition that 
there should be real solutions of dj», d,3 is that they should intersect, viz. K,< Ky 
with kK,>4fK,. With C,H, this condition is just satisfied, and there are four 
real points of intersection in close pairs; with C,D, the condition is just not 
satisfied, and one ellipse lies wholly within the other. It is possible to separate 
the variables d,., d,3, obtaining the following simple equations: 


dyo* — 2Rdyy — (RC — RC,’ +4C;')/(k-4f) =, | 


dys? + 4/3fdy3 + (Cy! —fCy' + 4f?Cr')/(k —4f) = 0. 

A shift of —6cm~ in the A, C-H stretching mode reduces d,, to zero, while 
having only a negligible effect on the other constants. For the moment, then, 
we will attach no significance to d,,. It should perhaps be mentioned here that 
part at least of the large increment (C,’#—(C,'?) may be attributable to Fermi 
resonance; there are a number of other Raman lines in the region of vgq whose 
intensity Herzberg ascribes to resonance with the A, C-H stretching mode. 
A shift of —28 cm7! would reduce C,'" to the order of C,'?. Such a large shift is 
extremely unlikely, but certainly one of 6cm™! seems possible. This uncertainty 
as to the reality of small interaction constants is of course quite general; many 
small values quoted in the literature are unaccompanied by any data on the 
frequency-shifts associated with their elimination. With reasonable corrections 
for liquid—vapour shifts to C,D, there are still no real roots when the B;,, constants 
with d=0 are carried over. The procedure adopted at this stage is to solve for 
increments Ak, Af and Ad,, with respect to solutions satisfying ethylene, so that 
the corrected constants satisfy the equations for C,D,, the ethylene solutions 
being taken from Figure 1 and corrected as explained in the last section. The 
following set of equations is obtained for the increments, from equations (19) of 
the paper on tetra-chloro-ethylene (‘Torkington 1950 b): 


(ental) 4(2e, +3) 4/3 Ak 
fd, + 8(ep + 2)f+ 44/3d,5} 4{dy,+2(ept2)k}  4(4/3Rk—2d)3) Af 
fay kd, = 2kd,s Ad, 


(2g —ep)(R + 4f) —(Cy7 —C,?)} | 
(8(ex ot ep )Rf— (Ce = Cx) 
(ee ee) | 


wher? €y p=Mc/My,p and Cy =}(Ay+AQ+A3), Cy =(AyAg + AjAg + AgAg)/Be, 
Cy =AjAqA3/128e", the superscripts H, D referring to C,H,,C,D,. 

Since it has been found a reasonable approximation to the factor for ethylene, 
the same constant k is used throughout; the value taken was 1-277, this being 
obtained from the A,/B;, fit, using a correction for kinetic coupling (Torkington 
1949d,e). The results, in units of 10° dyne/cm., are shown in the following 
Table; since the treatment is only approximate, though the approximation will 


du 8-67 9-48 9-60 9-95 10-45 

k 1:277 — = = te 

y 0-07787 0-07167 0-07090 0-06906 0-06722 

ee ee —0-139 —0-144 —0-163 —0-206 
Ak —0-0028  —0-0018 —0-0016 —0-0011  —0-0004 
Af 0-00103 0-00235 0-00252 0-00301 0-00387 


Ady. 0-057 000095 —0-0061  —0-0239  —0-0453 
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be quite good, the observed frequencies of C,D, were employed in the fitting. 
The increments obtained will of course be the anharmonicity increments intro- 
duced previously ; since no such increment is to be expected in d,,, this does not 
appear as a variable; d,, and d are neglected. 

An obvious first point in the table is the consistently negative increment Ak; a 
small correction for liquid-> vapour to C,D, would change the sign of this increment, 
but unless there were an additional positive correction to vgp, or a negative 
correction to vey, such as the Fermi resonance correction discussed previously, 
it would still not be sufficiently positive (compare the data given in a preliminary 
analysis of these molecules by orkington (1949b)). The trends of the 
anharmonicity increments for f and d,, suggest that the correct solution for the 
pair of molecules will be in the region of d,,=9-5, where d,, changes sign and 
becomes small and negative. d,, will not have a smaller value, for then |d,;?| 
would be less than | d,,"|, and it will not be as great as 9-6, for then the increment 
for d,, is as large proportionately as that for f; it is unlikely that this would be the 
case, since d,, 1s for an interaction involving a coordinate (C=C stretching) 
whose anharmonicity is invariable. The anharmonicity increment for f amounts 
to 3 to 4% (cf. 2% for the B;, factor), and the absolute values for f exceed the 
B,, constants in both molecules by about 0-003. The results are consistent 
within the degree of approximation to be expected for this variational treatment, 
and suggest that there is a real positive constant for interaction between 
deformations of CH, or CD, groups separated by a double bond, and that the 
correct values for d,, and d,; are near 9-5 and —0-14 respectively. Note also that 
there is probably less anharmonicity in the B;, deformation than in the corres- 
ponding A, mode—certainly not more—and so there must be much more 
anharmonicity in the B;,, stretching mode than in the corresponding A, mode, 
or a negative Fermi resonance correction to the latter, confirming previous 
conclusions. Correction to zero-order frequencies would probably not eliminate 
the constant f(A,)—/(B,,) for interaction between the deformations of the two 
CH, (CD), groups. 

So far, we have solved for ethylene using the constants from the B,,, factor, 
found d,, negligible in the first instance, obtained a series of approximate solutions 
for C,H, and, from the anharmonicity increments required in order that the 
constants should satisfy the equations for C,D,, deduced that f(A,)>/(B3,) and 
that the correct potential function representing the observed frequencies has 
di,29°5, digo — 0-14. Since any correction for anharmonicity is likely to be 
found arbitrary, it is convenient to obtain ranges of explicit solutions at this stage, 
in case it is subsequently found preferable to correct the constants directly, 
rather than the frequencies or atomic masses. The procedure adopted was as 
follows. The constant f can be expressed as a function of k only, d and d,»y being 
neglected, and d,, as a function of k and d,,: 


f=(1/A)[(Ze + I)k—C, + C,/k] — C3/R?, 
dys = (1/44/3)[2(€ + 2){Cy — (2e + 1)R} + {(2e + 3)/Ri(4C3/k — Co) — dul. 
With a given value of k, f may be obtained immediately and d,,; expressed as a 
function of d,,, of such form that the identity involving the constants and C; is 
satisfied for all values of d,,, but the identities for C, and Cy are only satisfied at 


the correct value of d,,. This latter is soon found. The calculations may be 
repeated for as many values of k as required. (This is a rapid and exact method 
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for solving equations from cubic factors.) The results are conveniently plotted 
against d,,; they are shown in Figures 4, 5 and 6, (a) for the observed frequencies 
from Table 1, curves (i) and (iii); (6) for C.D, with empirical corrections for 
liquid-vapour frequency shifts, curve (iv); (c) for C,H, with an empirical 
correction of —10cm™! to vey, for Fermi resonance, curve (ii), Figure 4; (d) for 
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Figure 4. Solutions for k in the Ag factors ; 


(i) C,H4, observed frequencies; 
(11) C.Hy, with Avca=—10 cm; 
(iii) C,D4, observed frequencies ; 
(iv) C,D,, corrected to gas. 
Coordinates : A,;, C=C stretching ; 
C-H (C-—D) stretching; 
deformation. 


Aa 


As; CH, (CD,) 


d, (dyne/cm x 10°) 


0-078 


0-076 P 
0:074 
0-072 
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F (dyne/cm. x 10°) 
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d.,, (dyne/cmx 10°) 


Solutions for f in the Ag factors: 
(i) C.Hy, observed frequencies; 

(ii) C.Hy, with Adcu, = —Avc=o= +10 cm}; 

(iii) C,ID4, observed frequencies; 

(iv) C,Dy,, corrected to gas. 


Figure 5. 


Coordinates: Aj, C=C stretching; Ag, 
C-H (CD) Nstretchina sy AGse Choe Ds) 
deformation. 


9-8 10:0 


d,, (dyne/cm. x 105) 


Figure 6. 


Solutions for dj, in the Ag factors; 


(i) C,Hy, observed frequencies; 

(ii) C.Hy, with Adcu,= —A vo=c= +10 cm-?; 
(iii) C,D4, observed frequencies; 
(iv) C,D,, corrected to gas. 


Coordinates : 


Ai, C=C stretching; A2, C-H (C-D) stretching; A3, CH, (CD.) deformation. 
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‘C,H, with empirical corrections of +10, — 10cm! to dc, and vg_g respectively, 
to show the effect of Fermi resonance between these modes, curve (ii), Figures 5 
and 6. Note should be taken of the relative magnitudes of k and the effect of 
the corrections (b) and (c)—see the previous discussion—and of the trends of the 
anharmonicity increments for f and d,3. For f the trend is always positive; hence 
there is no solution for the pair of molecules with common constants without 
correction for anharmonicity. For d,, there is a change of sign which can be 
used to select solutions. Further, the correction for Fermi resonance between 
voc and doy, in C,H, increases d,, and decreases f, although the frequency 
increments are in the opposite directions. This is due to the importance of the 
coupling term; |d,,| becomes larger. The results are summarized below, all 
values being in units of 10° dyne/cm. 


Correlation : 
C.H, C.D, C,H, C.D, C.H, C.D, C.H, C.D, 
(obs.) — (obs.) (obs.) — (corr.) (corr.) — (obs.) (Corr) mee (Comme) 
if 0:0717; 0-0740 0-:07130 0-0740 0:0708 0:0729, 0:07050 0-0729, 
dys 9-455 Ihe 9:69 9-75 
dis —0-1340 —0-136; —0-157 —(0-160 
“obs.” refers to the observed frequencies in Table 1, (C,H, gas, C.D, liquid); 


‘corr.’ refers to corrected frequencies, for ethylene with Avg_g=—10cm"1, 
Adoy,= +10cm4, and for C,D, with the following empirical corrections for 
liquid—vapour shifts: Avg p=+6°3, Aveo_¢=+2:9, Adgp,= +0-75cm!, these 
being deduced from those for C,H, on the assumption that for isotopic molecules 
such shifts will be proportional to the square of the frequency. k is omitted from 
the table. Correction of v¢_y by — 10cm has only a negligible effect on the fits 
for the remaining constants. It may be noted that the corrections to C,D, give 
shifts of the solutions as follows: Af= —0-00030 to —0-00045 for C,H,, zero for 
C,D,; Ad,,=+90-06, Ad,, = —0-0025to—0-0030; and that the corrections to 
C,H, give Af=—0-00080 to —0-00095 for C,H,, —90-00105 for C,D,; 
Ad,, = + 0:24, Ad,;= —0:023. The solutions with common values of d,3 have 
been given; in fact, this constant would be expected to have a small anharmonicity 
increment. It would probably not exceed 2%; the increment for Bs, deformation 
amounts to 1:7%. The shift of the solution for Ad,,= — 0-003 is Ad,, = + 0:06, 
Af = —90-00040 to —0-00030 for C,H,, —0-00030 to —0-00020 for C,D,. The 
anharmonicity increments for f are between 3 and 4%; all the absolute values 
exceed the B,,, constants, by 0-0027 to 0-0014 for C,H,, 0-0037 to 0-0027 for C,Dy. 
(These differences apply to the solutions in the table for common d,3; additional 
correction to d,, leaves the result essentially the same.) The values of k at 
d,,=9-45, with the observed frequencies, are 1-2772 (C,H,); 1:2759 (C,D,). 
The slopes (0k/dd,,) are of the order —0-002 to —0-0015 for C,H, —0-005 to 
—0-003 for C,D,. The increments for the corrections (a) of —10cm™ to vg_y; 
(b) of the liquid—vapour shifts to C,D,, are respectively —0-0086 and +0-0078. 
‘The anharmonicity increments at d,, =9-5 referred to the corrected C,Dy constant 
are 0-0063 and 0-0149 for the observed and corrected C,H, frequencies. For the 
B,,, stretching mode it is 0-0469. It would appear from these results that some 
correction for Fermi resonance to the A, stretching mode was necessary, but 
that the anharmonicity could not be as great as in the Bs, mode. Further, the 
increments to k required for coincidence with the B3,, stretching constants are 
respectively —0-0058, +0-0425 for C,H, C,D,, referred to the results for the 
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observed frequencies. The corrections available do not lead to consistent values, 
which suggests that small real interactions, either d or d,, or both, have been 
neglected. 

It is convenient to give here the derivatives determining the dependence of 
the frequencies on the force constants (Torkington 1949d, 1950b, Appendix IT). 
The solution used was the standard anharmonicity-corrected one with C =0-09307, 
obtained in § 6. 


(OA;/0d;;) 

k yy 1141 12 13 Dy 23 33 
1 48 C,H, 1:518 0-690 11:92 0:314 5-418 93°57 
ne C.D, 1-476 1-375 6:560 1-280 6-109 29°15. 
C.H, 0-140 —2-°517 0°852 45:19 —15-29 5-172 
2 C.D, 0:478 — 3-318 2:207 23-02 (5-311 10:17 
Cash 0-342 —0:173 — 5-840 0-087 2:944 99-61 
3 C.D; 0-046 —0-056 —1-839 0-070 2:273 74:13 


Coordinates: A,, C=C stretching; A,, C-H (C-D) stretching; A 3, CH2 (CD,) deformation. 
The derivatives are matrix elements; for off-diagonal force constants multiply by 2 for dj=dj_ 


§6. CORRECTION FOR ANHARMONICITY 

The only planar factor for which the vibration frequencies for both C,H, and 
C,D, in the gaseous state are certain is B,,; the anharmonicity constant in 
Tchang’s formula takes the value 0:09307. Using the data in Table 1, the B,, 
factor (for which only dcp, is uncertain) gives 0-089. These two values seem 
sufficiently close to justify accepting 0-09307 as the true planar anharmonicity 
constant. he procedure adopted was to fit the frequencies to this value, using 
the most probable increments. ‘The calculated values in cm~ are C,H,: B,, 3086, 
L036. Cee gps 2262-45 1518; 9.983.) Bye c2 3104 S745 ae ee eee 
value for 6(B,,) for C,D, may be compared with the reported maxima at 723-4 
and 727cm~'. ‘The liquid—vapour shifts for the A, C,D, frequencies from the 
above are Avg_p=11-4, Avg_g=3-0, Adep,=2-0cm™!. The last two must be: 
fair approximations, but Avg_p is probably too large; the data for ethylene give 
Ave_y=11cm™?. But the solution will not be very sensitive to v;_p except for k: 
and possibly d,,.; the point is considered later. The B,, factor is the most 
uncertain. The value obtained for C,H, by taking Ave_y=11 for liquid+vapour 
is 3086; a value of 1056 is obtained if the excess of (6(B,,)+5(Bg,)) over 
2047-0 (an infra-red combination of the gas) is assumed to be 4cm™!. The 
value 4cm™ is derived from that for 6(B,,,)+6(B,,), the infra-red combination 
band at 1889-6, assuming the liquid—vapour shift for 8(B,,) to be 1-4 cm. 
Similar arguments were used to fix 6(B,,) for C,D, at 874cm™!; the fitted value 
for v(B,,) then gives a liquid—vapour shift of 4-4 cm7!, which is probably about 
2cm! too small. 

With these assignments, and using C=0-09307 (equivalent to my’ =1-10158,. 
my =2:14683) the following solution is obtained : 


Ay Bre B, 

9-3429 0:2582 —0-1463 
1:3988 —0-0347 Here 0-100 ee 0-050 1:3647 —0-0347 ] 
0-08035 0-1103 ed | el 
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The constant d in the A, factor is carried over from the B3, factor; this seemed 


the best assumption to make. The dependence of the remaining constants on d 
is shown in Figure 7. The curves for k and f are symmetrical about d=0; fis 
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Figure 7. The dependence of the force constants of the Ag factor on d (=dy3), in the region of 
the exact solutions with corrections for anharmonicity in the motions of the hydrogen and 


deuterium atoms (C=0-09307). 
Coordinates; Aj, C=C stretching; A,, C—H (C—D) stretching; A3, CH2 (CD,) deformation. 


always greater than f(B;,), and k is always greater than k(B;,) for |d|<0-10. 
d,, and d,, have minimum values for small negative values of d, and are not very 
sensitive to any variation in d over this region; dj, varies appreciably, and changes 
sign at d= +0-030. 

It is of interest to compare the above results with those of chang (1938) 
who used effective masses my’ =1-088, mp’ =2-126, corresponding to C=0-08. 
In the A, factor he obtained a small positive constant in place of the negative d 
{d,,; of the above matrix) and a small negative constant for Arg_c.Aro_y 
interaction in place of the positive dj. In the B,, and B,,, factors his interaction 
constants are of the same signs as in the present solution; in the B,, factor the 
constant is of opposite sign. (Note that his angle-bending coordinates for the 
B,, and B,, factors have senses opposite to those used here.) Discrepancies are 
to be ascribed to wrong assignments for the B,, and By, factors, and to the 
molecular dimensions used. But they suggest that a careful investigation into 
possible variations of the small interaction constants should be made before the 
present results are accepted as final. This has been carried out, the steps being as 
follows. Using the exact solution for the B,,, factor of C,H, with zero interaction, 
the frequencies of C,D, were predicted to be 2190-1, 1082-9cm™. The error 


here is of the order of 0:5°9%. Next, the angle 9 =HCH required for an exact 
solution with zero interaction to be possible was calculated for the B;, and Bp, 
factors. The general condition for such a solution for a pair of isotopic second- 
‘degree factors is 
(Ajy’ gg — AyyAoe’)*ArA2/D) + {Aura + Az) — Ania’ +2')} 
x {Algg'(Ay + Az) ~ Asal’ + Az’) } =9, 
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where D=|A| and (A,A,/D) =(A;'Aq’/D’) is the isotope product constant; primes: 
refer to the isotopic species. If the angle @ is taken as a variable, the result for 
Coote CsD is 
cos? =[(hy + As)(e’ +1) — (hy +2s')(e + L)PIL{Qr + Aa) — Oa! +2’)? 

— Ble = ¢')*(AA2/D)I7, 
where B=256 for B,,, 64 for B,,. Substituting, we obtain 9@=112°33’ for 
B;,, 115°8’ for B,,. From the experimental evidence one would not on the 
whole be justified in taking the interaction constants in the second-degree 
equations as zero, and certainly not in changing their signs. 

In the above treatment of anharmonicity only the masses of hydrogen and 
deuterium are corrected; only one constant can be obtained. To investigate the 
effect of correcting also the carbon mass, the constant was assumed the same for 
all three atoms and the above calculations repeated for the A, and B3,, factors, 
giving 

AG B3y 
9-5691 0-2150 —0-1429 
1-4028 —0-0347 13737 wax 00347, 
0-08101 E se 


The value for C becomes 0:0988, somewhat higher than before. 

The additional correction has no effect on the B,, interaction constant which, 
as before, was.carried over to the A, factor. ‘The values of both pairs of constants 
k and f are raised somewhat, but the differences k(A,) — k(B3,) and f(A,) —f(B3y) 
remain essentially the same; d,, is decreased, but was certainly not due to neglect 
of carbon anharmonicity; d,3 is hardly affected, and d,, is increased by about 0-2. 
For comparison, the calculations for the A, factor were repeated with C =0-09307 
(for my and my only) and with an empirical correction Avg_g= + 25cm lin C,H, 
for carbon anharmonicity. ‘The increment for C,D, required for satisfaction of 
the isotope product rule is then + 26-4cm7!. he increase in kp, is about the 
same as that obtained with corrected carbon masses and observed frequencies. 


Lastly, the A, factor was solved (a) with Avg-g=—10cm™, (6) with 
Avgeg= —Adggq,= —10 cm in C,Hy, vg_p being adjusted to fit in both cases. 
The shifts represent possible Fermi resonance. 
C.H, 

1648-3. — -- — 3009-3 — 1613-3 — 1352-4 
C.D, 

1541-4 — — — 2258-0 — — 2265-3 — 


1:-4010 —0-0347 1:3876 —0-0347 1:3968 —0:0347 
0-08065 0-08003 0-07931 } 


Units are, as usual, in cm™ and 10° dyne/cm. 


9-5545 0-2870 —0-1384 le 0-2438 —0:1473 ee 0-2400 —0-1667 


Fundamental frequencies different from those used so far in this section are given 
in the appropriate columns of the force constant matrices. It will be seen that 
it appears to be impossible to eliminate the small positive constant d,. for 
Arcg_c-Arc_y interaction if the constant d is carried over from the B,,, factor; 
the dependence of the constants on d will remain essentially as in Figure 7. It 
remains to attempt to decide whether or not the carbon anharmonicity correction 
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is an improvement. ‘The only way available is to predict the zero-order 
frequencies from the uncorrected masses and the force constants obtained here. 
With C=0-09307, my and my corrected, and frequencies going with C,H, (obs.), 


ae B 1g B 2u B 3u 


C,H, 3147-0 1643-0 1390-4 3202-2 1117-2 3238-4 1035-7 3117-5 1502-6 
C.D, 2318-8 1530-1 1014-0 2361-5 915-9 2416-6 742-7 2263-5 1107-4 
With C=0-0988, my, mp and mg corrected, 

CHA, 3157-3 1660-9 1397-6 

a 2304s -1543°8>"1017-7 


The effect of the additional correction to the mass of the carbon atoms is to 
raise the frequencies somewhat, particularly vo_¢; the results with C=0-09307, 
applied to my and mp only, with the empirical correction for anharmonicity to 
Voc 1s to raise the zero-order frequency of the latter to about the same value as 
for C=0-0988 (the actual magnitudes are: 1667-4(C,H,), 1553-9(C,D,)), with 
the remaining frequencies practically as in the first solution. It might seem that 
the additional corrections for carbon anharmonicity were superfluous; an 
increment of 20cm? at 1,500 to 1,600 cm™!, obtained with the correction to my 
and mp, is about what would be expected (Herzberg 1945), while the 40 cm7! 
for C,H, with the extra Avp_j= +25 cm! seems definitely too large, though the 
possibility that some of this is derived from the 30°% CH, deformation in the 
‘ C=C stretching ’ mode must be considered (Figure 2). The increment for the 
(very nearly) pure CH, deformation in the B;, factor for C,H, is 59-1 cm, while 
in the A, factor it is 48-0 cm“; in the latter there is only about 75° CH, deforma- 
tion. The corresponding increments in C,D, are 31-0 (4,) and 29-5 (B;,); here 
there is no appreciable component from dgp, in the C=C stretching mode, but 
some contribution from vg_p, since the factoring of this coordinate produces 
considerable error. The increments to vgig in C,D, are 12-1 and 25-8cm 
for C=0-09307, 0-0988 respectively. ‘The carbon anharmonicity constant would 
not be the same as the constant for hydrogen and deuterium, though the results 
obtained suggest that it must be of the same order; it would certainly seem 
necessary to allow for one, but at present there is no alternative mode of approach. 

Finally, with a correction of —10cm™! to the A, stretching mode, the fitted 
C-D stretching frequency corresponds to a liquid—vapour shift of +7-0cm™, 
which must be very nearly correct. 


§7. INTERACTIONS IN THE ETHYLENE MOLECULE 

It is not intended to treat interactions in the ethylene molecule at any length 
here, but a few points may be summarized for comparison with C,Cl,. ‘The 
interactions between the stretchings of C-H (C—D) bonds* go with constants 
—(-0169 (adjacent), +0-0133 (cis-), +0-0549 (trans-), taking the values from 
the standard anharmonicity-corrected constants with C =0-09307 and using the 
previous table of relations (‘Torkington 1950b, Table 4). A correction of — 10 Cine 
to A, C-H stretching produces an increment — 0-0112 in k(A,), which increases 
the negative constant for interaction between adjacent bonds and decreases the 
cis interaction constant nearly to zero, but only reduces the trans interaction 


* Note added in proof. See also a recent paper (Torkington 1950 d) on the inclusion of 
van der Waals forces in molecular potential functions. 
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constant by approximately 20°. Hence there is probably a positive constant for 
trans interaction between the two ends of the molecule, but possibly no cis inter- 
action, and a negative constant for interaction between adjacent C—H bonds. 
The latter is to be expected (of that sign) from the changes in hybridization which 
occur during vibrations of the CH, groups (Heath and Linnett 1948 *) if these 
are of more importance than the repulsions between the hydrogen atoms, which 
would lead to a positive constant. The trans interaction may be related to 
conjugation. The positive d,., for Aroc_¢.Aro_y interaction, implies H... C 
repulsion. The negative interaction constant in the B;, factor, and therefore, 
presumably, also in the A, factor, can be explained if the positive constant for 
Aro_a-A8coq interaction (see below) is greater than that for Ar¢_y.A@oq, 
interaction, where the two angles are adjacent to the same bond; a positive sign 
is to be expected for such interactions, i.e. in triatomic systems, unless attraction 
between the two non-bonded atoms is important. The negative constant for 
Arco_c- Acq, interaction is equivalent to a positive constant for Arg_¢. Abc 
interaction (Torkington 1949f) ; this has the expected sign. Other interactions 
which seem definite are: Ar,,.(A@,;— A035), + 0:300; Ary. . (A@,4— A825), — 0-100; 
A65(A8o4— AOag), +9-09 (normalized values from the B,, and B,, factors) and 
A6,3.A64g, +0:0037. This last could be accounted for by H, ... Hy, but not by 
Fu. -C, repulsion. 

Ina recent paper (Linnett, Heath and Wheatley 1949) the vibration frequencies 
of C,H, and C,D, are accounted for by an orbital valency force field (Heath and 
Linnett 1948) with eight force constants and one van der Waals repulsion constant ; 
there are no interactions which are not derived from this latter. The value obtained 
for Ra_q was 7:16. ‘The constants can be derived by applying a quite simple 
transformation to the usual valence force constant matrix; the results of the 
present paper may therefore be used to obtain accurate values of orbital constants 
should these be required. ‘The reverse transformation applied to the orbital 
constants gives kg_~=8-97, which is somewhat different from the solution found 
here. ‘The interesting point of the analysis is the discovery of negative force 
constants for distortions of 7 orbitals which occur in the non-planar bending modes. 
These negative constants arise from the strengthening of the bonding caused by 
overlapping of the z orbitals during non-planar bending (cf. Wheatley and Linnett 
1949). This explains the apparently anomalous relation between electron- 
distribution and non-planar bending constant found by the author (Torkington 
1949a, 1950a). It is of interest that the Tchang anharmonicity constants in 
these modes are much smaller than the planar constants: B,,, 0-0171, By, 0-043; 
evidently the overlapping of the z orbitals reduces anharmonicity more in the B 


2 : ‘ 5 lu 
than in the B,, mode, corresponding to the larger negative B,,, force constant. 
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Note added in proof. After this paper had been completed, a new assignment 
for the ethylenes was published (Rank, Shull and Axford 1950, Arnett and 
Crawford 1950). ‘This is discussed in a recent note by the author (Torkington 
1950c). Any changes would affect only the B,, and B,,, factors. 


* But Coulson, Duchesne and Manneback (1948) appear to differ. 
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ABSTRACT. It is shown that the divergencies in classical electrodynamics may be 
removed by subtracting the self-energies of isolated electrons by means of a modification 
of the usual expression for the force density in an electromagnetic field. The results. 
obtained are equivalent to those of Dirac and Wentzel. 


S17 INTRODUCHRION 


T is well known that the classical theory of the point electron involves the 

| difficulty that the self-energy and the self-force of the electron diverge. 

A way out of this difficulty was first shown by Wentzel (1933, 1934) and 
Dirac (1938), and more recently other devices have been given for the same 
purpose (see, for instance, Fremberg 1947). However, we suggest another 
approach, which seems to be particularly simple. 

The present approach is based on the postulate that the self-energies of isolated 
electrons are physically meaningless, and our aim is to carry out a covariant 
subtraction of these self-energies in such a way that the other usual properties of 
the electrons are not affected. In particular we have to ensure that the amount of 
radiation emitted by an electron remains unchanged. It will now be shown that 
all this can be achieved very easily by introducing a ‘ counter-term ’ in the usual 
expression for the force density in an electromagnetic field, which will introduce a 
corresponding counter-term in the energy-momentum tensor. 


§2. BASIC EQUATIONS OF CLASSICAL ELECTRODYNAMICS 


Classical electrodynamics is based on the Maxwell equations 


OF _ 47 
Cue “cance 
Fw, OF OF 9 1 eee (L) 
Chae Mohan iene 
and the Lorentz equation 
1 
K,= me Terra efelens A (2) 


where J , 18 the current density, K,, is the force density, and F,,, is the electro- 
magnetic field tensor such that 


(Fos, #31, P12) =H, (Fin, Fas; Fi) =tiee ae (3) 


* State Scholar of the Government of India. 
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Ifn electrons are present in the field under consideration, then we have 


1 eee a3 A eae a a eR eS (4) 
J, pe Ri (5) 

: : dz, 
with jM=e [Se fe 8( 31 — x1) 8(24 — x) 8(z4 — x5) S(29— x9) dt, ...... (6) 


where /,,,,, -x is the external field, f°... is the retarded field of the rth electron, and 
j is the current density at the point x, due to the rth electron, whose coordinates 
are given by 2,. 

Using (1), we may put (2) in the form 


0 1p 
IS = Fe ee al Wied eer MeN aiexte sts (7) 
where 7’, is the energy-momentum tensor given by 
1 

Lo Ge (PuBin + 28 Pie"). eer. (8) 

In particular, (8) gives for the energy density w and the Poynting vector S 
1 
u=Tyy= 9 (E*+ H?), ease (9) 
SiGe rere) PE) ee oh a (10) 
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Having stated the basic results in classical electrodynamics, we now consider 
the elimination of divergencies. 


§3. ELIMINATION OF DIVERGENCIES 
As already pointed out, our aim is to carry out a subtraction of the self-energies 
of isolated electrons from the existing scheme of classical electrodynamics, 
keeping the amount of radiation emitted by an electron unchanged. For this we 
first introduce the following definitions : 


ee n= <a 3( me +f? sav) 
Oa rea 5( Ae ret — fi? adv Ar J 


so that we may write ee ps ened earn = 26 | LE ccoe a (12) 


where f”’,ay is the advanced field of the rth electron and fm, and fi) aq Will be 
called the symmetric and the radiation fields respectively. As shown by Dirac 
(1938), at the world-line of an electron its symmetric field diverges quadratically, 


while the radiation field is finite and given by* 


DEO aa - 
jen rad — B ct (v,2, ey O01); eh elelens.« (1 3) 


where dots denote differentiations with respect to the proper time s of the electron 
andv,=dz,/ds. (Differentiations with respect to the ordinary time ¢ will always be 
shown explicitly.) 

* Our definition of the radiation field differs from that of Dirac by a factor 3. Also note that 


we are taking the space-time coordinates as %, X2, X3, Ict. 
4-2 
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We now modify the expression (2) for the force density as 


“(F ol Eft som phe eee: (14) 


which immediately leads to a new energy-momentum tensor, for the symmetric 
field satisfies the Maxwell equations 


pa 
ane oe 
eee: (15) 
Of ora ee of! ie sym 4 Phe sym 6) 
Ox, OnE Ox, 
‘Therefore we obtain from (14), (1) and (15), in the usual way, 
ee Lae 
K’ /.= Ait pelinea kite een eae (16) 
where the new energy-momentum tensor T”,,, is given by 
f EE ene! Brive eg eae aN oe 5A (17) 


t\? .vm being obtained from (8) on replacing F,,, Dich ee 

It can be easily verified that the tensor (17) has the desired properties. For let 
us consider a single isolated electron moving with a uniform velocity. In this case, 
the retarded and the advanced fields of the electron being equal, and there being no 
external field, the total field is equal to the symmetric field, and therefore, according 
to (17), T’,,, vanishes everywhere. On the other hand, (17) gives for the modified 


Poynting vector 
PEGE c 
Ss = 7 [ExH]-27 [Ee xX HO a oe (18) 
(F 
which keeps the amount of radiation emitted by an electron unchanged, because it 
is well known that an electron does not lose any radiation by virtue of its symmetric 


field. In fact, by reversing the positive sense of the time axis, we immediately 


find 
{ [Eom yn dt= = | (En ecH di 0 eter (19) 


It only remains to be shown that (14) leads to a finite self-force, and the energy 
and the momenta of the field are always finite. Since at the world-line of the rth 
electron we have 


1p fors Ar, 


equation (14) gives for the force density at the world-line of this electron 
B= * (Ey —fOomdi®. tenes (20) 
From (20) and (6) we find 


. : 3 e ae < 
mace = | Bi dom EA (B, —famm). 
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Multiplying both sides by dt/ds, we get 


é€ 
= ) ‘ 
mv,‘ ~e C v,”) CR rad a x Hes hag sa )e ROMEO (21) 
sr 
Using (13), and remembering that 
v,? a C U,v, ae Os; 
we obtain 
3 De. 2 e ; é 
see pee 2(r) — 1 ) 
mv,” 3 3 0" ale 3 © uv, U, ©) ois C vu, x fe ret ae thes op we 5 Oras (22) 
SAT 


which is just the equation of motion obtained by Dirac. 
Further, we notice that at the world-line of the 7 th electron we may put 


i See as Ln fT | iE Me ot ecshe) (23) 
where f{).,., diverges quadratically and F.,, q, is finite. Thus, according to (8), 


mw, ew, 
T,,, may be written as 
T= 1 oat Lpmis 1 nite terms, 


wy, Sym. 


. . (r) . : 
where ¢,,, nix, Which contains the cross-terms of f| and F’,, an, diverges only 


Ww, SYM. 


quadratically. Hence, at the world-line of the r th electron we have 


Let 
ww 
Since T7,, has only a singularity of the second order at the world-line of any electron, 


its space-integral is always finite. 


Site COLINS Mae eo Nl eG ater oe (24) 


“uy, Mix 


§4. CONCLUSION 


We add a few remarks about the advanced field. ‘Though the advanced field 
has been in use for a long time, its meaning does not seem to be very clear. In our 
opinion the advanced field cannot be a measurable quantity. It is, therefore, 
significant to note that the advanced field enters into the present treatment in such a 
way that it can never be measured. For, according to (14), the force exerted by an 
electron on other charges placed in its field is always given by its retarded field. 
It may also be pointed out that though the modified energy momentum tensor 
involves the advanced field, this does not matter at all. For (17) satisfies the 
usual mathematical requirements of an energy momentum tensor, and it can be 
used in the usual way to calculate the amount of radiation emitted by an electron. 

Finally, it should be noted that the success of the present treatment depends on 
the fact that at the world-line of the electron f,,, ,,q is finite and, sym does not have 
a singularity of an order greater than two. Since these results hold even for the 
meson field (Majumdar and Apte 1948), this treatment is equally applicable there. 
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ABSTRACT. The subtraction procedure of Dyson is modified in order to eliminate 
a certain difficulty in the renormalization programme, namely the so-called “ b ’ divergencies. 
The conclusions of Dyson concerning the finiteness of the renormalized S matrix are 
confirmed. 


quantum electrodynamics manifest themselves in unobservable mass and 

charge renormalization effects alone (Dyson 1949). ‘The demonstration of 
this behaviour rested, however, on a correct, but hitherto unproved, hypothesis 
concerning the behaviour of divergencies associated with self-energy Feynman 
graphs, the ‘b’ divergencies in Dyson terminology. By a slight change of 
Dyson’s renormalization procedure, it is possible to avoid this difficulty, and the 
proof of the finiteness of the renormalized matrix is complete. We shall use the 
notation of Dyson, and refer the reader to his paper for the meaning of the notation 
used. 

The difficulties referred to above arise because it is possible to regard a 
reducible self-energy part of a graph as constructed from irreducible components 
in many different ways, according as vertex parts are considered to be inserted at 
one vertex or the other of the irreducible second order graph, and each possible 
way in fact contributes separately to the resulting divergencies. 

A method of avoiding the ambiguity of construction of reducible self-energy 


parts will therefore be explained here, which makes use of the formal identity 
(Ward 1950) 


| has been shown to be very probable that the divergencies which arise in 


If we define p*=pA+p'(1—A) where p’ is the energy-momentum vector of a free 
electron then 


2u*(p) = —2n i aX p,, —p, )A,(p*, p’) 


where the effects of the contribution of the mass renormalization term to =* have 
been taken into account. 

A,( V4, p’, p*) will now contain the infinities arising from inserted self-energy and 
vertex parts in the V’ as well as the divergence associated with the V’ themselves. 
However, there is now no overlapping of divergent parts, and the finite parts may 
be separated out unambiguously in sequence. 

An analogous procedure may be adopted for the reduction of photon S.E. 
parts. 

It is here necessary to define a new linearly divergent operator A_(t,,¢,) by the 
equation 

A. (ty, t2)= 2X A(X, ty, ty). 
all X? 
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X’is a graph obtained from a completely unreduced photon S.E. part by insertion 
of an external photon line in an electron line carrying the external photon variable 
of the S.E. part, and A(X’, t,, t,) the corresponding operator. + 


Then II* = —2n | dds,A, (At, At).t 


Again, A,(Az, At) need only be summed over the irreducible X’, the overlapping 
divergencies have disappeared, and the divergencies may be removed in sequence 
in the usual way. 

The resulting finite operators S,,, Dy,, [,,, obtained by this process of successive 
approximation are defined by the way in which they reproduce themselves when 
substituted into the integral equations of the theory, together with infinite parts. 
‘This is represented by the equations 


Sp, =Sp—2nSy [ AN—P.MPSPI- MPP) env (1) 
1 aA, 
R= D2—2,Ds | a, {4,08 At) —A,(0, 0) — (At) (sas), pe A 
(2) 
La Vet A Py Po) =A AP’, 2') veeeee(3) 


On the other hand the divergent operators [’,’, Sp’, Dy’ reproduce themselves 
identically when substituted into the integral equations, as is described by the 
equations 


SAC Cs a ee (4) 
De Daal Dee we lore (ater 
Deo ee Pee eh Seg (6) 


all irreducible V? 


fel , —L 
Dy — —2n | 9 (Pu Py VALS?) 
Il* = —27 iE ANE NE bo) 


It is now easy to prove that Sp =Zoh ri (41) 
Dy’ =Z3 Dy (1) C 
D/ =2Z,7T 1 (é) 


where e, = Z,1 Z, Z,!% e, provided that Z,, Z, and Z, are suitably chosen constants. 
This is seen by substituting these expressions for S,,’, Dy’, and [,’ into the set of 
equations B. 


Notes added in proof. = 

+ This prescription for the necessary function Aj(t, t) is correct only for photon S.E. parts having 
a single electron loop carrying the photon variable. This a ean pointed out to the author by 
P. T. Matthews. Avy is really to be defined by the equation— Eee = An, which is to be regarded 
as an implicit definition of Ay. It will be noted that Ay may occur once in the integral equations 
of the theory, since the definition of Ay implies that it may have to be inserted at one of the 
irreducible graphs defining Ay. | Hence one has to define a new divergent function Wu=2itu + et 
and a corresponding finite function Wy1=2itn.+Ayi. It is then straightforward to prove that 
Wu=Z;Wya(e,) with the value of Z; given in the text. 

+ The photon self-energy I1*(0) has been put equal to zero. 


56 F. C. Ward 


By counting up the number of vertices and lines of the irreducible graphs 
which define A,, and A,, one finds that 


E*(eq) = — 2m Zi |) AND, — Py!) A AP Ps) 


IG) == 2nZ ae | adt,, A, (At, At, ¢1) 


Bh 


with the Sp,, Dy,, I’, functions substituted into A,, and A... 
Equations B therefore become : 


Z, Sy (€,) = Sp — 27S y(p,— Py Zr 1Ze |, DA, (D’, P’)ArSy (ex) «----- (4) 
Za DiC) aD p= 2D at Ze ; OBA) Dy) vee (57) aBe 


Za (C1) =V pe het) 0 ee ee ee eee (6’) 


These have to be identical with equations A. Equations (3) and (6’) are the same 
if Z,y,=y,—A,(p’, p’). Substituting this result into (1), one finds that (1) and (4) 
are the same if Z,=Z, and finally (2) and (5) are the same if Z,=1+c where 


By if arse) 5) _ (Note A, (0, 0) =0 by invariance.) 
Hence the relations C are proved. The finiteness of the renormalized S matrix 
now follows in the manner indicated by Dyson. 

There may of course be doubt about the validity of the formal methods of 
manipulating infinite quantities used above. The author believes that it is 
possible to construct a consistent and unambiguous calculus in this way, as long 
as sufficient notice is taken of the requirements of Lorentz and gauge invariance at 
each stage. It is to be hoped that future developments will enable us to avoid 
dealing with unnormalized charges and masses, and the divergencies that appear 
here will then find their rightful place as abstractions from physical reality, 
namely the charges and masses of non-interacting particles. 
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ABSTRACT. ‘This paper contains a new derivation of the equations of motion of bodies 
moving slowly in their (weak) gravitational field, up to terms of the order (v/c)?.__ In contrast 
to the method developed previously by Einstein and co-workers, in which only the field 
outside the bodies has been considered, in the present paper a new method is used in which 
the field outside as well as inside the bodies, and the internal structure of the bodies, is taken 
into consideration. ‘The equations of motion derived by this new method are identical with 
those obtained by Einstein and co-workers, but the new method is essentially simpler and 
the necessary calculations much less lengthy. 


Soil 


T was proved by Einstein and Grommer (1927) that in general relativity the 
field equations contain the equations of motion of bodies inside the gravitational 


field they produce. Subsequently the problem of deriving these equations of 


motion from the field equations has been treated in detail, especially in the last 12. 
years. When one tries to survey the work done in these years one sees that two 
distinct methods have been used. The first is the method developed by 
Einstein and co-workers (1938, 1940, 1949), its characteristic feature being that 
ull the calculations refer to the space outside the bodies, for the interior of which 
only the assumption of spherical symmetry is used. ‘The second method has. 
seen introduced by Fock (1939) and, in contrast to the previous one, it also 
-onsiders the gravitational field inside the material bodies. 

Though far more developed than the second, the method of Einstein and co- 
vorkers is very restricted from the physical point of view, since it is only applicable 
o bodies with spherical symmetry. That this is a real restriction one can see 
mmediately by considering the motion of stellar bodies. When we treat this. 
problem according to general relativity, the first approximation confirms the results 
»f the Newtonian theory of gravitation, while the second approximation leads to the 
very small) relativistic effects. Now in the Newtonian theory the oblateness 
f, for example, a planet will cause a rotation of the line of apsides of a satellite. 
[his effect is usually much larger than the similar relativistic effect (i.e. the 
elativistic rotation of the line of apsides, which is independent of the oblateness 
f the central body). Thus in this case it will be indispensable to consider the 
lateness of the central body. But any departure from spherical symmetry can 
nly be treated by considering the internai structure of the bodies in detail, which 
as been excluded from the very beginning in the method of Einstein and 
o-workers. 

Besides the astronomical problem one might hope to use the equations of 
notion for discussing the motion of elementary particles in a gravitational field, 
nd this seems to be the ultimate aim of Einstein and co-workers. But the only 
ase of practical interest is the motion of photons in the gravitational field of a star ; 
nd for this case the method of these authors is entirely inadequate, since the 
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photon moves with the velocity v =c, while the method is based on the assumption 
v<c. A further argument might be mentioned here. The photon—and all 
other elementary particles—has a spin and therefore it is doubtful whether the 
use of a model having a spherical symmetry is justified even for a calculation 
in first approximation. 

In contrast to this Fock considers from the very beginning the astronomical 
problem only, and he also includes in the discussion the internal structure of the 
bodies. ‘This has the immediate consequence that the physical meaning of his 
calculations can be followed in detail. Additionally there is no longer any difficulty 
(except the more extensive calculations needed) in discussing bodies without 
spherical symmetry. However Fock has not reached the full efficiency of his 
method because he has insisted upon using the field equations alone. ‘Thus he 
finds the equations of motion to the first approximation as the integrability 
condition of the field equations in the second approximation by actually integrating 
them in this approximation, that is, wnen—as we shall see in the following—he has 
calculated nearly all quantities necessary for the determination of the equations of 
motion 1n second approximation. 


§ 2. 


We shall show in the following that one can give to Fock’s method a much more 
compact and effective form. We mention in anticipation that this will be possible 
when we use, together with the field equations, the well-known relation 


of,” vo O28 v6 
aur OO Oy 


thus allowing the simultaneous calculation of g,, and @,’. In the following calcu- 
lations this relation will play a role just as fundamental as that of the field equations. 
For this reason we shall give it a name (which will be justified later): we shall call 
it the dynamical equation of the theory of gravitation. 

It is not difficult to see why the dynamical equation has not been used in this 
problem up to now. Orthodox general relativity considers g,,, as the fundamental 
quantity, and this determines everything else. ‘The field equations are considered 
to be well defined only for the empty space, where they are 


= (iota Raa (2.1) 


Rye 0.0 a Rie ag eer (2.2a) 
In regions of space containing matter, i.e. where R,,,40, we put 
Sa fa goat a Ks etaleyolans (2.2) 


and this is then the equation of definition of 7.,,, which consequently appears as an 
auxiliary quantity. Under these conditions the dynamical equation becomes a 
mere mathematical identity and this is the reason why one has tried to derive the 
equations of motion from the field equations only. In its extreme form this point 
of view has been accepted in the papers of Einstein and co-workers who refuse to 
allow any considerations on 7”, and insist on deriving the equations of motion from 
the field equations for the empty space surrounding the material bodies (which they 
characteristically call singularities of the gravitational field). ‘To a certain degree 
this point of view has also influenced Fock, who, at first, considers T,,, merely to 
find out the general structure of g,,; then he calenlane the g,, independently of the 
Tf 


1» and only later does he determine T,,, from the already calculated g,,. 
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An entirely different situation arises in the more moderate interpretation of 
general relativity as a theory of gravitation in flat space, which has been suggested 
by Rosen (1940) and Papapetrou (1948). In this interpretation g,, is only the 
otential of gravitation (and no longer the metrical fundamental tensor) and 
consequently 7", will have at least an equally primary physical meaning as that of 
7w- che field equations will now be (2.2), reducing to the form (2.2a) in the 
“ase of empty space. Further, we must now consider equation (2.1) as an 
-quation having a primary physical meaning and not as a mathematical identity 
nly. A first argument for this point of view is the following. If we assume for T,,, 
he form corresponding to incoherent matter, 

CF =purw,, 
where uw“ =dx"/ds is the four-dimensional velocity of the element of matter, then 
>quation (2.1) leads to the conclusion that the orbits of these elements of matter are 
he geodesics of the field g,, (de Donder 1926). But this is a statement of such 
hysical importance that one cannot understand why it should be derived from a 
nathematical identity without physical meaning. We shall show in the following 
hat the use of the field equations (2.2) together with the dynamical equation (2.1) 
or the simultaneous calculation of g,, and T,,, will provide us with an extremely 
fective method for determining the equations of motion; this result may be 
-onsidered as a new argument in favour of the moderate interpretation of general 
-elativity. 
§ 3. 

In this section we shall summarize the physical and mathematical assumptions 
on. which the following calculations will be based. 

There will be two fundamental physical assumptions. We shall consider a 
system of bodies whose velocities are small compared with c and assume that the 
sravitational field is weak; i.e. g,,, will differ very little from the values 7, of the 
netrical tensor of special relativity. The difference g,,—7,, will be written as a 
ower series of asmall parameter. Following Fock we formally take as parameter 
he quantity 1/c. 

‘The next physical assumption which we shall use in the following is that the 
nutual distances of the bodies are large compared with their dimensions; i.e. if 
he distances are of the order of magnitude of a length L and the dimensions of /, 


FESS OL 0k pet ee (3.1) 


[his is an assumption which we introduce in order to avoid complicated calcu- 
ations. The reason is obvious. Stars consist of matter in the fluid state (this 
yeing a suitable assumption even for planets, though they have a thin solid crust) 
ind consequently, if the condition (3.1) was not fulfilled, the form of each star 
would depend on the relative positions of all the others; the problem of their 
notion would, then, become a very difficult one even in the Newtonian theory of 
yravitation. But when the condition (3.1) is fulfilled, the form of each star will 
Jepend only on the conditions existing on this star; and in the stationary case it 
will be a spheroidal one, the oblateness of the spheroid being due to the proper 
‘otation of the star. 

For a further simplification of the calculations, we shall treat the case of 
pherical bodies only. (The only difficulty in treating spheroidal bodies 
s that the calculations are more lengthy.) This will be true if the bodies have 
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no proper rotations; consequently the last simplification must be combined 
with the assumption that the motion of each body is a pure translation. . 

Mathematically, we first have to choose the independent variables (coordinates) 
and the unknown functions which we shall have to calculate. Following Fock 
we take as independent variables the space-coordinates «' (¢=1, 2, 3) and the time 
x°=t. As unknown functions we shall take the quantities g” and @,". Generally 
latin indices will take the values 1, 2, 3 only, while greek indices will take the values. 
1,2,3 and0. Asubscript comma will indicate ordinary differentiation : 


alee pS, 
fi= Gy Lo=oq0 = Oe 


We must now introduce another assumption, which according to Einstein and 
co-workers is merely a matter of mathematical convenience, but which in our 
opinion might have a deep physical meaning. We must choose the so-called 
coordinate condition, which from the mathematical point of view is necessary in 
order to restrict the variety of forms which a given solution of the field equations. 
can take (because these are tensor equations). But the physical aspect is very much 
different: Our problem is to derive equations of motion from general relativity, 
and to express them in terms of space and time of Newtonian mechanics. 
Obviously not every frame will be suitable for this purpose; the coordinate 
condition will have to exclude the unsuitable frames, or, in other words, to 
determine what we might call the inertial frames in general relativity. One must 
expect that the choice of this condition will influence the equations of motion, 
at least in the higher approximations. * 

In this paper we shall accept the coordinate condition used by Fock: 


W00.9 coe eae (3.2) 


which has several remarkable properties (see Papapetrou 1948). It is worth 
stressing that by using (3.2) we shall in the second approximation arrive at the 
same equations of motion as Einstein and co-workers did, though these authors. 
have used a different coordinate condition. ‘The reason for this coincidence and 
the deeper meaning of the coordinate condition should be discussed in more detail, 
since there seems to be an important problem behind them. If the equations of 
motion in second approximation really depend on the coordinate condition, then 
the observational data will give support to some special coordinate condition; but 
such a result would be a decisive argument in favour of the interpretation of general 
relativity as a theory of gravitation in flat space. 


§ 4. 

Before starting our calculations we still have to express the field equations and 
the dynamical equation in terms of g“” and @,"; additionally in the field equations 
we have to consider the condition (3.2). We shall not repeat here the calculations. 
necessary for the transformation of the field equations, but we take the final 

* Ina recent paper Einstein and Infeld (1949) have tried to prove that in the second approximation 


the equations of motion are independent of the choice of the cocrdinate condition. But their proof 


extends only to a special class of coordinate conditions; the initial condition being given by their 
equation (9.8), 


Y00,0 Y0n, n— 0, Ymn,n— 0, 
they keep it unchanged in the first approximation, allowing a change of this condition only in the 


second approximation. However, the important question is what will happen to the equations of 
motion in the second approximation if we change the coordinate condition for the first-order terms.. 
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ormulae from Fock.* With a slight difference in the notation g“” will in second 
ipproximation be of the form 


el 40 — 4,500 
ee Gm rpt oint 
. 4U% — 4,80 
ge ae oe eee eo Spee (4k) 
vale | 
ike 
q a oo 6 | 


The field equations up to terms seme nen to 1/c3 are 


KE j= — 5 ASH 4 U% 43-5 Eades ORL 


2 2 2 7 8 
-«T,o = - AU 4+ = AS00 — = U 9, — 5 UU , — = UA U0, | 
Z 2 2 6 
~KT,* = : AUSk 4 = ASOk — 2 U% oy + = UU ,, | 
) 
2 (eee ee : UA UK, (4.2) 
- Bp = — 2 ave, | 
c 
| 


withx =87G/c?. (Gis Newton’s gravitational constant.) These are the only 
results which we shall take over from Fock; all other calculations necessary for the 
derivation of the equations of motion will be given in this paper. 

It will be more convenient to express the dynamical equation in terms of g” and 
t,,.__Asimple transformation of ee 1) gives 


! 1 : 
©’, =320,"a" ge? a 2 Ct,’ A Ce) {vV/( Sak (4.3) 


{v(-2)}? . 
The quantities g,,, and \/ —g can be found by straightforward calculation from (4.1) 
and are 


4gu ] 
Boo =O? — t= | 
4U% = 480! | 
i See eee ee oP et ee ict Sen 4.4 
No; Cc C3 ) ( ) 
4 4U% 4500450 jpee 
fin = — Bike (1+ aot 7 | 
J 


ee =o oo + (5 S#— (0%) i eee (4.5) 


[t is convenient for the following to write equation (4.3) separately for »=0 and 


Lt: 

1 
7% T,,° =1€° oo : 4 (ied 
ee dee Uryepp oc va 


1 
aot Legere =e rs eee ey Ear) Pee (4.7) 
sik i,o-2%,9 (/( fae Go, 2 a/( 2) {V/( )} 


EER) eevee as (4.6) 


* For those readers who are unable to see Fock’s paper we summarize the main steps of these 
calculations in Appendix I. 
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Following Fock we shall denote the different bodies by a, b,. . . . (at the same time 
avoiding the use of these letters as tensor indices). If we integrate (4.7) over the 
volume of the body a, the term @;” ,, drops out (since €,” =0 outside the body) and 
we get 

1 


carte 1 
fl ed =| 4eege —— ——  aV |  ey 
a), AV =| Ate" ppp tent — | BE! Tay tv(-ah, 


We shall see in the next section an essential difference in the usefulness of (4.7) and. 
(4.8). 

We stress that, contrary to the method of Einstein and co-workers, in which the 
field equations are separated into equations for each approximation step—thus. 
making conditions of integrability necessary for each step—in our calculations for 
the solution of the equations in the /th approximation we shall take all terms of 
order /'</. Noconditions of integrability will be needed since the equations will 
be integrated directly. 

28, 

We shall first derive the equations of motion to the first approximation. This 
will show the efficiency of our method, and at the same time it will give us some 
quantities which we shall use later. We start by writing down the leading terms 
in the expressions of @, which we take to be of the same form as the components 
€,” of a perfect fluid in special relativity : ) 


f,° = pe, 
T,* =pev*, 
a ec (5.1) 


; 1 
CF = -£ DOr - 3" ; 


v' is the (3-dimensional) velocity of the fluid at the point which we are considering. 
Since our system consists of a number of discrete bodies, the mass density p will 
be of the form 

P= = Pas 0 Se eee (5.2) 


with p,A0 only inside the body a. Because of the spherical symmetry of the 
body a the density p, at the point r (coordinates x’) will be a function of the 
distance r,, between this point and the centre a (coordinates a‘) of the body a: 


Pata) (74 =|r —al). a aeoretiee (523) 


If /,, denotes the radius of the body a, then p, =0 forr,>1,. 
Since the motion of each body is a pure translation, the velocity of any element 
of the body a will be equal to the velocity of its centre: 


inside a a (5.4) 


~ 


It follows that vi =O (9 =a) a ee eee 6.5) 
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lhe pressure p will again be of the form 
PE i ee Weak ties (5.6) 
he partial pressure p, will be determined from the dynamical equation (4.7) 


ogether with the equations of motion. Using (5.2), (5.4) and (5.6) we can put 
5.1) into the form 


C,° = Upals ] 
a 
T,* i Up,qca*, 
a 


Laura ia 
Wiley tye ee ee a (3.7) 
ac 


ie Tae 1 
ieee: ye ink _§ ky — 
Th = —E pall HKD pa 


We now write down equation (4.6) for a point inside the body a; retaining 
only terms of order 1/c we obtain 


(PRCe irs ps0). = US ae, ee ie Beas (5.8) 
Since for any function f(x’) it is 
Oe de moj 1 dy. 
k pe HS Naked aps eae 
Fv +fo= ae dt at dt’ SO OOO (5.9) 
we can write (5.8) in the form 2: =H 5 WS ley, (5.10) 


i.e. in this approximation the mass density p, is constant in time. We then apply 
equation (4.7), now retaining only terms of order 1/c (there are no terms of order ¢ 
in this equation). ‘The result is 


dedtrey. | 1 
== 5 (pal!) = = Pai=— = Pa seers (5.11) 


We see that this relation contains the quantity U®. ‘This quantity together with 
U°* (which is not necessary for the present calculation, but will be needed later) 
can be calculated from the first two of the equations (4.2). “Taking these equations 
with the first order terms only and considering (5.1) we get 


AU = —4nGp, 
Fanner (5.12) 
AU% = —4nGpv*. 
Because of (5.2) and (5.4) the solutions of these equations are 
C27, 
ea eee eee (5:13) 
a 
with uw, fulfilling the equation Nava AmGp, ua wis (ox14) 


and being therefore the Newtonian potential of the body a, Using (5.10) and 
(5.13) we can write equation (5.11) in the form 


Pad! + Pai =Pal= Us,i t+ Ya,i)» eoccee (o.1'5) 
6 
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the notation © meaning summation over all bodies except the body a. (But X will 


b a 
mean summation over all bodies. This convention will be used throughout this 
paper.) 

Equation (5.15) can be split in two if we notice that u, ; is very nearly constant 
inside the body a because of the assumption (3.1). ‘These two equations will be 


CY (5.16) 


Pa,i=Pata,is eceee (5.17) 


the suffix a outside a bracket will denote the value of the quantity inside the bracket 
at the centre of the body a: 


(fio VaULeT OLS, “Ath anw 9 ba) enters (5.18) 
We see that equations (5.16) are the equations of motion to the first approximation 
—-identical with the corresponding Newtonian equations—while (5.17) is a 
condition determining p,. Since p, will be a function of r,, we can write (5.17) 
also in the form 


DoF hi ae ee (5.19) 
the accent meaning derivation with respect to 7,: 
’ ee asp 
Lira) = ae ere (5.20) 


Equation (5.19) is again identical with the corresponding Newtonian equation. 
Another way of arriving at the equations of motion (5.16) is to use equation 
(4.8), or, equivalently, to integrate (5.15) over the volume of the body a. ‘The 
term p,,; makes no contribution because p, =0 outside the body, and again the 
integral of p,u,,; will be zero since u, ;=u,’(x'—a')/r,. We therefore find 


a 


with ma= | pad ; ee (5.21) 


i.e. the equations of motion (5.16). We now see the difference between equations 
(4.7) and (4.8). ‘The latter gives the equations of motion alone, while the former 
gives an additional condition determining @," in the approximation which we 
are considering. We shall make use of this remark in the following calculations. 
In the second approximation we shall only be interested in the equations of motion 
and consequently it will be sufficient to use equation (4.8). (We had to make use of 
the full equation (4.7) in the first approximation, since the quantity p, will enter 
in the second order terms of €)".) 


§ 6. 


For the derivation of the equations of motion to a second approximation we 
shall have to follow a step further the procedure outlined in the previous section: 
calculate {," and gq” and then apply equation (4.8). Not all the second order terms 
in ©," and g’” will be needed in this calculation (as in the calculation of the previous 
section, where only some of the first order terms were used). What is really 
necessary will be seen when we write equation (4.8) up to terms of order Lics: 
Anticipating equation (8.1) we state here that the following quantities will be 
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needed: (i) from €,” the quantities T°, T," and T;* up to terms of order 1/c, and 
the quantity ° up to terms of order 1/c?; (ii) from g” the quantities U and 
U% already determined in the previous section and additionally the quantities S” 
and S°, 

We start with the calculation of €,. For @,° we provisionally keep the 
previous expression 


eee 1) iene as (6.1) 


including in p terms of first and second order. In @,” we introduce an additional 
second order term depending on the pressure p, identical with what we should have 
for a moving perfect fluid in special relativity : 


1 
T,* =pcv* + - DOR ca! ON et (6.2) 


‘The leading terms in @,° and {;* will be the same as in (5.1) but with the first order 
term of p instead of the total p. The second order terms in @,° will be deter- 
mined later. 

We first derive an interesting relation which is valid for bodies undergoing 
arbitrary deformations. For this purpose we combine (4.6) and (4.7) after multi- 
plying the second by vw’. Taking in account (5.9) we find 


(VE, + Ey") .t (VTP + T°) .—v' ,TF-v' ,T;? 


1 dq Ae ae) 
=1, “92” Be “eae RL 4T,° BS wean 6.3 
SW CaFe de 8 7-9) dt ee 
Introducing in this relation the values of @,” and q” and retaining terms up to 
the order 1/c we get a iS ag 


We now consider the volume element SV containing a given mass element pdV. 
After a time interval dt this mass element will occupy a volume 6V +(d/dt) (6V)dt. 
‘The derivative (d/dt) (6 a can be calculated easily and is 
£6 Vj=cV, divv= MANE eens (6.5) 
If therefore we multiply (6.4) by 6V we find 
d Cf Bie bed 
oF ERS = ae” 5 =) + 3poV = (5 :) + + 5a ) +e ZOV)=0. 


According to this equation the quantity (d/dt)(pdV) is of a higher order and 
consequently we can put (d/dt) (p5V/) =0 in the second order terms. We thus get 


d v2 U% 
riaae —5at a) | +; Shs EAS iy Ove ee (6.6) 


a relation which is valid for arbitrarily large deformations, as for example those 
encountered in pulsating stars. 
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In our case we have only translational motion for which (d/dt)(6V)=0. 
Equation (6.6) becomes therefore, if we consider the interior of the body @ 


(with v,? =a'a’) 
dT fy _ m8, U7 
di| P@\- Da 7 2 - 


This relation can be written in the form 


Va" [00 
pa( Soi “| Sy eee (6.7) 


with dp,*/dt=9. For reasons which will become evident later it is more 
convenient to introduce another quantity p,* by the relation 


Ka ate va, 
aye (1 +3). Paros (6.8) ; 


Since the Newtonian potential ~, is also constant in time, du,/dt=0, we shall | 
have dp,*/dt=0. It follows from (6.7) and (6.8) 


oes A ANGI S97) 
pa=po*(I ese or +54): 


In all the following calculations we shall use the quantity p,+ only. We can 
therefore omit the cross and write instead of (6.1) and (6.2) 


riers WE Nale 
= ae = a 
Ey =E pee(1 okay eee 35) , 
ee eat fa tage eae (6.10) 
k yh SN pt. 2 ye eg Sa ilk 
TE EPaee (i oy 22 C2 at: 34) +3 pas ? 
with 
apa 
So 0) eee (6.11) 


In agreement with (5.7) and (5.10) we shall consider p, as the first order term in 
the expression for the mass density. ‘Thus equations (6.10) give ,° and T,* 
separated in their first and second order parts. 


Since ©," is needed up to terms of order 1/c only, it will have the same value 
as in (5.7): 


Tapa 
Gt = —¥2 pyaiak — = 5B ou. Sane ieee (6.12) 


Finally we have to calculate T,° up to terms of order 1/c3. This can be done 
by using (6.10) as follows. It is 


1 - 1 
C77 yas ele — TT aaa (FB Lo? + 0 Gig G9 + 9G, Cy* + 9° a, T,/). 
iy (=) VT) 


Using (4.1) and (4.4) we see that the second and fourth term in the bracket do 
not give any terms of order 1/c*; taking into account (6.10) and (4.5) we finally find 


il : 4 = Ass 
(Ua =e Na ,' — <u pglt,(a'— B'). eos e278 (6.13) 
; ‘ 


2 3 
Oi 
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This formula displays a very interesting feature, expressed by its last term. 
Inside the body a the momentum density €° has terms proportional to the 
relative velocity of the body a against each one of the other bodies. 

We notice that the values of €,” given by (6.10), (6.12) and (6.13) are identical 
with those following from the values T”” given by Fock. 


Sls 


The next step is the calculation of S” and S®. From the last of equations (4.2) 
we findt 


Ke? 
A Sie “7 CO El ey ate ae (7.1) 


/ is to be calculated from (6.12), while for U as well as for U®* we have to retain 
the values (5.13). By taking into account the relation 


Ce = ACU ee ON ates (7:2) 
we can write (7.1) in the form 
AS" =A[7(U)? +2 u,0,7] —407Gd (3p, —4p ata) +20Gip,Uu,.  ...... (7.3) 
a a a b 
The solution of this equation will be 
= (Ol eee De (aes We ee ae (7.4) 
with 
INS Jo te G (SP, sulla) |e eee (7.5) 
and 
RV aD pe ys. 0 = ee Bes (7.6) 
b 


We first notice that the quantity s, has the property that it is a function of 7, 
which vanishes outside the body a: 


Sea) OUCSIGes te: DOGVAOst ea. geen Aarne ee (7.7) 


This is so because, as shown in Appendix II, the following relation holds: 


A | patdV = 3 { EL Va ats gil ho ae (7.8) 


Because of this property s, will play no part in the final calculations, as we shall 
see in the next section. "his is also the reason for the introduction of p,* instead 
of p,* in the previous section: if we retained p,* we should now have a more 
complicated function s,,. 


+ One might wonder whether, after having determined the quantities g/” from "by means of 
the field equations, the coordinate condition (3.2) would be satisfied automatically. It is easy to 
see that this will be true in our case. Combining equation (4.2) and taking into account the values 
T" already determined, we find (neglecting terms of higher order) 


A(U%, 9 +S x)=0, 
al U0 9+ U% E+ 4 (S 9+ S%, » | =0, 


In contrast to the calculations of Einstein and co-workers, the equation AA =0 leads to the conclusion 
A=0, since we are considering only continuously distributed sources (no point singularities). 
Therefore the above equations, together with (4.1), give &'%,.==0 (up to terms of order Vics) 
g°~ .=0 (up to terms of order 1/c’). 

57-2 
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In order to calculate W,, we develop u, in power series around the point q;: 
Up =(Up)o + (Uy, nal%* — a‘) + Seg cubes uses. (72) 


the omitted terms being of no importance in this approximation because of the 
assumption (3.1). We now get 


W = — 4g D (Uy) g— F(a" — a"), (Up nar see (7.10) 
b b 
the new function w,=w,(r,) fulfilling the equation 
Wy + = Wa = —47Gpy. 


The solution of this equation which fulfils the boundary condition is 


4nG/1 to A e 
=——(- 5 eee Ze 
Wy 3 (= I, Pala dr, a ie Pala ara) ( ) 
We finally write equation (7.4) in the form 
S#=1(U)? + [u,v,2 — buy 2 (Uy) +Sq—4(x* —a*)w, % (pj \els ena (7.12) 
a 


For S® we find from the first equation (4.2) 
3 
AS = — > To — FAU™ + U og + ; CO tt CNL ES a eee (7.13) 


Taking in account (6.10), (5.13), (5.14) and (7.2) we write this equation in the 
form 


AS® = —27G = Pal? — s Uy) + - MN (OY Se 1a cons? en BS (7.14) 
‘The solution of this equation will be, if we remember (7.6) and (7.10), 
one = U%)? + » [5UgVy? — ta (Up)a— (x — Ores (Ugly lal, sees (7:19) 
with AP [Susser (7.16) 
We still have to determine F,. A straightforward calculation gives 
eta ~ (=) Ge a’\(ak —ak)aiak — ne (x! — a')di +. <2 Aijie pee (7.17) 


‘This suggests putting 
Ee a (x! rs a')(x* ‘e a’)a'a*f(74) =e (a! = a’)i'b (Ta) aN aa'V4(1q)- SereiTevelie (7. 1 8) 


For the functions /,, 6, and v, we then get the equations 


Ta\la 
4 u’ 
| fail (nine sear vs Ae! 1 
Gl tad LM (7.19) 
Vo Vo 
2 1 
" fi es Gee, 
© G See ' Da ~~ 
a a 
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The solutions of these equations are 


u Sua 
= _& 2 
if a is | Ugh a ara, 
Ta (0) 


Laie , 
ba mk COs G Ugla ary. 
Va 0 


In addition we shall need the value of F, outside and at large distances from the 
body a. One can easily see that for this purpose in the integrals of (7.20) it is 
sufficient to use the value of u, outside the body a, even for the part of this 
integral referring to the interior of the body. Since this value is 


Fa CS i ean Sr (7.21) 
Va 
we find, after some simple calculations, 
aT. Amel be 
F,=4Gm, («ava Biiaiee a a] fot 7 le Vee 2a (7.22) 


In the final calculations we shall have to use the quantity S” +S —2(U%)?, 
We find from (7.12) and (7.15) 


SY +S — 2(U)? =¥ [8u,0,2 — ugg +5 + Fy—(xi—ai)w gg]... (7.23) 
with the additional notation 
(7.24) 


Ug= = (Ura, Uai=U (Ui 
b b 


§ 8. 
The final step for the derivation of the equations of motion in second 
| approximation is the application of equation (4.8), now up to terms of order 1/c’. 
Taking in account (4.1), (4.4) and (4.5) we find 


d dV 4 1 
at GdV = I. @| (tl ©) Ue, et: ao" Coe vos ToS? 1 S00 = 2( yay | M4 


where in the last two terms of the right-hand side we have to take T)* and T,° 
| in first order only (since the contribution of the second order terms will be of 


order 1/c®). ; 
The calculation of the integral on the left-hand side is straightforward. 


Considering the value of ,° given by (6.13) and (6.10) we find 


Oe AN bx... sep @ioe e, a 
od pe a (jf SS = d 
ik Gals c E BS oan ae Cc, 


with m, and ¢, defined by (5.21) and (7.8). | 
| The first integral on the right-hand side of (8.1) is, according to (6.10), (6.12): 


and .(5.13), 
dV 7 Pa Rite ye Ue Pa 

IE 7 CE a a = av|Pe(1 pen oe a) 

) 
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Since for any function f(r,) it will be 
| fuga =| fur’ 


while on the other hand we can in this approximation replace u, ; by (up, i)a, We 
get from (8.3) 


oy tee 
| ees ay Rut foe 


Using the power series (7.9) we find 


; : the Craig: 
| AV pala, iy = (Us, wa | dV pq, —# (x'—ai(x*—a¥). ...... (8.6) 


‘This can be simplified if we take into account the following relation holding for 
an arbitrary function f(r,): 


{ flrg)(2e! — a’)( 0" — a®) dV = 38; | fred. ees oa 


We then get 
d 
| S@i-aum, 


Mg Wat do Ra! ( ’ 
— — me(1 = 9) See 4) = (Uy ia a 58 = (Mo,i)a j . Pata r,4V. 


The second integral on the right-hand side of (8.1) gives immediately, if we 
consider (6.10) and (5.13), 


| aE AE KUM, = s mn fa SO * (tips. bas Ree (8.9) 
a b 
The last integral in (8.1) is, according to (6.10) and (7.23), 
dV 
~ | ayses $2 U%P| = — 5 | pudV a, Feat Us) 


Ss 3 Zee \" 
“ by (SU, — Uy) Sah Sy lg ee as 
b b 


culo = a" \ton | ae Url e108 yep peel eee ee (8.10) 


‘The first term in the bracket on the ee side of this equation gives zero 
because of (8.4). In the second we have to put (wu, ;), instead of uy ;: 


| padV 2 Uy, (20? — Up) =m, > (Up, i)a($%s°— Uy). +s. (8.11) 

a b b 

The integral of p,s,,; vanishes for the same reason as (8.4) and the integral of 
PaXv8p,i again vanishes because of (7.7). The term F,; gives, according to 


Appendix III, 
| pl dV = ted (8.12) 
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In the next term we have to put (F, ;),, instead of F,, ;: 
Pal Ky dV amgE (Fre vee (8.13) 
‘The next term gives, remembering (8.7), 
fe Paarl («* —a*)w,] dV =U, | Polat $0q'ta)dV. sears, (8.14) 


‘The last term in (8.10) makes no contribution in this approximation because 
of (3.1) since, according to (7.11), it is 

ELEPA aoa (Me saat ead Bis 
‘Thus we finally obtain 


-| * WY ©,0[St 45% 2(U)2] , 


Ma D My 
aos es asus. — Uae 208 0 73 2 (ras 
1 a 
Biz 7) Uai| Pal(@q ue 3, r,aV. she siiessss: (8.15) 


In order to find the equations of motion we have to introduce the results 
(8.2), (8.8), (8.9) and (8.15) into (8.1). We notice that there are three kinds of 
terms. Most of them have the factor m,, two have the factor e, and three are still 
in the form of integrals. ‘The terms with e, cancel, since because of (6.11) 
du,/dt =0 and therefore according to (7.8) de,/dt=0. ‘The terms in the form 
of integrals give the sum 


1 ~ 
eS eT V, 
yar (usa) Pata ‘ Pow, + We 7 akg 
which is also zero as shown in Appendix IV. We thus are left with the terms 


having the factor m, only. Remembering that dm q|at = 0 because of (6.11), 
and dividing by —m,/c we get the equations of motion in the form 


2 aa Os aa Uc 4 ; ; 
ata sata ae fery + 3a 5 (a) - gu bY uy)a — ub T (oa 
3 uv Nx 
=3 (tide + (52S — BE) EO de— Be EMH, 
Beveeue oU), 1 
+3 (,da(5 & - Gt) + al ae ce eeee (8 16) 


In the second order terms of (8.16) which contain second time derivatives 
we can introduce the first order value 


a => (uy, i)a: 


b 


We further write pee (Usa)jy oes ls Gane (8.17) 


denoting by a the sum over all bodies except a and 6. Similarly we write 


Up==(ue)o+ Wale tts (8.18) 
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Denoting by 7,, the distance between a and b, 


Tab = | a—b | , 
we get from (7.22) 


S My 6 “ 5 037 5 os cae Or ap 
(Poa ae G D (i b daioa*da! b daoa” iss eee 0 (8. 19) 
If we then take into account the relation 


Bray _ df _ (a'—bi)(a*—b%) 


daida® 7) Tap 
we find 
ener 1 1 (ile 
BE Me — BE Valea goa Cm saisak 


=— 326mm, - re ~ sa(=-) + tL Ba i)o(Us) a 
b 


» Oa’ OE be 
1 1 
— o> (H)o(Mo,a + Te — a") (Ue 1) o(Ma,i)ax 
er (8.20) 
Further, we find by direct calculation, 
Cx 
(tp, :)a=Gmp wails) ) 
d Pete Pe aay ho! 
a —b"/\(Uy, x)a = Gm,(a* —b ) ak Top aa Vie Was (8.21) 
5, (Ua) = ZS j (Hn) = AE (ty, 4a EBM Ur, i) 
p) 
And if we write 8 re of (7.24) in the form 
Oh =(Uy)q +2 (Ue)as 
we find 
4 4 4s: 
moo nll me Ae Say C2 a (Uy, ia(Mr)a a a 2 2 (Up, ia(Ucdax 26 COC (8.22) 


Introducing all this in (8.16) we get the equations of motion in their final form: 


ia Gm, 2 (— 4, Gis ahah + 3bhbE— Aakpe_4 CM _ Gm,\ 0 i 
See da'\ Fay ce ab da! 


Gy (gii-n_agigh + 3gige —ajipe) 2 (1 ses kpl rn 


.Gmm,~ 7 a (1 2.0/1 en 
ae 2c? Peis wel of Yue da' (=) -; Tae dai ic =| 


+ (bk = a*) = (-) - (—) | ee (8.23) 


If there would be two bodies only—the bodies a and b—then the double sum = = 


b 
is zero and the summation & on the right-hand side of (8.23) has only the one 


oy) 
term corresponding to 6. In this case equation (8.23) becomes identical with 
the equations of motion given by Einstein and co-workers (1938). 


EE —— 
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AP POND ex al 


After a series of transformations in which the following abbreviations are 
introduced 


1 
hi2cc= oe: OT po OT MO a llT p00 
2, (8 ie te yn a), 
i cet so) Weis ) 
Vi logy (—£)} a. Y= 29 63 
I gi? 
(Sra 
LS =, (Biol = He Seok Z Mt Sie °), 
{De eer a 


g being the determinant of g,,, (or g’”), one obtains the general formula 


lee 


v uy 1 00. v ee v 1 iy OT v y 
es Fe Se 2% Gio lee + ie CeCe eae. ey 


— DA DHy + Dy!) + $e (geT og + LY) 
In our calculations we assume (3.2), i.e. '=0 and consequently the last three 
terms on the right-hand side of this formula vanish. 
In order to find the field equations up to terms of order 1/c? one wants II 
up to terms of order 1/c* and all other II’ up to terms of order 1/c*. Considering 
(4.1) one gets by a straightforward calculation: 


Zz 
Weta sbi Ta OG) 
2 
[oe — [1.00 =— ei Oe ee Or) 
D, 
[]o-i* — = CUS tee) ree, 


re xo — 2 (U% .— Use) aes 


ct 
Tess 5 
the dots meaning terms of higher order. It follows that 
8 
11°10 = — m UU + ..., 


; 4 
Tél * = - U U0 + ..., 


nd: 8 
TT! f= “6 Oe en “a (U% ,— U% UM (+ woes 
yw 8 00 [00 
IT 2a eo oa ail Brae oerster 
Further, one finds 


i, So ees 


ah : 
yey = — me UA OM | ooee 
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Introducing these values into the expressions for R’”—4g’’R Fock obtains the 
result 

= 1 Z 

—$ (RO — $2°°R) = aa + Ag? 78 3 8 S00 — “a Bey Se 


- ; , 6 8 ; 
a Sarees =a" oot U7? ,U% , (Ue Ue 


Sg. 2 ! 
—F (R*—4g*R) a2 5, Ant ae ai YO [jo | — q5ikU™ U +.. 


From these equations by lowering the one index and considering (4.1), (4.4) 


and (4.5), we obtain our equation (4.2). 


AP PEN DEX Tt 


We have 3 { aor | Patqedr, =40 i d(par,2)—4n | 7p! dq: 
a 0 0 0 


The first integral vanishes because p, is finite at the centre of a and vanishes 
outside the body. ‘Taking into account (5.19) we find 


3 | padV = —40 | Tq°Palta 47g. 
a 0) 
And since according to (5.14) 


—47Gp, =U,_" + cae 
we ge : Ra = all Yo \ Ug Ug za r, 4 a 


oat ee Tere 
aes Syel\e ee oe 2, 72 
el ere le ie 
Again the first integral vanishes because u,’=0 for r,=0 and u,’~1/r,? for 
7,>#. Therefore 


3 | Pua = ie Te rer 


Similarly we find 
1 ee) 
= page ie Oe 
3) PatladV = xa). (1 Pty’ Ug) + = Bea plat 
‘The first integral vanishes a we finally get 
1 Tavis? 
Y aS ee ll WY Le) 
3 |, PatladV neue Yq Ug ao, 
which proves (7.8). 
ACP RON Dexa Lr 


In the integral | Pal’y,;4V only the second term of F,,, equation (7.18), will 


give a non-vanishing result. We find 


J Pali av =| Pat |(x* = a bal i av =| Pat" [9;"b, ctr * by'(xt— a’)(x* — a*)) dV, 
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And according to (8.7) 


| pakuidV =@ | paltat d4u'ra)AV. 
a a 
Now because of (7.20) 


ba 3 fu tA, t 
=S=_l ceo ae —( ono 
a rA x ala “la r 
et Ug 
and Oya 3 Pa gir 2 . 


‘Therefore, according to (7.8), 


| Pata dV = — a a 


a,t 3 a 


ACD AE ENED TX) Lev? 
We find from (7.11) 


, 4nG ("a 
De | 'Wenl RAe 
a 0 
i a es " 2 if 
and 3Wy t+ We 1q= 4nG| Palate — | (Ug" + —Ug )Ta@a 
"a "a Ta 


Sls | ; U(r Uy ste Uz) a Teale + Ug. 
a 
“Therefore 
| Paltta’ta + 300g +204! 7y)AV = | pa(2tta'Ta + Ua)dV. 
a 


a 


On the other hand 


( / 2 = wu" 2 , / 
| Pa2ta r4dV ae G | (Uy aaa lla )Ua rear 
a 40 Yq 


+00 1 ee) 1 - 00 ; 
a hee > UD \iee ees 2,, '2 neicaag PH 2, 2q . 
an GE | F AU Ta ) CG I, Vo Uy ary al Ya Ug Ya 


And since according to Appendix I 


| PallgdV ar = | Hee Pars, 


J 0 
we shall have | Palit ta + 300, + y't,)dV =O. 
a 
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ABSTRACT. Since the lower excited states of nuclei have very small widths (< lev.), 
resonant scattering of gamma-rays requires precise matching of the energy available from 
the gamma-ray with the energy necessary to excite the scattering nucleus. 

Resonant scattering should be observable if (1) the emitting and scattering nuclei are 
of identical type, (2) the gamma-transition goes to the ground state, and (3) the source 
and scatterer are given such a relative velocity that Doppler effect restores the energy lost 
by the gamma-ray to nuclear recoils. Thermal velocities of the emitting and scattering 
nuclei broaden and correspondingly weaken the resonant scattering peak, and the cross 
section at the optimum speed of 32E/A cm/sec. is 3-6 x 10-°(/T/E*)(A/T)*? cm?, where 
E and TI are the energy and intrinsic width of the excited state in electron volts, J the 
isotopic abundance of the resonantly scattering isotope, A its atomic weight and JT the 


absolute temperature. 
Preliminary experiments have been made with the 0-411 Mev. radiation from the 


nucleus !°8Hg, the source being carried by a high-speed rotor up to a speed of about 
7X 104 cm/sec. and the scatterer being liquid mercury (10% 1°8Hg). A small but apparently 
significant increase of scattering was found, corresponding to a width J] of the order 


of 10-° ev. 
No such increase was observed with 18’Ta gamma-rays scattered from tantalum carbide. 
The negative result for 1®’Ta and the positive result for 1°*Hg are consistent with the 
latest information about the life-times of the excited states concerned, viz. 1:1 x 107-8 sec.. 
for 181T'a and less than 2 x 10-!° sec. for 1°8Hg. 


TIN TRODU CianoON 


F a source of mass M emits a photon of energy £, the source will recoil with 

energy £?/2Mc?; an equal kinetic energy of recoil is involved if the photon 

is captured by a body of the same mass as the source. This does not prevent 
the optical excitation of one atom by another, because the widths of optical levels 
are large compared with amount of energy dissipated by recoil; but, owing to 
the high value of FE, it does prevent the emission and capture of a gamma-ray 
from being an effective means of transferring energy of excitation from one 
nucleus to another of identical type. ‘Thus, while the selective scattering of, 
for example, the mercury resonance line 12537 by mercury atoms is of quite 
spectacular prominence, the corresponding nuclear phenomenon has hitherto 
proved unobservable. 

Following Kuhn (1929), various workers have discussed the situation and 
have looked for the resonant scattering. For example, Pollard and Alburger 
(1948) have reported a search for resonant scattering of {Mg gamma-rays 
(£=2-3 Mev.) in magnesium. In this instance the energy dissipated in recoil 
amounts to about 90ev., while the width of the nuclear resonance is certainly 
less than 10-°ev. ‘Though the Doppier effect of thermal motions broadens the 
resonance, and though for heavier elements and less energetic gamma-rays the 
recoil energy can be of the order of | ev. only, the effective energy of the gamma-ray 
is always relatively far out in the low-energy wing of the resonance curve. 


eso.ant Nuclear Scattering of Gamma-Rays ie] 


& 


The present paper reports a theoretical and experimental study of the 
possibility of restoring the resonance with the aid of the Doppler effect, the 
source being made to move towards the scatterer with an appropriate velocity. 


See ETE O Rey. 
Electrodynamical theory* shows that the cross section for resonant scattering is 


Ae je 
° 82 (EB? + 40" ©) sie) eile 16) (1) 


where A(=hc/E) is the wavelength, I’ the width of the resonance, E the effective 
energy of the gamma-ray and £, the resonant energy +; statistical weights are 
here omitted and the transition is assumed to be between the ground state and 
the first excited state of the nucleus. With energies in electron volts, equation (1) 
becomes 
GP ate Vie 
ke (E£-—E,)? +40?" 


oO 


If the nuclei are initially at rest but free to recoil, E, — E = E?/ Mc? as indicated 
in §1; this is a simple consequence of the conservation of energy and momentum. 

Suppose the thermal velocities of the nuclei in the source and the scatterer 
to be distributed as in a gas at temperature JT. ‘The probability that a source 
nucleus and a scattering nucleus shall have thermal velocity components in the 
direction of the gamma-ray that differ by v cm/sec"! is then given by 


P(v) dv =(M/4nkT) 2 exp(—Mv?/4RT)dv. oe. (3) 


The quantity (4RT/M)'? is a measure of the thermal width in cm/sec; it is 
actually the half-width, at 1/e of maximum height, of the distribution curve for 
the relative velocity components. (The intrinsic gamma-ray width I" of 
equation (1) is, in accordance with the usual convention, the whole width of the 
nuclear resonance, measured at half maximum height.) 

The Doppler energy-displacement corresponding to visvE/c. ‘The scattering 
cross section is obtained from equation (1) by replacing E by E+ vE/c, multiplying 
by P(v) and integrating with respect to v. It is 


eee P(v)T? 
a ii: (E—E,+vE,c)? +41? 


ot 


| ae oe (4) 


where the superscript T is a reminder that we are now including the effect of 
thermal velocities. 

If the source as a whole is moving towards the scatterer with velocity u, the 
Doppler effect will add an amount uwE/c to the effective energy E of the gamma- 


radiation. Since E=£,—E?/Mc?, equation (4) now becomes 
ete P(o)l? 

T =o ——_______~_________. dy... =) 
= _o (E4/c*\(ut+uv—E/Me;*+il? z (5) 


* See, for example, W. Heitler, Quantum Theory of Radiation (Oxford : University Press, 1936), 
equation III, 12.17. By a misinterpretation of this equation, Pollard and Alburger obtained a 


value for o that is too great by a factor of 47°. : 
+ Except when the difference between them is involved, EF and Ey can be used interchangeably. 


78 P. B. Moon 


In this integral, the factor [?/[(£?/c?)(u+v—E/Mc)*+4I"] has a very sharp. 
maximum* around v=(E/Mc)—u, and little error will be made by supposing 
the thermal-distribution factor P(v) to remain constant at the value appropriate 
to v=(E/Mc)—u. Writing E/Mc=u,,, we thus find from (3) and (5) 


2 [? dv 
me 1/2 wi — 
oF = 5 (M/4nkT)" exp [— M(t, — with] | (B20) (u— ty, +0) +412 


The quantity beneath the sign of integration is now readily integrable to 27Ic/E, 
whence 


oy 


RD Me? \1/2 
oe Sone CED ie i) te ee (6) 


Using A= he/E, expressing F and I in electron volts, and letting A represent 
the conventional atomic weight and J the isotopic abundance of the resonantly 
scattering nuclear species among the nuclei in the scatterer, one finds 


1/2 
0.x 10-9 (5) exp | - e (ut) | ene Pe (7) 


The maximum scattering cross section occurs when 


w=4,,—H/ Mc=32E/A-cmsce— 


12 
and is Go a0 10° (=) Coe. oe eee (8) 


§3. DESIGN OF EXPERIMENT 


In the experimental arrangement envisaged (Figure 1), a radioactive source 
of gamma-rays moves on a circular path and irradiates (principally when 
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Figure 1. (Not to scale.) 


approaching) a scatterer containing nuclei identical in type with those from which 
the gamma-rays are emitted. A counter, shielded from direct radiation, records 
* The width of this maximum, measured at half height, is eI'/E ; with typical values such as 


=10% ev., T=10-* ev., it is of the order of a few centimetres per second and so much smaller than 
ne thermal width (4kT/M)12 which is of the order of the velocity of thermal agitation. 
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the scattered gamma-rays, and the rate of recording should increase as the 
velocity of the source becomes comparable with the optimum value 32E/A. 

The choice of the isotope to be used as source and scatterer is limited by 
several considerations. ‘The isotope must have substantial natural abundance 
for use as scatterer. It must also be available in the excited state (following the 
decay of its radioactive parent) in order to serve as the source of gamma-rays. 
The gamma-ray to be studied must go to the ground state; if not the only 
gamma-ray emitted by the source, it should preferably be the hardest of them, 
so that the others (together with the much more abundant but softer Compton- 
scattered radiation) can be removed from the scattered beain by absorbers placed 
round the counter. ‘The transition should have as large a width I as possible; 
since the width decreases rapidly as the multipole order increases, dipole or 
electric quadrupole transitions are to be preferred. The optimum speed, 32/4, 
must be within the range of mechanical devices—say, up to 10° cm/sec. Since, 
ceteris paribus, 1 increases with £, it is desirable for the atomic weight A to be 
large so that this last condition can be satisfied with a reasonably high value of E. 

A survey of known radioactive transitions yielded two which seemed 
distinctly the most promising: the 0-411 Mev. transition in !8Au and the 
0-48 Mev. transition in 1*!T'a. The former follows the 2-7 day f-decay of 1%8Au, 
the latter the 46 day decay in '*'Hf, passing through a 22 usec. metastable level 
in '§!'T’a. Both transitions are believed to go to the ground state, each is the 
most energetic radiation of appreciable intensity from the nucleus in question, 
and both are believed to be electric quadrupole. ‘The isotopic abundance of 
ePTigis 10%; of *'Ta,400%. The optimum speeds are 6:3 x 10* cm/sec. and 
8:5 x 104 cm/sec. respectively. ‘The thermal breadth of the response curve is in 
each case such that the intensity of resonant scattering should fall to half its 
maximum at 80°% of the maximum speed and should be negligible below half-speed. 

Selecting 18Au for detailed discussion, we find from equation (8) that the 
maximum cross section for resonant scattering is 4-6 x 10-71 [ cm?; it is here 
assumed that 7 =360°K., since the rotor which carries the source runs at about 
420°K., while the scatterer is at room temperature. It is clearly desirable to 
use large angles of scattering, so that the Compton photons may be the more 
easily filtered out. Let us therefore consider scattering through an angle of 115° 
(the average angle used in the experiments). ‘The angular distribution of the 
resonant scattering will depend upon the multipole order of the transition but 
is unlikely to be violently non-uniform, so the resonant scattering cross section 
per unit solid angle may be taken as equal to o/47 or 3-6 x 10-*% [.cm?. At 115°, 
Compton scattering and Rayleigh scattering have cross sections per unit solid 
angle of about 10-24 and 10-26 cm? respectively. If Compton-scattered radiation 
can be adequately filtered out, an intrinsic width P of 10-*ev. should give readily 
observable resonant scattering, since the Rayleigh scattering has been observed 
with no great difficulty (Moon 1950, Storruste 1950). 


§4. EXPERIMENTS WITH ‘Hg 


The tips of a doubly tapered steel rod were electroplated with gold, and the 
whole was irradiated for several days in the Harwell pile (BEPO). A few days 
after irradiation, the activity was of the order of 100 mc. and was mainly from 
the gold plating. ‘The rod was then spun in vacuum about an axis perpendicular 
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to its length, the speed of the tips being taken up to the limit of safety of about 
7 x 104 cm.sec-! and down again; the top speed of the centre of mass of the gold 
was 6x 104. Meanwhile, observations were made of the rate of counting of a 
Geiger-Miiller counter shielded from direct radiation but exposed (through an 
1-inch lead absorber) to gamma-rays scattered from a surrounding thin-walled 
iron-alloy cone containing liquid mercury (Figure 1). ‘The cone was placed so 
as to be exposed mainly to gamma-rays from the advancing tip of the rotor. 

Four complete experiments were made, each lasting for about 1} hours and 
each involving the registration of upwards of 250,000 gamma-rays; corrections 
(unimportant to the final result since acceleration and deceleration occupied 
about the same time) were made for the experimentally observed decay of the 
source (about 0-7°% per hour). ‘The first two runs were made as described above. 
In the third, the direction of rotation was reversed; a smaller effect would be 
expected owing to the less favourable position of the rotor tip when advancing 
towards the scatterer. ‘The fourth run was made in the forward direction with 
a scatterer of copper instead of mercury; any increase at high speed would in 
this case be due to extra-nuclear phenomena such as stretching of the rotor. 
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Figure 2. Effect of speed of source on rate of recording of scattered gamma-rays. Negative 
speeds indicate rotation in opposite sense to that shown in Figure 1. 


During a fifth run, with a double thickness of lead round the counter, a vacuum 
failure before full speed had been reached caused the rotor to strike the wall of 
the vacuum chamber, with catastrophic results to both. 

For purposes of illustration, the results for the second ‘ forward’ run and the 
‘reverse’ run, which were made on the same day, are plotted together in Figure 2. 
Each point represents the number of particles recorded in a ten-minute interval, 
and the mean speed during that interval; circles refer to readings taken during 
acceleration, crosses to readings taken during deceleration, while the heavy cross 
represents a reading taken at very low speed between the two runs. The vertical 


lines show the probable error, calculated from the number of particles observed: 


in each interval. ‘The broken lines show a possible analysis into background 
and resonant scattering, varying with speed in the expected manner and more 
intense (as it should be) with ‘forward’ than with ‘reverse’ rotation. Such an 
analysis might be over-ambitious and it is preferable to rely on the ratio of the 
mean counting rate at all speeds above 4 x 104 cm. sec~! to the mean rate at all 


} 
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lower speeds. The two ‘ forward’ runs gave values for this ratio of 1-007, + 0-008, 
and 1-015; +0-006, the probable errors being calculated from the experimental 
fluctuations of counting rate within each of the two speed ranges; since any 
genuine increase will vary with speed, the errors may be overestimated. The 
“reverse’ run gave a ratio of 1-005, + 0-005, and the ‘blank’ run, with a scattering 
cone of copper instead of mercury, gave a ratio of 0-998, + 0-005. 

The difference between the mean of the two ‘forward’ ratios and that for 
the blank experiment is 0-013+0-007. This result is distinctly suggestive of 
the presence of resonant scattering, and it seems worth while to deduce the 
corresponding values of I’ and of the half-life of the excited state. The figure 
of 0-013 represents, crudely, the ratio of counts due to resonant scattering to those 
from Compton scattering, both at a mean angle of 115°, but it must be corrected 
on account of their different chances of emergence from the thick scatterer, 
their different transmissions through the absorber surrounding the counter, 
and the different sensitivities of the counter itself to the two energies in question 
as well as for background of various origins. It has also to be remembered that 
only those gamma-rays that leave the source nearly in its direction of motion 
will receive the full Doppler hardening, and that the experimental ratio is an 
average over speeds ranging from 4 x 104 cm/sec. to 6 x 104 cm/sec. With these 
factors taken into account, I is found to be about 3 x 10-°ev., corresponding to 
a half-life of the order of 10-1! sec. for the 0-411 Mev. excited state of *Hg. 

Shortly after these experiments were completed (April 1949), this half-life 
was reported to be about 2x 10-® sec. on the basis of delayed-coincidence 
measurements (MacIntyre 1949). If this were so, resonant scattering would 
be about two thousand times less than the present work indicated. Because of 
this contradiction, plans were made to verify the scattering with a different 
experimental arrangement. ‘This has now been done with the help of 
Mr. A. Storruste and Mr. T. H. Bull; the effect has been qualitatively confirmed 
but the detailed analysis of the results, involving various auxiliary measurements, 
will take some time to complete. In the meantime, the contradiction has been 
removed by the work of Bell and Graham (1950), who find the life-time of the 
excited state to be shorter than the limit of resolution of their apparatus, which 
is12X10-" sec. 

§5. EXPERIMENT WI?PH Ta 

A similar experiment was made with '§'T'a as the emitting and scattering 
isotope. The source was about 8 mg. of Hf,Os, irradiated in the Harwell pile 
for two months to obtain about } mc. of the 46-day ''Hf. ‘This source was 
contained in small cup-like odes in the ends of a rotor which could withstand 
higher speeds, and the apparatus built for this experiment differed in other 
details from that used earlier for 18Hg. ‘The scatterer was tantalum carbide. 
In spite of the decreased statistical accuracy due to the much weaker source, the 
experiment should afford a more sensitive test of resonant scattering than was 
possible for !°8Hg, since 18!Ta has 100% abundance. ‘T'wo runs, in which the 
counting rates from 4 x 10*to 9 x 10* and from 0) to 4 x 104 cm/sec. were compared, 
gave ratios of 1-003 + 0-015 and 0-990 + 0-014, with a mean result of 0:9965 + 0-01. 
It is to be concluded that the 0-48 Mev. y-transition either does not go to the 
ground state or has a width less than 10-*ev. and hence a life-time greater than 
about 4x 10-1 sec. After this measurement had been made, a y-transition of 
life-time 1-1 x 10-8 sec. was reported (Barber 1950) which may plausibly be 
identified with the 0-48 Mev. transition in question. 
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Detection of u»-Mesons and other Fast Charged Particles in 
Cosmic Radiation, by the Cerenkov Effect in Distilled Water 
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ABSTRACT. The Cerenkov effect has been established for the passage of cosmic-ray 
p-mesons through distilled water. The light pulses were observed using a photomultiplier 
and fast amplifier; the experiments were carried out using the detector alone, and in a 
coincidence-telescope arrangement. The ratio of the light emitted in the downward and 
upward directions was measured as a function of discriminator bias. In addition measure- 
ments were carried out at various depths of water and with various orientations of the 
detector with respect to the vertical. Estimates of the absolute efficiency were made and 
other characteristics of the detector are discussed. 
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URING the course of investigations into the possibility of utilizing various. 
1) organic solutions as scintillation detectors in cosmic-ray applications, it 
was found that several liquids of commercial purity produced light pulses 
from the passage of single cosmic-ray particles. A study was then made, using 
distilled water, to ascertain whether some of the pulses could be attributed to the 
Cerenkov effect (Cerenkov 1934). ‘This effect, the theory of which was worked out 
by Frank and Tamm (1937), has been investigated by several workers in recent 
years, notably by Collins and Reiling (1938) who employed a photographic 
technique, and by Dicke (1947) who used a photomultiplier and an optical system 
proposed by Getting (1947). The conclusive evidence for the effect has been 
limited to detection of artifically accelerated electrons and secondary electrons. 
produced by y-rays. Among those who searched for the phenomenon with cosmic 
rays Weisz and Anderson (1947) found a small effect while Dicke obtained no effect 
that could be attributed definitely to Cerenkov radiation, 
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§2. EXPERIMENTAL TECHNIQUE 


The arrangement adopted in the experiments to be described is shown in 
Figure 1. A tray of Geiger counters of effective area 120cm? was mounted 
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Figure 1. Essential features of apparatus. 


vertically above a glass cylindrical container completely filled with distilled water. 
The container, of length 20 cm. and diameter 8-5 cm., was made of glass, as were 
also the end plates; the curved body of the container was silvered to reduce loss of 
light. Black paper over the upper end of the container was arranged to absorb 
light emerging upwards while, at the lower end, a small light-integrating cone 
coated with magnesium oxide was mounted to fit snugly round the cathode K of 
an 11-stage E.M.I. photomultiplier (type 5311). The water container, white 
cone and phototube were all mounted in an outer light-tight brass envelope E. 

The output pulses from the phototube, after passing through a low-gain 
pre-amplifier situated at the base of the tube, were amplified by a main amplifier, 
after which they passed through a discriminator to one channel of a coincidence 
mixer of resolving time r=1-8ysec. ‘The other channel was connected to the 
output of the counter tray via a cathode follower. ‘The amplifier was of special 
design having both integration and differentiation time-constants of 3:2 x 10-8 sec. ; 
this extremely narrow bandwidth resulted in a considerable reduction of the 
phototube noise relative to a pulse as short as that computed for the Cerenkov 
effect, in this case of the order of 10-°sec. ‘The phototube was operated at 1-9 kv. 
and the overall gain of the amplifier was about 100. 

Water was chosen for the light-producing medium for it was believed there 
would be a complete absence of the ‘scintillation type’ mechanism of light 
production. The cylindrical geometry was chosen as this was symmetrical about 
the centre S of the container when this and the phototube were inverted. 

The experiments were divided into two groups: (i) coincidence experiments, 


carried out under conditions of relatively low bias at a time when it was not known 
6-2 
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how large the light pulses would appear from this type of radiation ; (11) experiments 
conducted at various conditions of bias using the Cerenkov detector alone and the 
pulses counted directly at the output of the discriminator. 


§3. RESULTS 


(i) Coincidences between the Geiger- Miiller Tray and Cerenkov Detector 

A number of runs was carried out in which the coincidence rate was measured 
between the counts in the Geiger tray and those in the phototube ; the chance rates 
were calculated in the usual way and subtracted from the observed rates, for this 
purpose the single rates were monitored on scaling units and pen recorders. It was 
realized that a crucial test for Cerenkov radiation lies in its characteristic of emission 
in directions making angles of less than 41°* with that of the incoming particle (this 
figure refers to the case of water). Coincidence runs were carried out under two 
conditions, firstly with the phototube in the position shown in Figure 1 with its 
cathode at K below the container (this position is denoted by B), and secondly, 
with the phototube cathode at K’, above the container (denoted by position A). 
Position A was brought about by rotating the whole detector assembly through 180° 
in a vertical plane about the centre-point S of the container; the counter tray was 
not moved. Runs in the two positions were alternated and interlaced with runs 
in which the container was emptied, in order to allow for those coincidences 
arising from the passage of particles through the phototube alone in its two 
positions K and K’. In Table 1 are presented the averaged results of these runs, 
carried out at three bias settings on the discriminator, 2, 5 and 15 volts; the figures 
refer to genuine coincidence rates (i.e. chance subtracted) in counts/minute, and 
the symbols P, C and L refer to the phototube, container and liquid respectively. 


ables). 


Bias 2 volts 5 volts 15 volts 
Phototube B A B A B A 
position 
(P+C+L) 4:14+0:03 2:29+0-10 4:01+0-11 1:73+0-04 2:54+0-08 0-71+0-03 
(P+C) 0:79+0-03 1:26+0-02 0:47+0:01 0-98+0-03 0:15+0:01 0:-40-+-0-02 
L  3-35+0:04 1:03+0:10 3:54+0-11 0:75+0-05 2:39+0-:08 0-:31-+0-04 
Ratio [Lp/La] 3-25+0°32 4:72+0:°35 7-741:0 


The average value of the chance rate, only significant at the 2-volt bias setting, 
was approximately 0-2/minute; this was deduced from the counter-tray and. 
phototube rates which fluctuated around figures of 400 and 7,000 counts/minute 
respectively. From the above results the ratio between the counting rates in the B 
and A positions, Lz/L,, see Table 1, rises from 3-3 at 2 volts bias to 7-7 at 15 volts. 
These ratios are good evidence in favour of the Cerenkov effect taking place in 
water due to the passage of sea-level cosmic-ray particles. 

A control experiment was then carried out using a 0-5°, solution (by weight) of 
terphenyl in xylene (Reynolds et al. 1950) to replace the water as a source of 
fluorescence. Ifitis assumed that the fluorescence radiation is emitted isotropically 
from each point in the track of a particle one would expect a value for Ly/L, for the 

* In general, the angle of light emission @ is related to the refractive index m, and the ratio of the 


velocity of the particle to that of light in vacuo B, by cos 0=1/Bn. For water (n=1-33) and B > 1-0 
41° 2 
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solution of less than unity because of the larger flux of tracks near to the photo- 
cathode in position A (Figure 1). The ratio obtained for Ly/L, in this experiment 
at a bias of 2 volts was 1:26+0-06. It would have been expected that this ratio: 
would have been unity or less if the Cerenkov effect in this case were negligible ; 
however, the observed ratio is explicable if a comparable contribution is present 
from the Cerenkov effect in the xylene. 

That the ratio Ly/L, in water is not infinite is interpreted as being due to 
internal reflections within the container such that even in position A sufficient light 
is received in a fraction of cases to give acount. The increase of this ratio with 
discriminator bias confirms this hypothesis. Taking the accepted value for the 
total cosmic-ray intensity at sea level (Rossi 1948), the coincidence rate of a 
hypothetical detector 100% efficient and of area 57 cm? placed at S would be about 
7/minute. ‘This figure is to be compared with that of about 3-4/minute obtained 
experimentally for the phototube in position B at 2 volts bias. In the B position 
with the terpheny]l solution a liquid-alone rate L of about 3-5/minute was obtained 
which agrees with that for the water within the errors. 

It is interesting to note that at 2 volts bias a coincidence rate of 0-55 + 0-03 per 
minute was obtained for the phototube alone in position B which agrees very 
closely with that calculated, 0-6/minute, from the total vertical intensity and solid 
angle, when an area comparable with that of the photocathode is assumed to have 
an efficiency of 100°%. Which elements of the E.M.I. phototube are in fact 
responsible for this was not investigated. It could well be interpreted as due to the 
emergence of ‘ delta rays’ from the cathode surface during the passage of the 
particles through the glass envelope. : 

A test was carried out with a lead absorber 10 cm. thick placed above the counter 
tray to ascertain, by absorption of the soft component, that the observed effects. 
were due to p-mesons as well as to electrons. With this arrangement the 
coincidence rate L for the water alone in position B was measured to be 
2-67 +0-11/minute as compared with the value of 3-35 +0-04/minute with no. 
lead. It can therefore be concluded that approximately 80% of the effects 
obtained were due to p-mesons. 


(ii) Single Rate Cerenkov Detection 


It was found that above a bias of about 10 volts the single counting rate from the 
water detector alone depended markedly upon its orientation and that the phototube 
‘noise ’ then made only a small contribution to the total counting rate. A number 
of bias curves was then taken, see Figure 2, in which the counting rates of the 
Cerenkov detector alone were measured at values of bias up to 100 volts; the gain 
of the system and amplifier time-constants were maintained the same as in the 
coincidence experiments. ay 

Curve (a) was obtained for the phototube and empty container in both positions 
A and B; the counting rates obtained at each bias setting were the same, within the 
error limits, for the two positions. Curve (4) was obtained in position B when the 
container was filled with distilled water while curve (c)was obtained with the detector 
inverted, position A. Subtraction of curve (a) from curves (6) and (c) yielded bias. 
curves (d) and (e) for the water alone in positions B and A respectively. In the 
inset to Figure 2 is shown the ratio of the counting rates Ly/L ,obtained from these 
last two curves; it will be seen this ratio is as high as 70:1 at a bias of 30 volts. 
Curve (f) was obtained at half depth (i.e. 10 cm. of water) ; as expected this lay at a 
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lower level of counting rates. At full depth, curve (d), the extrapolated value of the 
counting rate at zero bias was approximately 170/minute. This is considerably 
larger than the rate calculated, of the order of 70/minute, for the total cosmic-ray 
radiation entering an area equal to the cross section of the cylinder integrated over 
the hemisphere ; however, it must be remembered that at all angles of incidence 
other than the vertical the projected area of the container is larger than 57 cm?; it 
rises in fact to 170 cm? for particles entering horizontally. 
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Figure 2. Single-rate bias curves; inset, the ratio LR/La. 
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The variation of counting rate with depth of water was carried out, at a fixed 
bias of 40 volts; a plot of these results is shown in Figure 3. 
The curve indicates, at this bias, a detection efficiency of less than unity since an 
asymptotic value of the counting rate even at 20 cm. depth is not attained. 
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Finally, the variation of counting rate of the detector was determined as a 
function of the angle 6 between the axis of the system and the vertical. These 
results, carried out at a bias of 30 volts, are shown in Figure 4. 

The general form of the curve agrees with that expected when it is remembered 
that the light radiated in the water by a relativistic particle is emitted at an angle of 
approximately 40° with respect to the track of the particle. This means that there 
is a considerable counting rate when the detector is lying in a horizontal plane, and, 
for angles of # such that 180°>0>90°, detection at reduced efficiency still occurs. 
For this and other reasons connected with the geometrical form of the detector and 
the bias conditions, it would not be expected that the form of the curve in Figure 4 
should follow the cos*@ distribution law, though for comparison purposes this has 
been included in the figure. 


§4. CONCLUSIONS 


The Cerenkov effect has been firmly established for the passage of cosmic-ray 
p-mesons in distilled water. It now opens up the possibility of developments of 
methods of measuring directly the velocities of singly charged particles with values 
of 8 lying between 0-75 and 1-0 (for water), as first suggested by Getting (1947). 
It is not easy to understand why the Cerenkov effect with cosmic-ray particles had 
not previously been obtained using the phototube technique. It may be that in 
the present experiments the phototube noise was less than in the tube used by 
Dicke in his earlier experiments, or that the extremely fast amplifier available to the 
author enabled a very much better signal/noise ratio to be obtained. 

Further characteristics of the new detector for cosmic-ray research are 
(i) extremely high speed of operation, (ii) high efficiency, (ii) negligible 
sensitivity to nucleon stars and (iv) a very small sensitivity for y-rays from radio- 
active sources such as the laboratory background. While properties (i) and (ii) 
hold for the best scintillation counters, (111) and (1v) are unique to this instrument. 
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LETTERS TO THE EDITOR 


The Angular Distribution of Protons Emitted in some 
F(«, p) Ne Resonance Reactions | 


The angular distributions of the protems emitted im some “Fi, py Ne resenaree 
reactions have been studied by using a Giumtegration camera and the photogeaphic emuliem 
technique to record the emitted particles. The destgm ef the camera has alread beem — 
desertbed by Jolley (199). A collimated beam of a-parteles from: 2 pelenmum source ef © 
about 3 milheuries streneth is ebtaimed by means of three separated coamal discs of cuppes. 
The energy of this beam on striking the target is adjusted te that of 2 kmown resemance level 
in the compound nucleus “Na by varying the distamee between source and target: ar at 
atmospheric pressure 8 used as absorber. The target of Guemme censists ef @ laver of 7 
evaporated calcrum fluoride om very thim gold foil, the Guoride laver having am esimard — 
thickness of 3 mm. air equivalent. 

With the energies of the a-particles set at already established levels of $1 and ST we 
about 2,600 and 900 protons were recorded respectively. The mtensty-eserey dsmmbunens 
give m each case two distinct peaks together with indications ef amether smaller peak> 


a Om ten ie ff Bocce aa mma * 


a tt 


hh 


these peaks correspond to transitions from: the level im the campeund aucieus te diffrent 
levels in the fimal nucleus. The energies of the peaks for the $-1 wev_ level are }-6—#28 = 


Number per Unit Solid Angle 


Nurher per Unit Solid Angle 


a aS ET 8S (ga 
Angiz(degrees}_ centre-af-mess coordinates 


Angie cegreeS. catrea es ordinetas 
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Angular distribution of protoms from the F(z, p)Ne reacnem (@) Ground seme cansiens- 
(@) Second exerted state mansitiers. 
Pigure 1, E,>=*-t ev; Figure 2. E.— tT wex 


and 3-67—0-13 wev_ with indications of another peak at about $4 Mev_- fae the $e 
level the peaks are at 4-64+0-10 wev_ and 3-22+@-11 wev_ with mctcsnens of amotter ; 

at 4-2utev. The spacings of these energy levels im the fimal nucieus “Ne are & G+ ang 
1-41 wev_ zhove the ground state. 

The angular distributions of the protons from: these reactions may be seem im Fizures ¥ 
and 2. The ground state transitions both show distibutions having 2 minimun: at W* amd 
@ mazimum at 0 to the direction of the abeam. The second eacited state tansitiens batik 
show mmimz2 at 0° and 90°- there is 2 maximoume at 30° mm the case of the £-1 ew resomane= 
and 2 maxmmure at 45° for the 3-7 mev_ resomance Nerher of the Gist excited sare 
transitions is sufficiently intense to allow 2 determination of the angular Gamcbutens. 
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The recorded distributions are in all cases consistent with symmetry about the equatorial 
plane but because of the geometrical limitations of the apparatus measurements can be made 
only between 0° and 135°. There is a striking similarity between the two corresponding 
angular distributions for the ground state transitions for both resonance levels. This 
similarity persists in the second excited state distributions and it may be that the two levels 
of the compound nucleus overlap so that they behave effectively as one level as far as the 
angular distribution is concerned. 

The relative yields for transitions to the ground, first and second excited states are in the 
ratios 1 : 0°35 : 4 and 1: 0-16 : 3-1 for the 4-1 and 3-7 Mev. resonances respectively. Early 
experiments of Chadwick and Constable (1932) gave the yield ratio as 1:4 for ground and 
first excited state transitions. Our value of 1 : 0-35 is in better agreement with results of 
May and Vaidyanathan (1936) who were unable to record any transitions to the first excited 
state presumably because of its low intensity. 

One of us (J.D.J.) wishes to thank the Department of Scientific and Industrial Research 
for a maintenance grant which made it possible to carry out this work. 


Wheatstone Laboratory, J. D. Jo.tey. 
King’s College, F. C. CHAMPION. 


University of London. 


Submitted as a paper 27th March 1950; 
resubmitted as a letter 15th September 1950. 


CHADWICK, J., and ConsTaBLE, J. E. R., 1932, Proc. Roy. Soc. A, 135, 48. 
Joey, J. D., 1949, Ph.D. Thesis, University of London. 
May, A. N., and VAIDYANATHAN, R., 1936, Proc. Roy. Soc. A, 155, 519. 


The Thermal Conductivity of Solid Helium 


We have measured the thermal conductivity of solid helium at temperatures between 
1-:2° and 2:5°x. A german silver tube of 6 mm. 1.D., 8 mm. 0.D. and 5 cm. long was filled 
through a capillary with liquid helium at a temperature and pressure corresponding to a 


0:218 gm/cm’ 
0-214 
0-208 
0-203 
0-194 


mMoOO> 


Watt/cm deg 


o 


0s 


20 
Temperature (°K.) 


Thermal conductivity of solid helium at constant density. 


point on the melting curve, where the density is known. The capillary was then closed off 
and the tube and its contents cooled down. Provided that the liquid had been sufficiently 
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compressed initially, the variation of the densities of solid and liquid helium is such that the 
tube will always be completely filled by the solid. Thus all the measurements have been 
made at constant densities whose magnitudes were fixed by the initial pressure and 
temperature of the liquid. One end of the tube was mounted on a cryostat the temperature 
of which could be varied, and an electric heater was used to supply heat to the other end. 
The conductivity was then found by measuring the temperature at each end with leaded 
phosphor-bronze thermometers. The corrections for the conduction of the metal tube, 
for heat flow down the leads and for end effects were only a few per cent. One set of 
measurements at a density of 0:208 was also made using a larger tube (9-9 mm. 1.D. and 
10 cm. long); the points so obtained are shown in the Figure and lie on the same curve as 
that obtained with the smaller tube. 

The Figure shows the thermal conductivity as a function of temperature for various 
densities; different symbols on the same curve correspond to experiments on different 
occasions. The only previous value known to us is that of Kikoin (1939), who gives a value 
of 3 x 10-4 watt units at a temperature of 1-8°x. but does not state the density; this agrees in 
magnitude with our results. 

Although for practical purpose it is convenient to use conductivities, it is better for 
the theoretical approach to consider their reciprocals, the thermal resistivities. "The total 
resistivity of a cylindrical crystal may be given approximately by 


W=AT%e—§2T 4. 6T+-1/CDT?. 


The first two terms (Peierls 1929) give the resistance of an infinite crystal due to 
interactions between the elastic waves, brought about by ‘umklapp processes’ and 
imperfections of the lattice, in terms of three constants A, v, b, and the Debye specific heat 
parameter 0. The last term (Casimir 1938) is due to the scattering of the waves by the 
boundaries of the crystal and is only important at low temperatures. C is a constant and D 
the diameter of the crystal. 

It appears that in the region of our measurements the principal contribution to the 
resistance is from the first term. No values of the constants A and wv are available save 
that wv is of the order unity, so we write to a first approximation W «e927, which gives a 
thermal conductivity K « e%/2T, As the value of @ for solid helium does not vary too much 
with temperature, we treated it as a constant and chose values of it so as to get the best fit 
with our experimental points. The appropriate values to fit the curves at densities of 
0:218 and 0-194 are 26° and 21° respectively; these are similar to those derived from 
specific heat measurements (Keesom and Keesom 1936) and vary with the density in the 
same way. 

At really low temperatures the only term of consequence in the expression for the 
thermal resistivity will be the first one, so that eventually K=CDT*. Thus a maximum 
must occur in the curves at some temperature lower than we have measured. To make an 
estimate of this temperature we put CDT*=ae%2?. Substituting in values for all the 
constants, we find that for the same sized specimen the maximum should be between 
0:5° and 1:0°.. Klemens (to be published shortly) has shown that the conductivity in the 
region of the maximum is reduced by a factor of about ten, due to interactions, so we would 
expect the conductivity there to have a value of the order of 2 watt units. We hope shortly 
to extend our measurements below 1°. to verify these predictions. 

Solid helium is very compressible and is, therefore, a good substance on which to study 
those pressure-dependent constants which cannot yet be evaluated theoretically. The 
variations in the densities of our specimen were obtained by varying the initial pressure 
over the liquid through the comparatively small range of 60 to 136 atmospheres, yet this 
was sufficient to vary the thermal conductivity by a factor ten. We intend to make measure- 
ments using much higher filling pressures, of the order of a few thousand atmospheres. 

In view of the high conductivity of solid helium it is interesting to review the suggestion 
made by Dr. C. A. Swenson that a capillary filled with helium might be used as a heat 
switch below 1°: liquid helium II in the tube would give a very big conductance, while 
on applying pressure it would solidify and have a much lower conductivity. In the 
course of measurements on the specific heat of solid helium below 1-0° we suspended a 
pill of iron alum at 0-3° by a fine capillary (0-13 mm. 1.D.) from a cryostat at 0:9°. With 
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solid helium in the capillary, the heat flow corresponded to a mean conductivity of about 
107! watt units, which is consistent with our other measurements (after allowing for size 


effect). On releasing the pressure, the liquid equalized the temperatures almost 


instantaneously, which is consistent with measurements of the heat flow in a capillary just 
above 1-0° (e.g. Allen and Ganz 1939). Thus a heat switch is possible in this region, but it 
must be emphasized that the conductivity of the solid helium is still considerable, and 
precautions need to be taken to prevent too great a heat flow down it. At much lower 
temperatures, say below 0-3°, it seems likely that the conductivities of both liquid and solid 
will be small and similar, both being limited by size effect. 

We both wish to express our appreciation to Professor F. E. Simon for his kind 
interest in the work. One of us (K.R.W.) is indebted to the Department of Scientific 


and Industrial Research for a maintenance grant. 


Clarendon Laboratory, K. R. WitkINson. 
Oxford. J. WILks. 
13th November 1950. 
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Factors Involved in the Accuracy and Reproducibility of Depth 
Measurements on Nuclear Research Emulsions 


During the early stages of our work with boron-loaded Nuclear Research emulsions, we 
observed considerable variations in the apparent thicknesses of the processed emulsions. 
‘The plates used were Ilford Type C2, of nominal unprocessed thickness 100 ~. Thicknesses 
of the 22 processed plates examined ranged from 33-504; over any one plate there were 
variations of as much as 10 pp (25°), and even measurements on the same area of a plate made 
at different times varied by as much as 44 (10%). Since these variations would apply to all 
depth measurements, the following factors, which we considered may have caused them, 
were investigated : (i) the reliability of the measuring instrument; (ii) the variation of the 
emulsion thickness with relative humidity; (ili) the effect of oil immersion on the emulsion. 

(i) The reliability of the measuring instrument. "Tine thickness of the emulsion was 
measured by focusing a microscope first on the uppermost grains in the emulsion and then 
on the lowermost grain in the same area, and reading the distance traversed on the depth- 
scale of the microscope. The reliability of two microscopes was tested by calibrating the 
depth-scales using an interference method. One of the microscopes was found to have a 
periodic error which gave rise to considerable variations in the depth measurements. The 
other microscope gave a consistently uniform calibration, i.e. 1:02; u per division +1%, and 
this instrument was used for the rest of this work on depth measurements. The uncertainty 
of the depth measurements when using this microscope was estimated to be +1 p, and 
variations outside this could be regarded as true variations. 

(ii) The variation of emulsion thickness with relative humidity. As considerable variations 
in thickness of emulsions still seemed to be present, we next examined the effect of relative 
humidity on emulsion thickness. "Three plates, two from the same batch A and one from a 
different batch B, were placed in turn in constant humidity baths ranging from 0-100% 
relative humidity and measurements were taken on four areas of each plate. All plates 
examined showed similar increases in emulsion thickness with humidity, the thickness more 
than doubling itself between 0 and 100%, relative humidity*; typical results are shown in 

* S. L. Martin (Aust. ¥. Sci. Res., 1949, 2, 389) refers to this change of thickness with humidity. 


He found that the shrinkage factor k=(thickness after processing)/(thickness before processing) 
changed from about 0:5 at 0% r.u. to about 1°3 at 100% r.u.; this is in agreement with cur results. 
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curves 1 and 2 of the Figure. In order to test the reproducibility of results, measurements 
were made on one area of each of six plates at the following relative humidities : 0, 43-5, 53, 
64, 76 and 88, and then on the same areas, and at the same humidities in the reverse order. 
Typical results given in curves 3 and 4 show that the measurements repeat themselves 
within the estimated experimental spread of 2 p. 
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Variation of emulsion thickness with relative humidity. 


Since the range of humidities encountered in the Laboratory over a period of five months 
was found to be 35-80% and the variation in thickness between these two humidities was of 
the order of 25°% for the plates measured, it is necessary to work at constant humidity. 
As a standard procedure we now place all plates in a constant humidity bath of 53% R.H. 
(this being approximately the mean R.H. in the Laboratory) for at least 24 hours before making 
any depth measurements. On removing the plate from the constant humidity bath, the 
area on which measurements are to be made is immediately covered with cedar wood oil; 
it has been found that this effectively seals the emulsion and prevents any changes due to: 
varying humidity in the room for normal periods of use, i.e. 5-6 hours. Check measure- 
ments at 53°4 R.H. after a period of two months confirmed that, provided the standard 


procedure outlined above is used, measurements repeat themselves within the estimated 
spread of +1 p. 


Since the emulsion normally adheres to a glass surface, it was anticipated that there 
would be no appreciable lateral expansion due to change of humidity. This was confirmed 
by observations at 0 and 88° R.H. on thorium stars produced in the same type of plates. 

(iii) Effect of oil immersion on the emulsion. During the course of the work it was 
suspected that when plates were left covered in oil for long periods, swelling was liable to 
take place due to the absorption of oil by the emulsion. To investigate this, two plates, 
one from each of the batches previously investigated, were placed in a constant humidity 
bath of 53%, R.H., one test area on each being covered in oil all the time and another area 
being used as a control. Measurements were taken over a period of nine days. It was 
found that in one plate, from batch B, both the test area and the control area showed no 
change in thickness, outside the experimental error of -|+1 4, over the whole period; whereas 
in the plate from batch A the thickness of the test area showed a tendency to increase slightly 
with time (the increase exceeding experimental limits after about 7 hours), while that of the 
control area remained constant. This result is rather inconclusive, but it does appear that 
plates are liable to absorb the oil if left for long periods, and it is considered advisable to: 
remove the oil from the plates as soon as measurements are finished. 

The work described above has been carried out as part of the research programme of the 
National Physical Laboratory, and this letter is published by permission of the Director 
of the Laboratory. ‘The authors desire to acknowledge the comments and suggestions they’ 
received from their colleagues in the Physics Division, National Physical Laboratory. 


Physics Division, J. M. M°ALIsTER.. 
National Physical Laboratory, D. W. Kream. 


Teddington, Middlesex. 
11th October 1950. 
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Measurement of Gamma-Ray Momenta and the Thresholds of 
the Photo-disintegrations '*C—- 3*He and 'O-> 4*He 


In a previous letter (Goward and Wilkins 1950 a) the accuracy obtained in measurements 
of various charged-particle photo-disintegration stars was described. 'The stars were 
produced in nuclear emulsions by irradiation with continuous gamma-ray spectra of energies 
up to 25 Mey. For each star, assuming its identity, the observed momentum unbalance 
vector A was calculated by the relation 


A (2D,)—Pyr e 2 OO EE eee (1) 


where p, (r=1, 2, ....) are the individual particle momenta, in (Mev. X mass number)}/? 
units, and p,, is the incident gamma-quantum momentum. The value of p, was obtained 
from ; 


Py USAT eh), eee ee ee eee (2) 


where Ep (Mev.) is the total energy of the charged particles and Fp is the reaction threshold 
calculated from the mass values. The accuracy of measurement, and the certainty with 
which the stars could be identified, was discussed in terms of | A] statistics. 

Momentum unbalance analysis of stars can be adapted to measurement of reaction 
thresholds Ep. From equation (1) it follows, since A should be zero, that &p,, resolved along 
the known gamma-ray direction, gives an independent value for p,, and hence Ep may be 
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Measurement of thresholds for photo-disintegrations. (a) *C—+34He; (6b) “O—44He. 


deduced from equation (2). The errors in the individual Ep values are large and a consider- 
able number of stars must be measured to give a reasonably accurate mean value. The 
method cannot be compared with those using monochromatic gamma-ray sources of known 
energy, but it is the only recourse when suitable monochromatic sources are not available 
(e.g. no available monochromatic source gives O->4*He stars); moreover it provides 
additional evidence on the correct identification of the stars. 

The Figure (a) shows a histogram of Ep values obtained, in the above manner, from 
275 2C + 34He stars observed in emulsions irradiated with the gamma-ray direction 
parallel to the emulsion surface. Momentum unbalance statistics for the stars indicated 
that each p, value would be subject to a probable error of about +0-19 and consequently 
each Ep value to a probable error of +8-2 Mev., since the errors in Ep are relatively small; 
this is confirmed by the fitted curve. The mean value of Ep for the stars of the Figure (a) is 
(7:3-£0-5) Mev., in agreement with the expected value, 7-3 Mev. 

The Figure (6) shows a histogram of Ep values for 111 “O— 4*He stars. Here the 
individual p, values had a probable error +-0-23 (see Goward and Wilkins 1950 b) the 
increase arising mainly from the additional track to be measured in each star. ‘The fitted 


curve thus corresponds to a probable error of —+9°9 Mev., and the mean value of E’p obtained 


from the histogram is (14:8-L1-0) Mev., consistent with the expected value, 14:5 Mev. 
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It is interesting to observe that the above techniques, applied in a reverse way, provide 
a test of the equation p=hv/c or, in the present units, p,=0:0232E,, where E,, is the gamma- 
quantum energy. Assuming the reaction threshold Ep, Ey, is derived from E,=Ey+ Ep. 
and p, is derived from the momentum balance analysis. For both the *C + 3*He and the 
16Q -» 44He stars discussed above, the probable error of the individual (p,/E,) values is 
about +0:01. The combined statistics give py/E, =0-0232 +0-0005. 


Atomic Energy Research Establishment, F. K. Gowarb. 
Harwell, Didcot, Berks. J. J. WIdLkINs. 


18th September 1950. 


Gowarp, F. K., and Witxrns, J. J., 1950 a, Proc, Phys, Soc. A, 63, 662; 1950 b, Ibid., 63571716 


The Photo-disintegration of Oxygen into Two “Be Nuclei 


Continuing the work reported in a recent note (Goward and Wilkins 1950), measure- 
ments have now been made on some 200 1®O - 4 *He photo-disintegration stars. As already 
reported, analysis of the first 100 stars showed the predominant modes of disintegration 
to be via intermediate !2C nuclei, no definite evidence being obtained for direct disintegration 
into two §Be nuclei. That such direct disintegrations do occur, however, has been established 
by the larger number of stars, which include six clearly formed by initial disintegration into 
two ground state *Be nuclei, i.e. by the reaction *O-+hv > 2*Be. These stars are readily 


Photo-disintegration ®O-—>2 §Be. 


recognized by the characteristic ‘ V’ tracks of the alpha-particles finally produced (Wilkins 
and Goward 1950), identification being confirmed by analysis; a typical example is shown 
in the photomicrograph. 

It is, of course, possible for such a track formation to arise by chance in some other mode 
of disintegration. A rough estimate indicates that, after analysis, only about 1 star in 1,000 
would be misinterpreted in this way in the present experiments. 


Atomic Energy Research Establishment, F. K. Gowarp. 


Harwell, Didcot, Berks. 
6th October 1950. J. J. WILKINs. 


Gowarb, F. K., and Wixxins. J. J., 1950, Proc. Ph 
5 5 : Jills ' . Phys. Soc. A, 63, : 
Wikins, J. J., and Gowarp, K. F., 1950, Proc. Phys. Soc. A, 63, se 
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Determination of Fission and Neutron Yields, and the Average 
Neutron Energy in the Photo-disintegration of Uranium 


Although the fission and neutron yields in the photo-disintegration of uranium have been 
measured (Baldwin and Klaiber 1947, Price and Kerst 1950, Baldwin and Elder 1950), the 
interpretation of the results is obscured by two difficulties ; (a) the gamma-ray intensities 
were determined by ionization chambers of somewhat arbitrary design which were not 
intercalibrated or related adequately to numbers of incident quanta, (6) the energy of the 
neutrons emitted in the photo-disintegration was not determined and no estimate could 
therefore be made of the likely errors in the neutron yield measurements. 

To avoid these difficulties the yields have been remeasured (a) using an ionization cham- 
ber, with aluminium walls of variable thickness, whose calibration has been studied both 
practically (Lawson 1950) and theoretically (Flowers, Fossey and Roberts 1950), and (b) by 
determining approximately the average energy of the neutrons. In these remeasurements 
23 Mev. bremsstrahlung from a synchrotron (Fry, Gallop, Goward and Dain 1948) was used. 

The fission yield was determined by two ion chambers placed back-to-back, one of 
which contained a uranium film; this system reduced the size of gamma-ray ionization pulses. 
Only electrons were collected, and, in addition, the chambers were shallow and filled with 
hydrogen so that only the first portions of the fission-fragment and alpha-particle tracks 
were recorded (Pontecorvo and West, private communication); these expedients increased 
the discrimination of fission-particle pulses against alpha-particle pulses and microphony. 
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Measured ‘C’ neutron distributions (shown by dots and triangles) compared with theoretical results 
from a two-group diffusion theory (full curves): (a) for Na—y—D neutrons, (b) for Ra—«—RBe 
neutrons, (c) for photoneutrons from uranium. 

The theoretical curves are calculated for average energies of 0:25, 4:0 and 1:8 Mev. 
respectively. Air-water interfaces are indicated on the graphs by dotted vertical lines. 


The neutron yield was determined by irradiating about 300 gm. of uranium and slowing 
down the neutrons in a large distrene tank (68 cm. x 47 cm. x 44 cm.) containing MnSO, 
in aqueous solution. The activity induced in the stirred solution was compared with that 
induced by a standard Ra—a—Be neutron source. 

The average neutron energy was measured with the same tank, but the distribution of 
C neutrons with distance from the uranium target was investigated by activating rhodium 
foils. ‘The uranium target was separated from the solution in the tank by an air space, so 
that a rigorous analysis of the slowing down process was not possible. A simple ‘ two-group’ 
diffusion theory did, however, give reasonable agreement between experiment and theory 
for Ra-a—Be and Na-y-D sources (of known neutron energy) placed in the same position 
as the target, so this theory was used to determine the energy of the photoneutrons from 
uranium by interpolation. The degree of fit between experiment and theory is shown 
in the Figure. 

The results obtained in the experiments were : 

(1) 1 mg. of uranium gives 24:2-+-2 fissions/réntgen. 
(2) 1 mg. of uranium gives 255425 neutrons/réntgen. 


23 , 5 
(3) i odE=1-5-40-3 mev. barn for the fission yield. 


(4) ee odE=11:5+3 mev. barn for the neutron yield. 


(5) The mean energy of the neutrons emitted is 1°8--0°5 Mev. 
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In results (1) and (2) the réntgen is defined as the value obtained by extrapolating to 
zero wall thickness the straight portion of a curve for log (ionization) against aluminium 
wall thickness. Result (2) is high, compared with previous results (Price and Kerst 1950, 
Baldwin and Elder 1950), by a factor of at least 1-2, and the ratio of neutron to fission yield, 
10-5-42, from results (1) and (2), is distinctly at variance with the ratio 3-4 determined 
by less direct methods (Price and Kerst 1950). Result (3) is approximately 4:5 times that 
of the previous determination (Baldwin and Klaiber 1947). The energy of photo- 
disintegration neutrons has not been measured previously, so that result (5) is important 
not only in relation to the yield determination but also in elucidating the mechanism of 
photo-disintegration. ‘The mean energy is approximately that expected from a normal 
compound nucleus theory. 

Details of the experiments on which this letter is based will be published. 


Atomic Energy Research Establishment, F. K. Gowarp. 
Harwell, Didcot, Berks. E. J. JONEs. 
8th September 1950. H. H. H. Watson. 
D. J. Less. 
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The Ionization of Cosmic-Ray Particles 


It has been observed in photographic emulsions that all particles of momentum higher 
than that of minimum ionization produce tracks of about the same density. Further, some 
tracks of twice minimum ionization have been reported (Bradt and Peters 1950). In some 
cases it has been possible to separate such tracks into two, further along their length, the 
components each having minimum ionization. ‘These observations contradict the original 
ionization theory of Bethe and Bloch, as well as the modifications of Fermi (1940) and of 
Halpern and Hall (1948) which predict a rise of 50% for mesons of 10! ev. It has been 
suggested that this rise might merge into the background of the emulsion because it is due to 
the relativistic broadening of the track. In view of the interest in this question, we present 
here some preliminary results obtained with a proportional counter, which measures the 
whole of the broadened track. 

The proportional counter, filled to a pressure of 86 cm. Hg with a mixture of argon and 
6°% ethylene, was inserted in the cosmic-ray spectrometer briefly described elsewhere (Caro, 
Parry and Rathgeber 1950), between the Geiger counter tray immediately above the magnet 
and the magnet itself. In a run lasting 14 days we measured the ionization of some 400 
particles capable of penetrating 10cm. of lead and having momenta greater than 
2:4 108 ev/c. Our results averaged over a series of momentum ranges are shown in the 
Figure, together with those of Corson and Brode (1938) and Sen Gupta (1943), both of 
which were obtained by counting droplets in delayed cloud chamber photographs. Also 
shown, in plain curves, is the total energy loss in argon for electrons, 4-mesons and protons 
according to Halpern and Hall. Each of the two measurements on electrons has been 
multiplied by a factor which brings the points near minimum ionization into coincidence with 
the curve. Our measurements have been brought to the same level as those of Sen Gupta 
on mesons, which have been multiplied by the same factor as his electron measurements. 

It is evident that whilst the measurements on electrons agree with the theory of Halpern 
and Hall, those on heavy cosmic-ray particles do not. The agreement of the cloud chamber 
and proportional counter results for heavy particles shows that this discrepancy cannot be 
due to experimental errors. Moreover, Sen Gupta’s results for electrons and heavy 
particles were obtained during the same experiment. The following speculations might be 
made to account for this difference between theory and experiment. We shall consider 
the following suggestions : (i) a breakdown of the laws of electrodynamics which would 


Letters to the Editor Q7 


make the ionization constant at high momenta; (ii) that the correction of the Bloch—Bethe 
theory for polarization has been underestimated; (iii) an increasing difference between 
probable ionization and total energy loss with increasing momentum; (iv) that the proportion 
of cosmic-ray particles heavier than p- or 7-mesons increases with increasing momentum; 
(v) that in the ionization process the ju-meson does not behave like a simple charged particle. 

The first two of these effects depend only on the velocity of the particle. As the theory 
gives correct results for electrons it might be expected to apply also to mesons of the same 
velocity. The third factor was examined by taking, in the theoretical relations, 1,000 ev. 
as the maximum transferable energy to an electron. This value corresponds to a cluster of 
30 ion pairs. The fit for electrons appears to be better than for total ionization, but for 
mesons the difference remains, especially as the maximum energy transfer recorded by the 
proportional counter is estimated at 10° ev. 

As to the fourth point, we know that there is little difference between the ionization of 
4- and 7-mesons and so an admixture of 7-mesons would not become apparent in the results. 
It is possible to bring the theoretical ionization into agreement with the results at 2 x 10° ev/c. 
by assuming that all particles of this energy are protons. However even this extreme 
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me? 10° 10’ ev/c 10° 10° 10° 10” 
Specific ionization of cosmic-ray particles as a function of momentum. 
V Electrons (Corson and Brode 1938). @ Mesons (Sen Gupta 1943). 
© Electrons (Sen Gupta 1943). @ Mesons (see present note). 


—___—__—_—— total energy loss in argon. on 
— — — — probable ionization in argon normalized at minimum for electrons. 
<= = aaah = probable ionization in photographic emulsions. 


assumption does not fit the experimental results over the whole momentum range. In 
addition it is known that at least up to 101° ev/c. the majority of particles do not have the 
properties expected of 7-mesons or protons, and that practically all particles in the range are 
p-mesons. We have thus to conclude that p.-mesons differ from electrons in some respect 
other than mass. 

To compare the above conclusions with the grain density measurements in. photographic 
emulsions we have calculated the probable ionization in silver bromide. The results show 
that over the momentum range of the double tracks (Carlson, Hooper and King 1950), 
i.e. above 10’ ev/c., the probable ionization changes by less than 10°%, which agrees with 
observations. It is also evident that, on the basis of our calculations, the difference between 
electrons and mesons observed in argon would be more difficult to detect in photographic 


emulsions. 


Physics Department, P. GoopMaNn. 
University of Melbourne. K. P. NICHOLSON. 
25th September 1950. H. D. RATHGEBER. 
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A New Ultra-Violet Band-System of SiF 


In the course of a search for the spectrum of AgF, two double double-headed bands, 
degraded to shorter wavelengths, were observed with heads at 2109 and 2102 a. which do 
not appear to have been recorded hitherto. Bands already ascribed to SiF (Johnson and 
Jenkins 1927, Asundi and Samuel 1936, Eyster 1937) were strongly developed on the same 
plates, and since the separation of the new bands corresponds closely with the ground-state 
interval *IT/.—°IT,/, in SiF, it seemed likely that they formed part of a new system of this 
molecule. 


Table 1. Deslandres Scheme for Q, and P, Heads 
(The sensitivity varies so steeply with wavelength in this region that estimates of 


intensity have little meaning, but the 0,0 band appears much the strongest, and 
the 1,0 band is probably stronger than the 0,1 band.) 


48952-4 
4 158°1 
48794 °3 
48823-1 
3 161°5 
48661 -6 
48689-8 47025 -4 
2 158i 158-2 
48531-7 46867 -2 
976°7 979°2 
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NE ae re 
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+ Ge) 
co (ee) 
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19.0) 
co [ee 
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i 47408 9 By 46561-7 
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Table 2 
New system a-system 
AG“o4 845-6 cm-t 845-3 cm= 
AG",,» 837-1 836-0 
INGIe 829-5 827-9 
Doublet interval! 160-0 160°8 


Confirmation of this view was obtained when the system was found to be excited in a 
high-frequency electrodeless discharge through SiF, gas prepared by heating BaSiF, to 
about 400° c. in a small furnace attached to the discharge tube. Spectrograms were taken 
with Hilger Medium and Small quartz instruments on Iiford Q.1 plates. "The measurements 
of the inner heads, probably P, and Q,, are summarized in the Deslandres scheme given 
in Table 1. In Table 2, the vibrational intervals in the lower state and the mean doublet 
separation, from Table 1, are compared with the values for the ground-state of SiF taken 
from Eyster (1937). The agreement is good and confirms the identification of the system, 


Letters to the Editor 99 


The constants of the new upper state derived are : 
Ve= 47495 cm! ; G=1005-5(v+ 4) —5:2,(v-+4)?. 


Since the doublet separation is so close to the *IJ, p— Il), separation, the upper state must 
be near to Hund’s case (6). We can therefore proceed to estimate B’ from the head—head 
separations. For both °P;,.—P, and P,—Q, separations, neglecting spin-splitting, A-type 
doubling, and the influence of D terms, we have A= 2B’ Béy/( Boge — B’), where Bog indicates 
the appropriate effective value for *II,). or *II3/.. With Eyster’s value Bj=0-5795, 
Bog=B’(1+B’/A), and an assumed value of «’, equal to 0-004, obtained by Bomosricon 
with similar molecules, the values of B’ given in Table 3 were obtained. Simple 
averaging gives B,=0:621—0-00, (v+-4). 


Table 3. Head-separations and Values of B’(cm™) 


S 
9 


v’, v” PP» B’ Q,-P, B’ ae 27d B’ QP, B’ 
0, 0 17-1 0-624 15-4 0-624 iPS — — 8-9 (0-647) 
Ol 16-1 0-622 15-5 0-619 Dit 20:5 0-612 15-7 0-618 
0572 — = 15-4 0-615 ops = = 12-2 0-622 
1, 0 19-2 0-619 16-9 0-620 Ef 15-0 0-620 19-0 0-605 
et 17-9 0-617 17-1 0-614 4, 3 = —~ 14-1 0-614 
a 15-1 0-620 14-7 0-617 


Violet-degraded bands with heads at 48325-0, 48311-3, 48163-5, 48149-6, 47799-1, 
47784-9 and 47639-0 cm7 were also observed. They do not appear to fit into the present 
system : they may belong to a further system of the same molecule. 


Physical Chemistry Laboratory, W. H. Dove Lt. 
University of Oxford. R. F. Barrow. 
30th September 1950. 


AsunpI, R. K., and SAMUEL, R., 1936, Proc. Indian Acad. Sci. A, 3, 346. 
Eyster, E. H., 1937, Phys. Rev., 51, 1078. 
Jounson, R. C., and Jenkins, H. G., 1927, Proc. Roy. Soc. A, 116, 327. 


Caiculation of the Oscillator Strengths for certain 
Band-Systems of N, and C, 


In this Letter the oscillator strengths (f-values) for the molecular electronic transitions 
associated with the second Positive system of nitrogen and the Swan and Deslandres bands 
of carbon are calculated using the hydrogen-ion approximation. 

The constants of the transitions are given below in Table 1, where R represents the mean 
internuclear distance for the two states in Angstrém units and ». the wavelength appropriate 
to the electronic transition in the same units. 


Table 1 
‘Transition Ae R(A.) R(ao) 
C, 3ITy—Il¢ 5180 1-29 2-44 
Le 1Ty—1Te 3853 1-29 2-44 
N2 3TIu—Te 3371 1-18 2-23 


Expressing these transitions in the form of molecular orbitals, it is immediately seen 
(Mulliken 1932) that they reduce to the form (oy2s)—(eg2p). ‘The wave functions of the 
(oy2s) and (cg2p) molecular orbitals may be approximately represented by a suitable linear 
combination of atomic wave functions as in the case of the hydrogen ion. 

It follows directly that the required wave functions may be written as 
(2,0, 0A|r)—¥2,0,0B)r) ay 

21/244 —S(2, 0, 0)}/? ; 
W(2, 1, 0 A|r)—(2, 1, 0 BIr) Q) 


¥(o,2p|r)= —— mG] — SG, 1,0) }= ee oe 


Y(o42s |r) = 
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where the y’s are normalized atomic wave functions with the principal, azimuthal and 

magnetic quantum numbers n, /, m as indicated, and with the nucleus to which they apply 

denoted by A or B in the usual way. The electronic position vector is denoted by r. 
The overlap integrals are given by 


Sta lan) | Kon Loe Al Oat, eB |e) dy, eee (3) 


the integration being over all space T. 
Defining the dipole moment factor D as 


| ViaeV (x) dr], 


and making use of the expressions for ‘’;(r) and ‘f’s(r) given by equations (1) and (2), it is 
found that 


mt tP(2, 0, 0; ip 16 Oy eZ, 0, 0; 2, i O)}? 
_ {1—S(2, 0, 0)}{1—S(2, 1, 0)} ; 


where PQ;0, 0; 2,1, = [ (2, 0, 0 Alr) rb(2, 1, 0 Ax) dr 
=| (2, 0, 0 B|r)rb(2, 1, 0 B]x) dr, 
and T(2, 0, 0; 2,1, 0)= [ (2, 0,0 Al x) rb(2, 1, 0 BI) dr 


=| (2, 0, 0 B] x) rys(2, 1, 0 Afr) dr. 


The calculation of these integrals was carried out using Slater type atomic wave functions 
with the appropriate effective nuclear charge. In the case of the S and T integrals it was 
found necessary to introduce elliptic coordinates in the integration. Table 2 gives the 
values of the effective nuclear charge for each system together with the values of the integrals 
in atomic units, the unit of length being equal to the radius of the first Bohr orbit. 


Table 2 
Z=3:08e Z=—s3-13e 
C,@IIu—g)  C,@lu—Wg) =N,@Tu—*Ig) 
P(2, 0,0; 2, 1,0) 0-94 0-94 0:77 
T (2, 0,0; 2, 1,0) 0:37 0:37 0-25 
S (2, 0, 0) 0-50 0:50 0-45 
Si (2aeten) —0-30 —0-30 —0-33 


Using equation (4) and the familiar expression f=304D/A for the f-value in terms of 
the dipole moment factor D (in atomic units) where A is the wavelength in a., we finally 
obtain the following results 


Table 3 
‘Transition A (A.) D f-value Life-time (sec.) 
C, @llu—®Ilg) 5180 0-50 0-029 ok Ose 
C, CM —*LHe) 3853 0-50 0-039 Sr <a Oss 
N, @llu—*I¢) 3371 0:37 0-033 SBS CNO-E 


The f-value obtained here for the C, triplet system is in close agreement with the 
experimental value 0-024 obtained by Lyddane and Rogers (1941). 
It is to be noted that in these calculations the hybridization of the orbits has been ignored. 


Department of Physics, G. STEPHENSON. 
Imperial College, 
London, S.W.7. 
20th October 1950. 


MULLIKEN, R. S., 1932, Rev. Mod. Phys., 4, t. 
Lyppang, R. H., and Rocerrs, F. T., 1941, Phys. Rev., 60, 281. 
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REVIEWS OF BOOKS 


Methodes de calcul dans des problémes de mécanique. Pp. 102. (Paris: Centre 
National de la Recherche Scientifique, 1949.) No price. 


In the spring of 1948, two international conferences on methods of calculation in 
problems of mechanics were organized at Marseilles and Paris by the C.N.R.S.; the 
papers contributed to these conferences have been collated and published in this volume. 
The contributors and topics are as follows: J. Valensi (Marseilles), ‘‘The role of applied 
mathematics in engineering’”’; D. N. de G. Allen (London), “‘ Relaxation methods and 
problems of frameworks’’, ‘“‘ Relaxation methods and the solution of differential equations’’, 
together with ‘“‘ Supplementary notes on the application of relaxation methods’; Th. Vogel 
(Marseilles), “‘ The escalator method for the calculation of eigenfrequencies ’’?; M. Picone 
(Rome), ““ New points of view in harmonic analysis’; L. Couffignal (Paris), ‘“‘ Calculating 
machines for harmonic analysis according to the method of H. and Y. Labrouste ”’; 
J. M. Burgers (Delft), ‘‘ Problems related to the theory of turbulence ”’; L. Malavard (Paris), 
“Some recent applications of the method of electrical analogies’; A. van Vijngaarden 
(Amsterdam), “‘ Potential flow about a solid of revolution’’?; F. H. van den Dungen 
(Brussels), “‘ The application of the calculus of variations in fluid mechanics ’’; L. Couffignal, 
“The role of numerical calculation in scientific and technical research ’’. 

The papers naturally vary in style and standard but they are, without exception, lucidly 
written and the book can be thoroughly recommended to a physicist or engineer who wishes 
to acquire some idea of modern techniques used in the realm of classical mechanics. R. M. D. 


Vorstufe zur theoretischen Physik, by R. Becker. Pp. vii+172. (Berlin: 
Springer-Verlag, 1950.) 7.50 DM. 


‘This book is interesting from a didactic point of view. Each section can be described 
as a sort of verbatim rendering of a (very good) lecture for beginners. Every detail is 
carefully written down, including observations of a methodical character, paternal advice, 
disclosure of useful tricks, and even humorous illustrations with which every teacher 
intersperses his lectures, but which he usually leaves out when he writes a textbook. ‘he 
subjects treated are restricted, apart from some slight digressions, to the field of mechanics 
and thermal phenomena (including kinetic theory). Within this range the selection is the 
fairly obvious one which, I think, every one of us would make in order to illustrate the 
application of the most fundamental methods of classical theory. The treatment seems 
to me somewhat uneven: some sections are extremely elementary, while neighbouring 
ones are on a decidedly higher level. The discussion of some questions of principle 
(e.g. the fundamental equation of dynamics) is, in my opinion, too cursory and even 
slovenly; that of the thermoelectric phenomena is definitely in error. Moreover, there 
is no reference anywhere to the historical development of the subjects; this is admittedly 
better than the uncritical repetition of historical mis-statements in which too many authors 
indulge, but just for pedagogical purposes a judicious use of historical material would be 
of the greatest help, and would materially enliven the picture. Becker’s book will certainly 
convey to the student a sound understanding of the attitude of mind of the theoretical 


physicist, but it is a bit too dry to give him even a hint of the beauty of his quest. 
L. ROSENFELD. 
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THE INSTITUTION OF ELECTRICAL ENGINEERS 
Proof and Reprint Service 


The following details of the Proof and Reprint Service of 'The Institution of Electrical 
Engineers, London, have been received, with a request that they should be brought to the 
attention of our members. Formerly the proofs of a paper allocated for reading before the 
Institution of Electrical Engineers were made available to members only 10 days before the 
meeting; thus, if accepted at a time when the current programme of meetings was complete, 
authors’ manuscripts might sometimes lie dormant for several months. Since October 
1948 ‘early proofs ”’ of all papers have been printed without avoidable delay. 

As soon as the early proofs of a paper are available, an announcement to that effect, 
together with a synopsis, appears in the monthly Yournal of the Institution. Application 
for such proofs can then be made, and those who avail themselves of this service are given 
the opportunity of applying for a reprint, from the Proceedings of the Institution, of the 
paper in its final form, together with a report of the discussion. The object of the reprint 
service is to ensure that the permanent record available to its users shall not consist of an 
early proof which, however carefully prepared, may contain minor errors and will certainly 
not include the discussion. 

The inclusive charge for an early proof and its associated reprint is 2s. 6d. (post free). 
This charge is made whether or not the reprint is required, but a reprint alone can be 
supplied at the price of 1s. 6d. (post free) if ordered shortly after the paper has been read in 
London. 

Requests for proofs should be sent, accompanied by a remittance, to The Secretary, 
The Institution of Electrical Engineers, Savoy Place, London W.C.2, England. 
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ABSTRACTS FOR SECTION B 


Some Further Investigations of Ionospheric Cross-Modulation, by 1. J. Suaw. 


ABSTRACT. ‘The work in this paper forms the continuation and development of the 
experiments described in a paper by Ratcliffe and Shaw published in 1948. Measurements 
of the phase and amplitude of the modulation transferred to an unmodulated wave by a 
modulated disturbing wave in the ionosphere are used to provide information about that 
part of the ionosphere in which the cross-modulation occurs. 

By measuring almost simultaneously the height of reflection of the wanted wave and the 
characteristics of the transferred modulation, a value of 1-4 x 10° per second for the collision 
frequency of electrons with molecules at a height of 92 kilometres has been deduced. This 
measurement has been made possible by the use of a very sensitive automatic apparatus 
which is briefly described. 

Special attention has been paid to the ‘ gyro-resonance effect ’ in cross-modulation first 
pointed out in a theoretical paper by Bailey in 1937. The experiments appear to show that 
the predicted enhancement of the cross-modulation when the disturbing station radiates 
on the local gyro-magnetic frequency does not exist. 

A modification of the original theory is given which is of great use in the calculation of the 
cross-modulation to be expected for a given pair of transmitters. This modified theory is 
used in the discussion of the gyro-resonance effect and of the observed decrease of the cross- 
modulation at dawn. ‘The agreement between the calculated and observed magnitudes of 
cross-modulation is stressed in view of the fact that the calculations involve a knowledge of 
the absorption of the wanted wave incident obliquely upon the ionosphere. 


Some Meteorological Aspects of Radio Duct Formation, by E. KNIGHTING. 


ABSTRACT. Three problems of meteorological interest concerning the formation of radio 
ducts are treated. The first concerns the rate of growth of ducts, with special reference to 
ducts over the sea; it is shown that ducts form quickly at first and then more slowly as the 
maximum duct width is approached. The second problem concerns the formation of ducts 
under conditions of nocturnal cooling; it is shown that there is a limit to the height which the 
duct may attain. The third problem concerns the width of a radio duct necessary to enclose 
the track width of a mode associated with a propagated signal of wavelength A assuming a 
power law of modified refractive index with height; it is shown that the necessary width is 
proportional to \?, and the factor of proportionality, which depends upon the index in the 
assumed power law, is calculated. 

The three problems have a common meteorological basis, for each is concerned with 
turbulent motion near the earth’s surface. The law of turbulent diffusion assumed here is 
that the coefficient of eddy diffusion is proportional to a power of the height, following the 
theory developed by O. G. Sutton and others. 


The Properties of Spark Counters of the Rosenblum Type, by R. D. CoNNor. 


ABSTRACT. The design and operation of Rosenblum-type spark counters having 
plane, concave and convex cathodes is described, and discussed with reference to the 
counting characteristics, useful life and background counting rate. ‘The absolute efficiency 
of the counter in air has been determined for a-particles, knock-on protons, and /-particles, 
giving values of about 90%, 5% and 0:0002%, respectively. A marked dependence of 
efficiency on specific ionization and air pressure is evident. The counter is directional in 
all its forms. The effect of employing gases other than air is also investigated and a 
qualitative description of the counting characteristics is given. Alpha-particle ranges may 
be determined using either a single counter or two counters used in anticoincidence as is 
described in the text, 
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A Coawial-Cylinder Elastoviscometer, by J. G. OLDRoyD, D. J. STRAWBRIDGE and 
B. A. Toms. 


ABSTRACT. An instrument is described which was designed to measure the elastic and 
viscous properties of liquids, primarily of solutions of linear high polymers whose elastic 
properties had already been observed qualitatively. The liquid is contained in a narrow 
annulus between two long vertical cylinders. The outer cylinder can be driven in steady 
rotation at a known speed, or in harmonic angular oscillations of known amplitude and 
frequency. In steady rotation, the torque on the stationary inner cylinder is measured at 
different speeds, and interpreted to give a viscosity coefficient as a function of shear stress. 
In oscillatory motion, the amplitude of oscillation of the inner cylinder, when constrained by 
a torsion wire, is measured and plotted against the frequency. The shape of the resulting 
graph can be interpreted so as to distinguish and measure different types of elastic behaviour. 
The annular gap, the restoring constant of the torsion wire and the amplitude of oscillation 
of the outer cylinder can be varied in discrete steps. The angular velocity of the outer 
cylinder can be varied continuously from 0-3 to 3,000 rev/min. in steady rotation, and the 
frequency can be varied continuously from 0:25 to 25 c/s. in steady oscillation. The 
performance of the instrument in both steady-state and oscillatory experiments is assessed 
by comparing observed and predicted results for Newtonian liquids of known viscosity. 
An example is given of the behaviour of a typical polymer solution in both kinds of experi- 
ment, together with a suggested quantitative interpretation of the results. 


The Evaporation of Titanium, by L. G. CARPENTER and W. N. Maire. 


ABSTRACT. Data are presented on the evaporation rate of 8 titanium over the tem- 
perature range 1,650-1,810° k. 
These are represented by 


e242 
log 1M = ae -+9-79, 


where E is in gm.cm~? sec~ and T is in degrees Kelvin. 

The latent heat of vaporization is deduced both from the measured values of evaporation 
rate, together with known thermodynamic functions of the solid and gas, and also from the 
temperature coefficient of evaporation rate. 

An estimate is made of the variation of total radiation with temperature, and of the vapour 
pressure at the melting point. 

Experimental details are given. 


An Electrostatic Focusing System and its Application to a Fine Focus X-Ray 
Tube, by W. EHRENBERG and W. E. Spear. 


ABSTRACT. ‘The description is given of an electrostatic focusing system by which 
electrons emitted from the tip of a hairpin filament are reunited at a distance of 2 to 3 cm. 
No intermediate potentials are required. An x-ray method is used for the experimental 
investigation of its properties. The electron spot is a reduced image of the tip of the 
filament and has a diameter of about 40 «. On the basis of this system, a fine focus x-ray 
tube is designed, with a peak specific target loading of about 11 kw/mm2. Some pe1formance 
data are given. 


Low Reflectance Coatings « .d the Sensitivity Curve of the Eye, by J. H. vAN DER 
Merwe and H. VERLEGER. 


ABSTRACT.  Heterochromatic light is treated and the spectral sensitivity of the human 
eye is taken into account. In spite of certain simplifying assumptions, conclusions of a 
general character may be drawn, the most important being that the optical properties 
required for low reflectance coatings can be obtained by treating the incident light as 
monochromatic, 
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Shower C Shower D 


(See back of plate for description.) 


DESCRIPTION OF PLATE 


Showers A, B, C and D were observed under 25 cm. of paraffin wax using a magnetic 
field of 7,500 gauss. 


Shower A. 


This photograph shows at least two related nuclear interactions in the lead plate induced 
by a neutral primary. ‘T'rack 1 is associated with the main interaction; the momentum is 
of the order ot 10° ev/c. and the track appears to emerge at an angle of 10° frorn the core of 
the main event. The track appears to start in the core and not in the plate to the left of the 
core. Thus track 1 may have been deflected from the core through an angle of 10° and the 
process may be similar to that described by Rochester and Butler (1947), that is, the spon- 
taneous decay of a meson with a mass between that of the 7-meson and the proton. 


Shower B. 


A penetrating particle (track 1) is accompanied by a neutral-induced secondary inter- 
action. This event is unusual since two distinct electron cascades were formed; the angle 
between the cores of these cascades is about 18°. The mean directions are shown by 
arrows 2 and 3. Unfortunately shower particles are mixed with the electrons, and therefore 
it is impossible to determine the total energy of each cascade. The cascades may have been 
produced by the decay of a neutral 7-meson. 


Shower C. 


This event is a simple shower of three collimated penetrating particles, tracks 1, 2 and 3. 
Track 1 is scattered through 31° in the plate. It is too short for measurement above the 
plate, but below it is negative with a momentum of 10° ev/c. 


Shower D. 


Tracks 1 and 2 are almost parallel and appear to produce electron cascades of very 
different energies in the lead plate. ‘Track 1 has amomentum which is too high for measure- 
ment and, in addition to the large cascade, produces track 3, which is negative with a 
momentum of 8x 10% ev/c. This particle was projected at an angle of 40° to the direction 
of track 1 and so it is very unlikely to be an electron and must be a meson. ‘Tracks 1 and 3 
intersect in the top half of the lead plate and, if the cascade was formed in this region, the 
total energy in the electronic component is likely to be about 101° ev. 

Track 2 is positive with a momentum of 1:5 x 10° ev/c. and is likely to be due to an 
electron which produced the cascade of three electrons below the plate. 

Tracks 1 and 2 are unlikely to be part of an air shower. In order to produce the large 
cascade, track 1 must have an energy greater than 101° ev., and such an electron would be 
very close to the core of an air shower. 

If the right-hand event is a nuclear interaction, then it is of a very unusual type, and in all 
probability there are mesons and protons, other than track 3, which must be located in the 
core of the cascade and so be unobserved. If this is so, the whole event must be of exception- 
ally high energy. 
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